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INTRODUCTION 

The  General  Aviation  Dynamics  (GAD)  model,  Is  implemented  in 
NUCLEUS*,  a  computer  software  system  developed  at  Battelle.  GAD  cur¬ 
rently  resides  on  both  the  Battelle  computer  and  the  United  Computing 
Systems  (UCS)  computer;  it  can  be  accessed  on  either  machine  via  tel¬ 
ephone  from  anywhere  in  the  U.S.,  using  either  a  remote  batch  terminal 
or  an  interactive  terminal.  This  GAD  systems  manual  contains  a  complete 
description  of  the  computer  program.  Appendix  A  of  this  manual  consists 
of  a  complete  listing  of  the  nucleus  GAD  program,  Appendix  B 


*NUCLEUS:  Numerical  CLassif ication  and  Evaluation  jSystem. 
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contains  a  full  batch  run  of  the  model.  Appendix  C  is  a  user  handbook 
for  the  Nucleus  language,  and  Appendix  D  is  a  user's  guide  for  running 
the  GAD  model  interactively. 

PROGRAM  DESCRIPTION 


NUCLEUS  is  a  Fortran  based  computer  methodology  featuring:  a 
structured  programming  language,  dynamic  simulation  modeling,  data  base 
management,  report  generation,  conversational  programs  and  batch  or  on¬ 
line  access.  The  GAD  model  is  a  dynamic  simulation  model  written  in 
NUCLEUS  and  with  a  conversational  dialogue  superimposed  on  it.  This 
conversational  dialogue  simplifies  user  interaction  with  the  system. 
The  user  is  led  through  the  steps  of  the  model,  one  step  at  a  time, 
selecting  input  and  output  options  by  answering  self-explanatory 
questions  written  in  English.  Thus,  even  a  user  unfamiliar  with  com¬ 
puters  should  experience  no  difficulty  in  using  GAD.  The  user-system 
interaction  introduced  by  the  conversational  dialogue  also  allows  direct 
examination  and  modification  of  the  input  data  and  assumptions  during 
on-line  use  of  the  system.  Thus,  the  effect  on  the  general  aviation 
forecasts  of  different  scenarios  may  be  investigated  by  the  user  during 
a  single  on-line  session  with  the  model. 

In  more  detail,  the  GAD  system  is  a  modularized  simulation 
system  with  a  hierarchical  or  tree  structure  imposed  on  the  individual 
segments  of  the  program.  The  need  for  modularizaton  or  segmentation 
arises  from  the  fact  that,  in  its  entirety,  the  system  occupies  too  much 
computer  memory  to  permit  interactive  use.  With  segmentation,  only 
portions  or  segments  of  the  program  are  resident  in  core  at  any  time. 
The  entire  program,  in  the  form  produced  by  the  NUCLEUS  interpreter,  re¬ 
sides  on  a  file  from  which  the  different  segments  are  loaded  into  core 
as  required.  The  segmentation  diagram  in  Figure  1  shows  the  tree  struc¬ 
ture  of  the  segments.  A  description  of  the  function  of  each  segment  is 
given  in  Table  1 .  There  are  ten  linked  segments  organized  into  three 
levels  of  hierarchy  or  dominance.  The  linkages  between  the  segments  in 
the  diagram  reflect  decreasing  dominance  in  a  top-to-bottom  direction. 
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1ST 


FIGURE  1.  THE  SEGMENT  DIAGRAM  OK  THE  GAD  SYSTEM 
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TABLE  1.  SEGMENT  STRUCTURE  OF  THE  GAD  SYSTEM 


Segment 

Level 

Description 

GAD 

0 

Controls  the  entire  GAD  system. 

INITIAL 

1 

Allows  the  user  to  select  the  ending 
year  for  the  simulation,  and  to  choose 
between  running  the  simulation  with 
the  initial  assumptions  unchanged  or 
changing  the  initial  assumptions 
of  the  forecast. 

START 

1 

Writes  the  initial  assumptions  for 
data  out  to  the  data  base. 

CHNGE 

1 

Allows  the  user  to  change  the  initial 
assumptions  for  the  exogenous  data. 

DYN 

I 

The  dynamic  simulation  segment  of  the 
model. 

SHOW 

I 

Allows  the  user  to  select  output 
tables  to  be  printed. 

SHWPL 

1 

Allows  the  user  to  select  plots  of 
the  results  of  the  forecast  to  be 
displayed. 

SNTVT 

1 

Controls  the  display  of  sensitivity 
results  for  comparing  one  simulation 
run  with  another  baseline  simulation. 

SHWSA 

2 

Allows  the  user  to  select  sensitivity 
tables  to  be  printed. 

SHPLS 

2 

Allows  the  user  to  select  sensitivity 
plots  to  be  displayed. 

I 


I 
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For  example  the  master  control  segment  GAD  (at  the  top  of  the  diagram) 
has  dominance  over  all  other  segments  of  the  system.  Information, 
available  or  computed,  in  the  GAD  segment  is  available  to  all  the  other 
segments.  Segment  GAD  is  resident  in  core  for  the  entire  period  of  user 
interaction  with  the  GAD  system.  Segments  SHWSA  and  SHPLS  at  level  2 
are  dominated  by  SNTVT  at  level  1  and  GAD  at  level  0.  Only  one  of  SHWSA 
and  SHPLS  may  be  in  core  at  any  one  time;  if  either  one  is  in  core,  then 
SNTVT  and  GAD  will  also  be  resident  in  core. 

To  explain  further  the  concept  of  segment  hierarchy,  Figure  2 
shows  a  primitive  segment  diagram.  In  the  context  of  this  diagram  the 
following  general  statements  apply: 


Level 


FIGURE  2.  A  PRIMITIVE  SEGMENT  DIAGRAM 
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•  A  segment  may  or  may  not  be  dominated  by  another 
segment.  For  example,  segments  B  and  C  are 
dominated  by  A,  but  B  does  not  dominate  C,  and  C 
does  not  dominate  B. 

•  If  segment  B  is  dominated  by  A,  then  segment  A 
cannot  be  dominated  by  B. 

•  If  segment  A  dominates  segment  B,  then  all 
information  in  segment  A  is  available  to 
segment  B. 

•  The  information  in  segment  B  is  not  available  to 
segment  C,  and  vice  versa,  except  for  any  overlap 
through  segment  A.  There  is  no  information  flow 
between  segments  B  and  C,  except  through  the 
dominant  segment  A. 

•  From  the  user  standpoint,  the  information  in 
segment  B  is  available  only  through  segment 
A.  Similarly,  the  information  in  segment  C 
is  available  only  through  segment  A. 

•  Operationally,  segment  dominance  is  manifested  in 

terms  of  the  following  computer  memory  (core) 
configurations:  segment  A  may  exist  in  core 

alone;  segment  B  can  exist  in  core  only  with 

segment  A  and  in  the  order  "A  followed  by  B"; 

segment  C  can  exist  in  core  only  with  segment  A 
and  in  the  order  "A  followed  by  C”;  segments  B  and 
C  cannot  coexist  in  core.  Thus  the  only  possible 
core  configurations  are:  A,  A+B,  and  A+C. 

Appendix  A  contains  the  full  Nucleus  listing  of  the  GAD 

program.  In  this  program,  the  segments  given  above  are  further 
structured  into  models  which  are  groups  of  equations  and  instructions 
that  perform  one  or  more  functions.  Table  2  lists  the  models  contained 
with  the  different  segments  and  gives  a  brief  description  of  their 
functions.  All  the  models  in  the  system,  with  the  exception  of  DYNAM, 
are  static,  performing  functions  that  do  not  vary  with  time.  DYNAM  is  a 


TABLE  2.  MODELS  CONTAINED  IN  THE  SEGMENTS  OF  THE  GAD  SYSTEM 


Segment 

Level 

Description 

GAD 

MAIN 

This  is 
tion  of 

the 

the 

model  that  controls  execu- 
entire  GAD  system. 

INITIAL 

ASKAGAIN 

Allows 

the 

user  to  view  the  different 

steps  involved  in  the  system  and  to 
choose  a  forecast  using  the  initial 
assumptions  unchanged  or  a  forecast 
in  which  the  initial  assumptions  are 
modified  from  the  base  line. 

ENYR  Allows  the  user  to  select  the  ending 

year  for  the  simulation. 

FIRST  Prints  the  heading  at  the  start  of 

the  GAD  run  and,  if  one  simulation  has 
already  been  run,  permits  the  user  the 
option  of  continuing  with  another 
simulation. 

START  BUILD  Writes  the  initial  assumptions  for  the 

exogenous  data  on  the  data  base. 

CHNGE  PRCH  Allows  the  user  to  change  the  initial 

assumptions  for  the  pilot  variables. 

CHVAR  Allows  the  user  to  change  the  values 

of  the  economic  variables  yearly  at  a 
constant  rate,  or  to  enter  them  year 
by  year. 

ECH  Allows  the  user  to  change  the  initial 

assumptions  for  the  economic  data. 

YRCHK  Allows  the  user  to  select  the  year  in 

which  changes  to  the  initial  assumptions 
for  cost  variables  are  to  begin. 

FOOST  Allows  the  user  to  change  the  initial 

assumptions  for  the  fixed  cost. 

VCOST  Allows  the  user  to  change  the  initial 

assumptions  for  the  variable  cost. 
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TABLE  2.  (Continued) 


Segment  Level 


Description 


CHANGE 

Requests  the  user  to  select  the  variables 
prior  to  calling  the  CHANGE  model. 

DYN 

GADYN 

Contains  the  equations  for  the  dynamic 
simulation. 

DYNAM 

Control  model  for  GADYN. 

SHOW 

PRINT 

Allows  the  user  to  select  variables  to 
be  tabulated. 

SHWPL 

SELC 

Requests  the  user  to  select  an  aircraft 
user  category  for  a  plot  of  a  multi- 
dimensioned  variable. 

SELT 

Requests  the  user  to  select  an  aircraft 
type  for  a  plot  of  a  multi-dimensioned 
variable. 

XAXIS 

Allows  the  user  to  select  the  independent 
variable  for  a  plot. 

PLOT 

Allows  the  user  to  select  the  dependent 
variable  for  a  plot. 

SNTVT 

LOADD 

Controls  the  output  display  of  sensiti¬ 
vity  variables. 

SHWSA 

PRINS 

Allows  the  user  to  select  sensitivity 
variables  to  be  tabulated. 

SHPLS 

SELC 

Requests  the  user  to  select  an  aircraft 
user  category  for  a  sensitivity  plot  of 
a  multi-dimensioned  variable. 

SELT 

Requests  the  user  to  select  an  aircraft 
type  for  a  sensitivity  plot  of  a  multi- 
dimensioned  variable. 

PLOTS 

Allows  the  user  to  select  a  sensitivity 

variable  to  be  plotted. 
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dynamic  model  containing  equations  that  describe  the  continuous  eval¬ 
uation  through  time  of  interacting,  time-dependent  variables.  T:  is  the 
model  DYNAM  that  generates  the  simulation  results  through  time. 

Appendix  B  contains  a  sample  batch  run  of  the  GAD  svstem 
including  a  baseline  simulation,  a  simulation  based  on  an  increase  in 
the  federal  fuel  tax  commencing  in  1979,  and  sensitivity  results  com¬ 
paring  the  results  of  the  two  sets  of  forecasts.  The  data  cards  re¬ 
quired  to  run  the  simulations  are  simply  the  responses  to  the  questions 
that  are  printed.  Appendix  C  contains  a  description  of  the  Nucleus 
language  and  details  of  the  commands  available  in  Nucleus.  Appendix  D 
is  a  user's  guide  to  exercising  the  GAD  model  on  either  the  Battel  le  or 
UCS  interactive  computer  systems. 
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A  SAMPLE  BATCH  RUN  01  GAD 
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VOU  ARE  ENTERING  THE  GENERAL  AVIATION  OVNAHICS  HOTEL 
CREATEO  AT  OATTELLE  COLUMBUS  LABORATORIES! 

WRITTEN  IN  THE  MOOELING  LANGUAGE  NUCLEUS. 
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1.0  INTRODUCTION' 

The  NUCLEUS  ays  cam  cons  is cs  of  a  user  language ,  an  interpreter 
which  interprets  and  performs  the  operations  specified  in  the  user  language, 
and  an  interface  Co  FORTRAN.  Hereafter,  unless  other-vise  specified,  the 
word  "system"  will  mean  the  NUCLEUS  system. 

The  user  language  of  the  system  is  a  structured  programming  lan¬ 
guage  according  to  whose  rules  and  conventions  the  user  writes  and  organizes 
his  instructions  for  entry  into  the  system  either  via  keypunched  cards  cr  bv 
typing  entries  on  the  keyboard  of  an  on-line  computer  terminal.  As  a 
structured  programming  language  it  has  the  following  two  characteris tics : 

(i)  it  has  no  statement  labels  and  thus  no  GO  TO  statements,  and  (ii)  major 
system  functions  —  such  as  report  generation,  data  base  access,  secondarv 
analysis,  dialogue  control,  etc.  —  are  accessed  via  single  statements.  The 
complete  set  of  user  instructions  specifying  a  job  to  be  done  by  the  system 
is  called  a  source  program.  The  component  instructions  of  a  source  program 
are  called  statements .  The  user  language  is  especially  useful  for  writing 
programs  which  construct,  manipulate,  and  display  arrays  of  numeric  data. 

The  interpreter  of  the  system  is  a  set  of  ANS*  (American  National 
Standard)  FORTRAN  programs  which  interpret  the  statements  of  a  source 
program  and  generate  interconnected  cables  of  internal  instructions  which 
later  direct  the  execution  of  these  statements. 


*  The  single  exception  to  this  involves  the  direct  access  routines  which 
though  they  are  not  in  ANS  FORTRAN  have  equivalents  in  all  major  computers. 


The  interface  to  FORTRAN  greatly  expands  the  potential  capabil¬ 
ities  of  the  user.  The  user  may  use  the  subroutines  of  the  system  in  other 
independently  written  FORTRAN  programs.  This  is  especially  important  in 
so  far  as  the  data  base  capabilities  of  the  system  are  concerned.  The 
call  to  the  FORTRAN  subroutine  which  manipulates  the  data  base,  is  des¬ 
cribed  in  Chapter  5  of  this  section. 

The  machine  independence  of  the  NUCLEUS  system  is  a  very  important 
feature:  because  the  interpreter  of  the  system  is  machine  independent  (i.e., 
written  in  ANS  FORTRAN) ,  any  application  program  written  in  the  user  language 
of  the  system  will  also  be  machine  independent,  provided  the  user-supplied 
FORTRAN  subroutines,  if  any,  that  are  introduced  via  the  FORTRAN  interface 
are  machine  independent. 

An  overview  of  the  User  Lnaguage  is  given  in  Chapter  2.0.  This 
chapter  includes  both  a  discussion  of  the  structural  elements  and  a  descrip¬ 
tion  of  the  general  syntax  rules  of  the  language. 

Chapter  3.0  contains  explicit  descriptions  of  each  statement  of 
the  User  Language.  For  easy  reference,  the  primary  word  of  each  statement 
is  placed  in  capital  letters  on  the  upper  left-hand  corners  of  all  pages 
describing  that  statement. 

Chapter  4.0  describes  the  error  messages  of  the  User  Language, 
and  Chapter  5.0  describes  the  interface  of  the  system  to  FORTRAN. 

Appendix  A  contains  sample  programs  -  batch  and  interactive  - 
written  in  the  User  Language  of  the  system. 


2.0  OVERVIE',;  OF 


The  s tat aments  of  a  source  program  are  the  means  by  which  one  user 
specifies  to  the  system  the  desired  operations  of  a  job  or  application.  A 
statement  may  specify  one  of  three  types  of  operations: 

(1)  Procrar,  control,  i.e.,  it  may  command  the  system  to  perform 
some  overall  program  control  function,  such  as  starting  a 
new  program,  ending  a  program,  opening  a  (fata  rile,  saving 
a  program,  loading  a  previously  saved  program,  or  adding  an 
array  to  a  data  file. 

(2)  Program  structure  and  identifier  definition,  i.e.,  it  may 
define  the  identifier  for  a  structural  element  of  the  pro¬ 
gram  and  allocate  central  memory  for  the  associated  struc¬ 
ture  and  information  for  later  reference  and  manipulation 
within  the  program. 

(3)  Data  manipulation,  i.e.,  it  may  instruct  the  system  to  make 
reference  to  or  manipulate  the  information  associated  with 
previously  defined  identifiers. 

The  basic  commands ,  definitions ,  and  instructions  of  the  user 
language  are  listed  in  Table  1-1.  The  first  word  of  any  user  language  state 
ment  must  be  one  of  the  capitalized  statement  words  in  this  table,  except 
in  the  following  three  situations: 

•  Equation  statements 

•  Model  execution  statements 

•  Table  execution  statements. 

In  an  equation  statement  the  first  word  of  the  statement  must  be  the  identi¬ 
fier  of  a  previously  defined  dataset,  A  model  execution  s tatement  consists 
of  the  identifier  of  a  previously  defined  model  and  causes  execution  of  the 
instructions  in  the  model.  A  table  execution  statement  consists  of  the 
identifier  of  a  previously  defined  table  and  causes  the  printing  of  the  tabl 
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TABLE  1-1.  THE  USER  LANGUAGE  STATEMENTS 


Commands  -  Program  Control 


SYSTEM 

START 

STOP 

OPEN 

SAVE 

LOAD 

TIME 

ADD 


Secs  various  system-control  parameters 

Starts  a  new  program 

Stops  execution;  ends  program 

Opens  a  data  file 

Saves  a  program 

Loads  a  previously  saved  program 

Defines  the  time  parameters  to  control  dynamic  simulation 
Adds  an  array  to  a  data  file 


Definitions  -  Program  Structure 

INDEX  Defines  an  index 

DATA  Defines  a  dataset 

GRAPH  Defines  a  graph 

MODEL  Introduces  a  model 

TABLE  Defines  an  output  table 

END  Ends  the  definition  of  a  model,  an  ASK  statement,  a  con¬ 

ditional  DO  or  IF  statement  or  the  definition  of 
multiple  indexes,  or  datasets,  or  graphs,  or  tables 


Instructions  -  Data  Manipulation 


READ 

SET 

SHOW 

equation 

model 

table 

RATE 

LEVEL 

WRITE 

INPUT 

DROP 

PUT 

IF 

DO 

TELL 

ASK 


ELSE 


Reads  index  descriptors,  or  dataset  values,  or  graph  values 
Sets  Index  size  and  ordering  of  index  entries 
Shows  values  of  a  dataset  in  table  or  plot  form 
Computes  the  values  of  a  dataset  or  quantity 
Executes  a  model 

Displays  a  structured  tabular  report 

Introduces  rate  equations  in  a  dynamic  simulation  model 
Introduces  level  equations  in  a  dynamic  simulation  model 
Writes  values  to  a  table  on  a  direct  access  file 
Inputs  values  from  a  table  on  a  direct  access  file 
Drops  data  set  values  from  working  storage 
Puts  data  set  values  into  working  storage 
Executes  user  instructions  if  some  condition  is  met 
Executes  user  instructions  while  or  until  some  condition  is  met 
Tells  the  user  something,  l.e.,  issues  a  message 
Asks  che  user  a  question;  based  on  the  response,  the  system 
executes  one  and  only  one  of  many  possible  sets  of 
alternative  instructions 

Defines  an  execution  branch  of  instructions  within  an  IF 
or  an  ASK  statement 
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The  basic  structural  ingredients  of  elements  of  the  user  language 
and  the  system  are  listed  in  Table  1-2.  Every  application  or  problem  is  con¬ 
ceptualized  and  formulated  in  terms  of  these  elements  before  being  coded 
into  a  program  of  statements. 


TABLE  1-2.  THE  STRUCTURAL  ELEMENTS  OF 
THE  USER  LANGUAGE 


Identifier  A  name  for  an  index,  a  dataset,  a  model,  a  graph,  a 

table,  or  a  program  module. 

Index  A  classification  (or  subscripting)  scheme  for  the  values 

of  structured  (nonscalar)  datasets.  Indexes  are  sub¬ 
scripts  of  datasets. 

An  n-dimensional  storage  array  containing  numeric  infor¬ 
mation,  where  n=0 , 1 , . . . , 10 .  Datasets  are  the  structural 
components  of  equations,  i.e.,  the  variables  of  a  problem 
or  application. 

An  ordered  sequence  of  instructions  or  equations  that 
may  be  compiled,  i.e.,  defined,  as  a  unit  for  later 
execution  as  a  unit. 

An  ordered  set  of  pairs  of  values  defining  the  x-  and 
y-values  of  the  points  of  a  function,  y*f(x). 

Table  An  ordered  sequence  of  datasets  defined  as  a  unit  for 

later  execution  as  a  titled  tabular  display  of  the  values 
of  these  datasets. 

Program  Module  An  ordered  sec  of  user  language  statements  that  may  be 

'  saved  for  later  loading  and  execution. 

Equation  An  instruction  containing  the  ■  sign  and  setting  the 

value (s)  of  the  dataset  at  the  left-hand  side  of  the 
■  sign  equal  to  the  value (s)  of  the  right-hand  side. 

Data  File  A  file  of  information  stored  in  mass  storage  outside  the 

working  space  (allocated  core)  of  a  program.  Data  files 
may  contain  either  program  modules  or  arrays  of  numeric 
data . 

TMS  Array  A  n-dimensional  storage  array  containing  numeric  infor¬ 

mation  and  located  on  a  TMS  data  file,  n-0,1, . . . ,10.  The 
acronym  TMS  stands  for  Table  Management  System  and  reflects 
the  manner  in  which  the  data  on  the  file  is  brought  into 
working  storage. 


Dataset 


Model 


Graph 


The  Structural 


ener.cs 


of  the  Is e r  Languaae 


The  structural  elements  of  the  user  language  are  listed  in  Table 
1-2.  The  semantics ,  i.e.,  the  meanings,  of  these  elements  are  further 
described  in  this  chapter.  The  svntax  of  the  language,  ie. ,  the  manner  in 
which  these  elements  are  specified  in  the  statements  of  the  user  language, 
is  described  in  Chapter  3.0. 

2.1.1  Identifier 


An  identifier  is  a  user-supplied  name  or  abbreviation  for  some 
structural  component  of  a  program  -  such  as  an  index,  or  a  dataset,  or  a 
model,  or  a  graph,  or  a  table,  or  a  program  module,  or  a  TMS  array.  The 
identifier  of  a  component  i3  used  to  define  and  to  refer  to  that  component 
later  in  the  program.  All  identifiers  except  those  for  TMS  arrays  must  con¬ 
tain  no  more  than  five  alphabetic  and/or  numeric  characters  and  their  first 
character  must  be  alphabetic.  They  must  contain  no  blanks  or  any  other 
special  characters.  To  avoid  ambiguities,  identifiers  must  be  unique  within 
a  given  program.  TMS  arrays  are  identified  via  an  ordered  pair  of  integer 
values.  The  first  integer  in  this  pair  defines  the  particular  data  file 
within  the  data  file  sec  being  used  which  contains  the  array.  Within  this 
version  of  the  system  a  given  data  file  set  may  contain  a  maximum  of  three 
files;  thus,  the  data  file  identifier  must  be  a  number  between  one  and  three. 
The  second  integer  in  the  TMS  array  defines  the  sequence  number  of  the 
desired  array  from  among  the  arrays  contained  on  the  file;  thus,  it  is  an 
integer  value  between  1  and  n.  The  value  of  n  is  unique  for  each  file  and 
is  established  by  the 'user  when  that  file  is  first  defined. 

The  identifier  of  a  structural  element  does  not  represent  a  single 
machine  or  file  address;  rather,  it  refers  to  the  structured  complex  of 
related  information  contained  in  the  element.  When  the  identifier  appears  in 
a  statement,  the  system  automatically  references  the  particular  portion  of 
the  element  that  is  relevant  to  the  context  of  use  of  that  statement. 
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2.1.2  Index 


The  index  is  one  of  the  two  most  basic  concepts  in  the  system;  the 
other  being  the  dataset  concept.  The  index  is  the  organizational  classifi¬ 
cation  scheme  in  terms  of  which  the  internal  structure  of  the  datasets  in  a 
program  is  defined.  Indexes  serve  as  generalized  subscripts  of  datasets. 

For  example,  age  and  sex  may  be  two  indexes  classifying  the  age  and  sex 
structure  of  a  population  dataset.  Each  index  has  two  mandator/  characteristics. 
These  are: 

•  Its  identifier 

•  Its  size. 

The  index  identifier  is  used  to  reference  the  index  later  in  the 
program,  where  it  is  used  primarily  as  a  subscript  of  the  dataset(s)  that 
are  classified  by  it.  The  identifier  of  the  sex  index,  for  example,  may  be 
SEX. 

The  index  size  is  a  non-negative  integer  which  defines  the  maximum 
number  of  entries  (or  categories)  that  may  occur  within  the  classification 
scheme  of  the  index.  For  example,  the  size  of  the  SEX  index  is  2  for  the  two 
sex  classifications,  male  and  female. 

The  index  descriptors  are  the  user-supplied  descriptors  for  the 
index  entries.  For  example,  the  two  entries  of  the  SEX  index  may  be  defined 
as  MALE  and  FEMALE  or  as  SEX1  and  SEX2  or  in  some  other  way,  as  desired  by 
the  user.  The  descriptors  of  an  index  are  used  as  stubs  (row  descriptors) 
and/or  spanners  (column  headings)  of  tabular  displays  of  the  datasets 
classified  by  the  index.  The  descriptors  of  an  index  are  not  mandatory  for 
the  definition  of  the  index.  However,  they  are  indispensable  in  defining 
the  rows,  columns,  or  subtitles  of  the  tabular  displays  of  the  datasets  of 
the  program. 

The  entries  of  an  index  are  sequentially  ordered  from  1  to  n, 
where  n  is  the  size  of  the  index,  according  to  the  order  in  which  they  are 
read  into  the  program.  This  ordering  may  be  modified  by  the  SET  index 
statement  to  allow  selective  data  manipulation  or  data  retrieval  and  dis¬ 
play  operations. 


2.1.3 


Dataset: 


The  concept  of  the  dataset  is  central  in  the  system.  It  is  in  the 
datasets  of  the  program  for  an  application  that  the  numeric  information  is 
stored,  and  it  is  through  the  use  of  equations  and  other  instructions 
which  manipulate  these  datasets  that  the  logical  structure  of  the  application 
is  defined. 

A  dataset  must  have  a  unique  identifier  so  that  it  may  be  referenced 
in  subsequent  instructions  -  particularly  eatan' ore  -  within  the  procram. 

A  dataset  may  contain  only  a  single  numeric  value  -  in  which  case 
it  is  referred  to  as  a  scalar  -  or  it  may  have  multi-dimensional  structure. 

A  multi- dimensional  dataset  is  an  array  of  storage  cells  containing  numeric 
information  which  is  classified  (i.e.  ,  subscripted)  by  one  or  more  indexes. 

The  user  may  use  as  many  (ilO)  different  indexes  as  are  needed  to 
classify  the  values  of  a  given  dataset.  A  one-dimensional  array  -  sometimes 
referred  to  as  a  vector  -  is  a  sec  of  values  classified  by  a  single  index. 

For  example,  the  total  population  values  by  age  for  a  given  state  and  a  given 
year  may  be  stored  in  a  single  one-dimensional  dataset  that  is  classified 
by  an  age  index.  The  population  values  by  age  and  sex  for  a  given  state  and 
year  may  be  stored  in  a  single  two-dimensional  dataset  whose  two  dimensions 
are  Che  age  and  sex  indexes.  The  population  values  by  age,  sex,  and  state 
for  the  entire  U.S.  in  a  given  year  may  be  scored  in  a  single  three- 
dimensional  dataset.  The  third  dimension  in  this  dataset  is  the  state  index 
whose  size  is  50.  A  four-dimensional  population  dataset  may  nave  year  as  the 
fourth  index,  and  so  forth.  The  almost  unlimited  dimensionality  of  datasets 
is  an  important  property  of  the  system  as  compared  to  other  systems  which 
allow  only  three  dimensions  for  subscripted  arrays.  The  sizes  of  the  indexes 
of  multi-dimensional  datasets  are  not  restricted  by  the  user  language;  they 
are  restricted  only  by  the  size  of  central  computer  memory  available. 
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A  dataset  has  the  following  characteris tics : 

•  A  unique  identifier  (mandatory) 

•  Previously  defined  index(es)  as  subscripts.  Some  datasets 
(scalars)  have  no  indexes.  Up  to  10  unique  indexes  may 
classify  a  dataset.  The  order  of  dataset  indexes  is  impor¬ 
tant  in  READ  and  SHOW  statements  as  well  as  in  equations 
involving  the  dataset. 

•  An  optional  label,  which  appears  as  the  title  of  the  dataset 
in  later  SHOW  operations . 

•  An  array  of  values,  equal  in  number  to  the  product  of  the 
sizes  of  the  Indexes  classifying  the  dataset.  Initially, 
these  values  are  set  to  zero. 

In  addition  to  the  above,  there  is  one  other  very  important 
characteristic  of  the  dataset  array  of  values.  This  is  the  location  type 
for  those  values.  Datasets  are  said  to  be  fixed,  scratch,  or  IMS  depending 
upon  the  manner  in  which  the  dataset  values  are  assigned  a  location.  Fixed 
datasets  are  assigned  to  fixed  location  within  the  program  module  at  the  time 
the  dataset  is  defined.  Any  value  associated  with  that  dataset  will  remain 
constant  unless  explicitly  changed  by  the  user  even  if  the  program  module  is 
saved  and  loaded.  Scratch  data  sets  are  assigned  temporary  locations  when 
they  are  assigned  values.  They  maintain  these  locations  and  values  only 
during  the  current  execution  of  the  program  module.  Their  values  are  not 
maintained  when  their  program  module  is  saved;  thus,  scratch  datasets  are 
used  for  scratch  data  storage.  TMS  datasets  in  essence  have  two  locations  - 
one  within  a  TMS  array  on  some  data  file  and  a  temporary  one  within 
the  program  module.  TMS  datasets  have  all  of  the  advantages  of  fixed  data¬ 
sets  in  that  they  maintain  their  values  at  all  times  and  they  have  the 
advantages  of  scratch  datasets  in  that  they  do  not  require  a  fixed  amount  of 
space  within  the  program  module.  It  is  through  the  use  of  IMS  datasets  that 
large,  data  intensive  models  can  be  written. 


An  additional  benefit  of  the  use  of  IMS  datasets  is  that  the  sysca 
automatically  performs  all  space  allocations  required  to  bring  the  IMS  array 
values  in  and  out  of  the  program  module.  The  user  of  the  system  need  never 
concern  himself  with  the  logistics  of  the  operation^,  in  the  user  language. 
This  fact  makes  simple  the  writing  of  programs  within  this  user  language. 


2.1.4  Model 


The  model  is  a  group  of  instructions  that  are  compiled  as  a  unit 
under  a  single,  unique  identifier  for  later  reference  and  execution  as  a 
unit.  Model  definition  or  compilation  is  initiated  with  the  MODEL  statement 
and  is  terminated  by  the  END  statement.  Model  execution  is  initiated  by- 
entering  the  model  identifier.  The  model  statements  are  executed  sequential 
in  the  same  order  as  they  are  defined,  i.e.,  from  top  to  bottom. 

In  a  general  sense,  a  model  is  a  symbolic  or  mathematical  represen 
tation  of  an  actual  system.  The  model  instructions  are  usually  the 
equations  defining  the  logical  relationships  between  the  variables  (datasets 
of  the  system.  These  relationships  may  or  may  not  express  the  interaction 
of  variables  through  time . 

There  are  two  types  of  models  in  the  system: 

•  Static  models 

•  Dynamic  models . 

The  equations  in  static  models  do  not  contain  any  time  dependence.  The 
equations  in  a  dynamic  model ,  on  the  other  hand,  describe  the  evolution 
through  time  of  the  system  variables,  and  thus  of  the  actual  system  being 
modeled. 

Every  dynami c  model  has  four  implicit  parameters  associated  with 

time: 

•  The  time  parameter  itself 

•  DT,  Che  time  increment  parameter 

•  Beginning  time 

•  Ending  time. 
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The  tine  parameter  is  a  quasi-continuous  variable,  designated  by  the  key¬ 
word  TIME,  which  is  incremented  from  beginning  tine  to  ending  tine  in  steps 
of  DT,  the  time  increment  parameter.  DT  is  also  a  system  keyword. 

The  equations  in  a  dynamic  model  are  divided  into  three  sections: 
the  initial  section,  the  rate  section ,  and  the  level  section .  The  initial 
equations  define  the  values  of  the  model  variables  (or  levels )  at  beginning 
time .  The  rate  equations  define  the  values  of  those  variables  defining  the 
rates  of  growth  of  the  system  levels.  The  level  equations  define  the  values 
of  the  system  levels  at  an  arbitrary  time  point  in  the  interval  between 
. beginning  time  and  ending  time. 

The  specifics  of  dynamic  models  are  further  discussed  in  a  later 
part  of  this  handbook. 

2.1.5  Graph 

The  graph  is  defined  as  a  finite  set  of  ordered  pairs  of  numbers. 
These  number  pairs  may  be  thought  of  as  the  points  on  a  plane  representing 
some  functional  relationship,  y-f(x),  such  as  the  one  shown  below: 


This  graph  is  defined  as  the  set  of  the  six  x- values  and  the  six  v-values 
defining  the  x-  and  y-coordinates  of  the  six  points  denoted  by  dots  in  the 
above  schematic. 


A  graph  may  be  referenced  ir.  rhe  rignt-hanc  s 
wr.ere  it  ta.<es  various  values  depending  on  rhe  value  of 
rime  or  equation  execution.  The  value  //-value)  of  a  g 
argument  (x-vaiue)  is  automatically  computed  by  the  s-.-i 
intern olat ion  cr  extrapolation  between  the  points  derm 
example,  the  value  at  x  of  the  graph  shown  above  is  eq 
y-coordir.ate  corresponding  to  the  point  ?  on  the  crarh. 


rr  arv 


2.1.6  Table 


The  output  of  a  program  may  be  organized  into  tao.es  ,  or  tr.e 
user  may  display  each  program  dataset  directly  by  using  t.-.=  5 ho*  statement. 

A  table  is  first  defined  by  means  of  a  unique  identirier,  a  title,  a  sice 
specification  to  define  the  stub  (row  descriptors)  and  column  widens ,  ar.c 
a  list  of  the  datasets  whose  values  are  to  be  displayed  in  the  tabxe.  A 
defined  table  is  executed,  i.e.,  displayed  or  printed  out,  by  simply  entering 
its  identifier.  Table  execution  results  in  a  centered,  titled,  and  otherwise 
properly  structured  tabular  display  of  the  values  of  the  datasets  contained 
in  the  table.  The  tabular  displays  of  the  system  are  designed  to  fit  a 
standard  8-1/2  x  11  inch  page.  Large  tables  are  appropriately  segmented  into 
pages.  Table  paging  is  automatic.  To  satisfy  various  terminal  screen  cr 
printer  requirements,  the  length  and  width  of  cables  may  be  adjusted  by  the 
user. 

2.1.7  Program  Module 

A  Program  module  is  an  ordered  sec  of  user  language  statements  chat 
may  be  saved  for  later  loading  and  execution. 

A  program  module  is  bounded  by  a  START  or  a  LOAD  statement  at  its 
beginning  and  a  SAVE  or  a  STOP  statement  at  its  end. 


2.1.3  Equation 


An  equation  is  any  statement  that  begins  with  a'  previously 
defined  dataset  identifier,  is  followed  by  an  equal  (■)  sign,  and  contains 
on  the  right-hand  side  (RHS)  an  arithmetic  expression. 

Arithmetic  expressions  are  ordered  sequences  of  numerical  con¬ 
stants,  previously  defined  identifiers,  mathematical  operators,  and 
parentheses . 

Two  types  of  mathematical  operators  are  available  in  the  system: 

•  Arithmetic  operators 

•  Functional  operators. 


The  arithmetic  operators  of  the  system  are: 


Operator 

+ 


* 

/ 

** 


Operation 

Addition 

Subtraction. 

Multiplication 

Division 

Exponentiation . 


The  functional  operators  of  the  system  are: 


Operator 

SUM 

MAX 

MIN 

ABS 

LOG 

LN 

EXP 


Function 

Summation  over  subscripted  datasets  or 
arithmetic  expressions  with  subscripted 
datasets 

Finding  the  maximum  value  in  a  vector 
dataset 

Finding  the  minimum  value  in  a  vector 
dataset 

The  absolute  value  of  a  dataset 
The  base  10  logarithm  of  a  dataset 
The  natural  logarithm  of  a  dataset 
The  exponential  of  a  dataset. 
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.  execution,  the  val  ues  :  -c  t  -  ■  int 

KHS  of  ;on  must  be  well  defined,  i.e.  .  t:.c  ..'.jic.j  ^r.La;;:^  .r. 

the  expression  must  have  definite  values.  Upon  execution,  the  values  con¬ 
tained  in  the  storage  cells  of  the  dataset  on  the  left-hand  side  i LHS ,  are 
replaced  by  che  values  of  the  arithmetic  expression  on  the  RKS.  A  single 
identifier  nay  be  associated  with  different  numeric  values  at  difterent 
points  ir.  the  program  or  at  difterent  tines  during  execution. 

Lef c-co-Right  Hierarchy.  In  the  absence  of  parentheses,  arithme¬ 
tic  expressions  are  compiled  and  evaluated  from  left  to  right  with  each 
arithmetic  operator  using  the  preceding  result  as  an  operand.  All  arithme¬ 
tic  operators  are  of  equivalent  rank  in  the  processing  of  equations  by  the 
system. 

Note .  This  hierarchy  is  different  from  the  standard  algebraic 
hierarchy  used  in  FORTRAN. 

To  alter  the  implied  lef t-to- right  hierarchy  in  which  arithmetic 
operations  are  performed,  che  user  may  employ  parentheses.  Parentheses 
allow  the  user  to  specify  his  own  order  for  the  interpretation  and  execution 
of  arithmetic  expressions.  An  arithmetic  expression  (or  a  portion  of  one) 
beginning  with  a  left  parenthesis  and  ending  with  a  right  parenthesis  is 
called  a  parenthetical  group .  A  parenthetical  group  is  said  to  be  nes  ted 
if  it  is  contained  within  another  parenthetical  group. 

The  operations  of  the  most  deeply  nested  parenthetical  group  of 
an  arithmetic  expression  are  executed  first.  Operations  within  a  paren¬ 
thetical  group  are  executed  from  left  to  right.  Parenthetical  groups  of 
equivalent  nesting  are  processed  from  left  to  right. 
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la  general,  the  rules  for  constructing  well  formed  equations 

are : 

(1)  No  two  operators  may  directly  follow  one  another. 

(2)  No  two  constants  or  variables  may  directly  follow 
one  another. 

(3)  The  right-hand  side  of  an  equation  must  begin  and 

end  with  a  variable,  constant,  or  parenthetical  group. 

(4)  A  parenthetical  group  may  appear  only  on  the  right- 
hand  side  of  an  equation. 

(5)  Right  and  left  parentheses  must  always  be  balanced, 
i.e.,  be 'equal  in  number. 

(6)  A  left  parenthesis,  (,  must  be  preceded  by  another 

(,  an  operator,  or  an  equal  sign,  and  must  be  followed 
by  another  (,  an  identifier,  or  a  constant. 

(7)  A  right  parenthesis,  ),  must  be  preceded  by  another 
),  an  identifier,  or  a  constant,  and  be  followed  by 
another  ) ,  or  an  operator . 

(8)  An  equation  may  be  continued  to  the  succeeding  line 

as  long  as  an  operator  or  a  comma  is  the  last  character 
on  the  line(s)  being  continued. 

(9)  Equations  are  executed  from  left  to  right. 

2.1.9  Data  File 


The  data  file  makes  it  possible  for  the  user  to  extend  the  storage 
available  to  him  from  the  actual  central  memory  of  the  computer  to  off-line 
direct  storage  devices  such  as  the  disk.  There  axe  two  very  important 
reasons  why  this  extension  is  needed  for  large-scale  modeling  applications. 


In  the  first  place,  ic  is  required  that  models  written  by  one  group  : : 
researchers  be  accessed  and  used  by  others  at  some  later  point  in  time. 
Off-line  storage  is  needed  to  perform  this  saving  operation.  In  the  second 
place,  large-scale  models  are  usually  data  intensive,  often  manipulating 
millions  of  numeric  values  during  a  single  execution  cycle.  It  would  be 
impossible  to  store  all  of  these  values  within  the  central  memory  of  the 
machine.  Off-line  storage  is  needed  to  extend  the  numeric  value  storage 
area  required  by  the  execution  of  the  models. 

As  is  implied  in  the  above,  there  are  two  basically  different  types 
of  data  files  used  by  the  system.  These  are  referred  to  as  save  files  and 
tab le  management  svstem  (or  TMS )  files.  Save  files  are  used  tor  the  storage 
of  program  modules.  They  are  needed  primarily  so  that  previously  written 
program  modules  can  be  saved  and  loaded  for  later  use.  Table  management 
system  files  are  used  for  the  storage  of  TMS  arrays.  They  are  needed 
primarily  to  contain  the  values  of  TMS  datasets  and  thus  to  extend  the 
numeric  value  storage  available  to  the  executing  models.  TMS  files  also,  of 
course,  allow  the  saving  of  these  numeric  values  for  later  use. 

It  is  through  the  combined  use  of  save  files  and  table  management 
system  files  that  this  system  achieves  its  powerful  data  base  management 
capabilities.  These  data  base  management  capabilities  make  it  possible  to 
easily  structure  and  implement  large  complex  models  and  their  associated 
information  systems  in  such  a  manner  that  they  are  easy  to  use. 

Associated  with  each  data  file  there  is  a  physical  unit  number . 

This  physical  unit  number  is  used  to  identify  the  particular  file  to  the 
operating  system.  Before  a  given  file  can  be  used,  it  must  be  declared 
to  the  operating  system  via  its  physical  unit  number,  using  the  directives 
of  the  operating  system.  Prior  to  using  this  system  the  user  should  ootain 
information  from  his  computer  staff  as  to  the  required  form  of  these 
directives . 
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In  addition  co  the  physical  unit  number  there  is  also  a  size 
associated  with  each  data  file.  Prior  to  creating  any  new  data  file,  the 
user  must  first  calculate  the  total  amount  of  space  which  he  will  require. 
Based  on  this  calculation  he  can  then  either  initially  declare  the  file  so 
that  it  has  the  appropriate  si2e  when  he  defines  it  to  the  operating 
system  or  he  can  choose  the  appropriate  physical  unit  number  for  a  physical 
file  which  has  the  appropriate  size.  The  particular  manner  in  which  either 
of  these  is  accomplished  depends  upon  the  operating  system.  The  user  should 
consult  with  his  computer  staff  before  beginning  the  creation  of  a  new  data 
file. 


2.1.10  Table  Management  System  Array 

The  table  management  system  (IMS)  array  is  the  basic  unit  on  a  TMS 
data  file.  It  may  be  thought  of  as  an  n-dimensional  (0<n<10)  array  of 
numeric  values.  In  this  sense  it  is  precisely  like  a  dataset.  The  primary 
purpose  of  the  TMS  array  is  to  contain  values  for  TMS  datasets. 

A  TMS  array  has  the  following  characteristics: 

•  A  unique  identifier  consisting  of  a  data  base  number 
and  an  array  number  as  described  in  Subpart  2.1.1. 

•  A  size  specification  which  is  a  sequence  of  from  zero 
to  ten  integer  values  which  define  the  number  of 
dimensions  in  the  array  and  the  number  of  entries  in 
each  dimension.  An  array  with  zero  dimensions  consists 
of  a  single  value.  The  total  number  of  values  associated 
with  a  multi-dimensional  array  equals  the  product  of  the 
number  of  entries  in  each  dimension. 
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2.2  General  Svma::  Rules 


The  User  Language  statements  are  listed  in  Table  1-1.  This  parr 
describes  the  syntax  of  the  User  Language,  i.e.,  the  manner  in  which  the 
statement  words  in  this  table  are  combined  with  other  keywords,  user- 
supplied  identifiers,  strings,  and  special  characters  to  form  instructions 
that  are  acceptable  for  processing  by  the  system. 

Keywords  are  reserved  words  in  the  User  Language  that  are  used  to 
specify  more  completely  the  statement  beginning  with  one  of  the  statement 
words  in  Table  1-1.  For  example,  the  word  TITLE  is  a  keyword  used  in  the 
definition  of  a  TABLE  statement  and  it  allows  the  user  to  specify  a  title 
for  the  table  being  defined.  Similarly,  the  word  LABEL  is  a  keyword  used 
in  dataset  definitions  to  allow  the  user  to  specify  dataset  labels  or 
descriptors . 

2.2.1  Character  Set 


The  system  interpreter  recognizes  the  47  characters  listed  in 
Table  1-3.  Each  character  is  associated  with  a  character  sequence  number. 
Some  familiarity  with  these  sequence  numbers  is  helpful  for  understanding 
how  the  system  interprets  input. 

The  characters  in  Table  1-3  may  be  divided  into  the  six  groups 
defined  in  Table  1-4. 

Alphabetic  characters  are  used  to  define  and  reference  all 
identifiers  in  a  program.  Each  identifier  must  have  an  alphabetic  character 
as  its  first  character. 

Numeric  characters  are  used  to  formulate  all  numeric  constants  of 
the  program  -  values  of  datasets  and  graphs,  index  sizes,  index  settings, 
dataset  subscript  values,  format  specifications,  etc.  They  may  also  be 
used  as  secondary  characters  in  the  formulation  of  identifiers. 

The  sign  characters,  +  and  are  used  to  define  the  arithmetic 
operations  of  addition  and  subtraction,  respectively.  They  are  also  used 
to  define  the  sign,  positive  or  negative,  of  numeric  constants  or  cr it.. me tic 
expressions . 


TABLE  1-3.  THE  CHARACTER  SET  OF  THE  SYSTEM 
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Seq 

Char 

Seq 

Char 

Seq 
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1 

A 

13 

M 

25 

Y 

37 

+ 

2 

B 

14 

N 

26 

Z 

38 

- 

3 

C 

15 
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27 

0 

39 

it 

4 

D 

16 

P 

28 

1 

40 
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18 

R 

30 

3 

42 

) 

7 

G 

19 

S 

31 

4 

43 

$ 

8 

H 

20 

T 

32 

5 

44 

- 

9 

1 

21 

U 

33 

6 

45 
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10 

J 

22 

V 

34 

7 

46 

» 

11 

K 

23 

w 

35 

8 

47 

• 

12 

L 

24 

X 

36 

9 

TABLE  1-4.  CHARACTER  GROUPINGS 


Group 

Type 

Characters 

Sequence  Numbers 

1 

Alphabetic 

A  to  Z 

1 

to 

26 

2 

Numeric 

0  to  9 

27 

to 

36 

3 

Sign 

+  - 

37 

to 

38 

4 

Special 

*  /  (  )  ! 

?  m 

39 

to 

44 

5 

Blank 

blank  , 

45 

to 

46 

6 

Decimal 

• 

47 

operators . 


oe  usea 


The  special  characters  x  ar.c  arc  nuitip-i  cation  ar.c  divisi : 
In  addition,  for  batch  use  of  t.-.e  sv?  cam,  may  be  used  to  skip 
to  the  next  page  in  the  program  listing  while  an  input  line  having  *  in  the 
first  column  indicates  a  line  of  comment  to  be  printed  but  not  processed  by 
the  interpreter.  The  character  sequence  **  is  che  exponentiation  operator. 
The  characters  (  and  )  are  used  to  enclose  program  descriptors,  dataset 
subscripts,  graph  arguments,  parenthetical  groups  in  aritnnetic  operations, 
etc.  The  =  character  is  used  in  equation  statements. 

All  special  characters  may  be  used  freely  in  defining  descriptor 
strings.  Special  characters  may  not  be  used  in  identifiers. 

The  two  blank  characters,  blank  and  comma,  serve  as  TTZ-breakers , 
i.e.,  they  separate  the  Textual  I_nput  Elements  (TIEs)  in  user  language 
statements.  In  this  context,  they  are  very  important  ingredients  of  the 
interpreter  of  the  system.  The  comma  is  also  used  as  a  continuation 
character  to  continue  a  line  of  input  to  the  next  line. 

The  decimal  character  is  used  to  specify  real  constants. 


2.2.2  TIE  -  Textual  Input  Element 


All  statements  are  recognized  and  interpreted  by  the  system  as 
ordered  sequences  of  Textual  Input  Elements  (TIEs)  which  are  separated 
from  each  other  by  TIE-breakers . 

A  TIE  is  any  sequence  of  characters  which  are  in  some  sense 
logically  related.  Thus,  the  sequence  of  characters  FRED,  preceded  and 
followed  by  blanks,  may  be  thought  of  as  a  single  TIE  which  names  something. 
Similarly,  the  sequence  of  characters  3.14159,  preceded  and  followed  by 
blanks,  may  be  thought  of  as  a  single  TIE  whose  value  approximates  the  value 
of  IT . 

Consider,  however,  the  sequence  SHOW  POPT.  Is  this  sequence  one 
single  TIE?  Or,  is  it  cwo  TIEs  separated  by  a  blank?  The  second  inter¬ 
pretation  appears  to  be  more  correct,  because  the  sequence  SHOW  is  a  request 
to  show  something  and  POPT  is  the  name  of  whatever  is  to  be  shown.  ~iis 
example,  demonstrates  the  need  for  the  notion  of  a  TIE-breaker.  A  1IE- 
breaker  is  any  aspect  of  the  input  stream  which  causes  a  TIE  to  end.  Thus, 
the  TIEs  FRED,  3.14159,  SHOW,  and  POPT  are  broken  by  the  blanks  r r  commas 
following  them. 
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The  notions  of  TIE  and  TIE-breaker  cannot  be  considered  inde¬ 
pendently  of  each  other.  Both  are  needed  to  understand  the  formulation  of 
system  statements.  Any  system  statement  is  interpreted  from  left  to  right 
as  a  string  of  recognizable  TIES  -  words,  numbers,  etc.  -  comprised  of  alpha 
betic,  numeric,  and  special  characters,  separated  from  each  other  by  TIE- 
breakers,  and  ordered  according  to  the  rules  of  syntax  spelled  out  in  this 
handbook.  It  is  the  recognition  of  distinct  TIEs  which  is  crucial  to  the 
understanding  of  sequences  of  characters,  not  just  the  recognition  of  the 
characters  themselves. 

TIEs  may  be  classified  into  four  types  depending  on  the  types  of 
individual  characters  that  make  them  up.  The  four  types  of  TIEs  may  be 
defined  in  terms  of  the  six  groups  of  characters  defined  in  Table  1-4. 

A  Type  _1  TIE  is  either  a  single  special  character,  or  a  sign 
character  which  is  not  followed  by  a  numeric  character,  or  a  decimal  charac¬ 
ter  which  is  neither  followed  nor  preceded  by  a  numeric  character.  A  Type  1 
TIE  consists  always  of  a  single  character. 

A  Type  2_  TIE  is  a  sequence  of  only  alphabetic  and  numeric  charac¬ 
ters,  i.e.,  it  is  used  as  a  keyword  or  as  a  name.  The  first  character  of 
a  Type  2  TIE  must  be  alphabetic.  The  first  nonblank  TIE  of  any  statement 
is  a  Type  2  TIE.  Also,  all  index,  dataset,  and  model  identifiers  must  be 
unique  Type  2  TIEs. 

A  Type  3  TIE  is  an  integer,  i.e.,  it  is  a  sequence  of  numeric 
characters  which  may  be  preceded  by  a  single  sign  character  or  by  a  sign 
character  followed  by  blank  characters.  A  sequence  of  numeric  characters 
preceded  by  alphabetic  characters  is,  or  course,  a  Type  2  TIE,  and  not 
Type  3 . 

A  Type  4  TIE  is  a  real  number.  It  consists  of  a  sequence  of 
numeric  characters  containing  one  and  only  one  decimal  point.  A  Type  4 
TIE  sequence  may  be  preceded  by  a  single  sign  character  or  by  a  single 
sign  character  followed  by  blank  characters. 

A  statement  is  converted  into  a  sequence  of  values  corresponding 
to  the  various  types  of  TIEs  comprising  it,  except  for  blanks  and  commas. 
Blanks  and  commas  have  no  internal  values .  Blank  characters  are  either 
ignored  (except  in  descriptors)  in  an  input  sequence  or  they  serve  as  TIE- 
breakers . 


By  definition,  a  TIE  of  a  given  type  is  net  ended  unless  it  is 
followed  by  a  blank  character  or  by  a  TIE  of  a  different  type.  A  Tin- 
breaker,  therefore,  is  either  a  blank  character  or  the  beginning  character 
of  a  TIE  of  a  different  type. 

The  TIE  and  TIE-breaker  notions  are  not  only  central  to  the 
interpreter  of  the  system;  they  are  also  linked  to  the  error  processing 
operations  of  the  interpreter.  Usually,  the  error  diagnostics  issued  by 
the  interpreter  are  expressed  in  terms  of  the  Ills  in  the  input  statement 
that  are  in  error. 

2.2.3  Strings  and  Descriptors 


In  addition  to  TIEs,  the  other  basic  syntactic  element  in  terms 
of  which  the  system  recognizes  and  interprets  statements  is  the  string. 

Every  statement  in  a  program  is  viewed  as  an  ordered  sequence  of  TIEs  and 
strings.  A  string  is  an  arbitrary  sequence  of  characters.  To  qualify  as 
a  string,  a  sequence  of  characters  must  be  enclosed  in  parentheses. 

The  Interpretation  of  a  sequence  of  characters  enclosed  in 
parentheses  as  a  TIE  or  a  string  is  context  sensitive,  i.e.,  it  depends  on 
the  relative  location  within  the  type  of  statement  in  which  it  occurs. 

All  structural  elements  of  the  language  have  mandatory  identifiers, 
which  are  mnemonics  of  up  to  5  alphanumeric  characters.  These  mnemonics  are 
abbreviated  names  of  the  structural  elements.  For  additional  program 
documentation  these  identifiers  may  be  further  described  by  explicit 
descriptors.  For  example,  the  identifier  POPM  for  a  dataset  of  population 
values  may  be  described  by  the  label  or  descriptor  MALE  POPULATION'. 

All  labels  or  descriptors  of  the  language  are  character  strings 
and  are  specified  as  part  of  the  definitions  of  the  structural  elements 
which  they  define. 


/ 
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2.2.4  Line  Format 


The  User  Language  is  generally  a  free-form  language.  Language 
statements  may  be  begun  anywhere' on  the  input  line  or  card,  and  as  many 
blanks  or  commas  as  desired  may  be  inserted  between  the  various  TIEs  to 
improve  statement  readability. 

2.2.5  Continuation 


If  a  given  statement  will  not  fit  on  a  single  card  or  input  line, 
then  continuation  of  that  statement  may  be  indicated  in  one  of  three  ways, 
depending  upon  the  context. 

(1)  If  the  system  is  currently  reading  a  character 
string,  then  continuation  is  automatic.  The 
first  character  of  the  following  card  or  input  line 
is  concatenated  directly  behind  the  last  character 
of  the  preceding  card  or  input  line. 

(2)  If  the  system  is  currently  reading  an  equation, 
then  continuation  to  the  next  cards  or  input  lines 
will  be  assumed  if  and  only  if  the  last  nonblank 
character  of  the  current  card  or  line  is  a  comma  or 
an  arithmetic  operator. 

(3)  If  the  system  is  currently  reading  neither  a  string 
nor  an  equation,  then  the  continuation  to  the  next 
card  or  input  line  will  be  assumed  if  and  only  if 
the  last  nonblank  character  on  the  current  card  or 
input  line  is  a  comma. 

2.2.6  Line  Length 

The  length  assumed  by  the  system  for  each  input  line  depends 
entirely  on  the  source  of  those  lines.  If  input  is  coming  from  any  source 
other  than  the  standard  input  file  connected  to  the  user  terminal  in 


interactive  node,  then  each  input  line  is  assumed  to  contain  exactly  30 
characters.  If  a  physical  input  line  contains  more  than  30  characters,  it 
will  be  truncated.  If  it  contains  less  than  SO  characters,  it  '..'ill  be 
padded  out  to  SO  with  blanks. 

If  the  input  is  coming  from  che  user  terminal,  the  situation  is 
more  complex.  Basically,  each  input  line  is  as  long  as  the  number  of  charac 
ters  typed  by  the  user  after  the  system  prompt  and  before  entering  carriage 
return.  Thus,  the  line  length  assumed  by  the  system  is  independent  of  the 
actual  line  spacing  being  generated  at  the  particular  terminal  being  used. 
The  user  must  be  careful  to  make  any  continuation  provisions  before  he 
enters  carriage  return  which  is  not  necessarily  before  he  exhausts  a  given 
line  at  the  terminal.  The  maximum  length  of  any  line  entered  at  the  termina 
is  80  characters. 

2.2.7  Prompting  and  Echo  Printing 

The  decision  as  to  whether  to  prompt  the  user  for  the  next  input 
line  or  to  echo  print  an  input  line  after  it  has  been  read  depends  entirely 
on  the  source  of  this  line.  If  input  is  coming  from  any  source  other  than 
the  standard  input  file  connected  to  the  user  terminal  in  interactive  mode, 
the  system  will  echo  print  that  input.  If  the  input  is  coming  from  the 
user  terminal,  then  the  system  will  prompt  for  that  input. 

For  the  interactive  or  on-line  user,  the  system  issues  a  prompt 
before  each  instruction.  The  interactive  prompt  for  a  language  statement 
has  the  following  format: 

READY/ 

The  interactive  prompt  for  input  data  following  a  READ  statement 
is : 

DATA/ 

For  the  batch  user,  the  system  echo  prints  all  the  instructions 
in  the  program  following  the  very  first  instruction 

READY/ SYSTEM  USER-OFF 

which  initiates  batch  mode.  This  instruction  is  not  echo  print  d. 


Each  batch  instruction  is  preceded  by  two  decimal  integers  speci 
fying  the  starting  address  in  core  where  this  instruction  is  stored  and 
the  last  available  address  in  core  that  can  be  used  for  the  storage  of 
values  associated  with  scratch  and  TMS  data  sets.  All  echo-printed  input 
data  cards  following  READ  instructions  are  preceded  by  the  word  DATA. 


The  5tatsner.es  or  the  User  Language  are  listed  tr.  Tatis  1-1.  Lost 
of  these  statements  are  used  in  the  formulation  or  me  sample  programs  shown 
in  Chapter  5.0.  In  this  chapter,  the  User  Language  statements  are  describe: 
in  terms  of  their  general  syntax  in  the  same  order  as  shown  in  Table  1-1. 

Each  statement  description  contains  a  summary  description  of  me 
basic  function  of  the  statement,  the  general  syntax  for  the  statement,  ar.u 
descriptions  of  all  options,  secondary  keywords ,  and  user-supplied  identifie 
descriptors,  and  parameters  used  ir.  the  general  statement  syntax. 

The  general  statement  syntax  consists  of  the  following: 

•  The  statement  word  -  always  the  firs:  word  of  the 
statement 

•  Secondary  keywords  or  options 

•  Special  characters 

•  User-supplied  identifiers,  descriptors,  and 
parameters . 

The  statement  words  of  the  language  are  listed  in  Table  1-1.  All 
language  statements  begin  either  with  one  of  the  statement  words  in  Table  I- 
or  with  a  previously  defined  identifier  in  the  case  of  equations,  model 
execution  statements,  or  table  execution  statements.  In  this  sense  stacemen 
words  are  primary  keywords  of  the  language. 

The  secondary  keywords  of  a  statement  are  reserved  words  in  the 
language  that  allow  the  user  to  specify  various  options  or  parameters 
associated  with  the  general  function  of  the  statement. 

The  meaning  of  secondary  keywords  is  context  sensitive  and  always 
relative  to  the  statement  to  which  they  belong.  For  example,  the  word  TITLE 
is  a  keyword  used  with  the  TABLE  statement  to  allow  the  user  to  specify  a 
title  for  the  table  being  defined  by  the  TABLE  statement. 

Some  secondary  keywords  are  mandatory  for  the  complete  or  correct 
specification  of  a  statement  and  some  are  optional. 
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Special  characters  -  such  as  blanks,  commas,  periods,  parentheses, 
asterisks,  etc.  -  are  used  in  the  formulation  of  a  statement  to  separate  or 
link  various  components  of  the  statement  or  to  define  specific  statement 
options . 

The  meaning  of  a  special  character  is  also  context  sensitive  and 
it  is  related  to  the  statement  to  which  it  belongs.  For  example,  the  use 
of  the  asterisk,  *,  in  an  equation  means  "multiply",  while  its  use  in  a  SET 
index  statement  means  "restore  index  setting  to  its  original  size  and 
ordering" . 

Even  within  a  statement,  a  special  character  may  have  different 
meanings.  In  the  following  example  defining  the  dataset  AGR 

DATA  AGR,  LABEL  (ANNUAL  GROWTH  RATE  (S/YEAR)) 

the  first  set  of  parentheses  enclosing  the  label  is  mandatory  and  part  of 
che  syntax  of  the  DATA  statement,  while  the  second  set  of  parentheses 
enclosing  5/YEAR  is  optional  and  part  of  the  user-supplied  descriptor.  This 
descriptor  could  have  been  stated  without  these  parentheses 

DATA  AGE,  LABEL  (ANNUAL  GROWTH  RATE  IN  5 /YEAR) 

Some  special  characters  are  mandatory  and  some  are  optional  for 
the  complete  or  correct  specification  of  a  statement. 

In  this  manual,  the  syntax  of  a  statement  or  an  option  is  shown 
enclosed  in  a  box. 

The  statement  word  and  all  secondary  keywords  defining  various 
statement  options  are  capitalized  in  the  syntax  of  each  statement. 

All  user-supplied  identifiers,  descriptors,  and  parameters  are 
in  lower  case. 

Identifiers  are  always  alphanumeric  and  may  contain  up  to  5 
characters,  the  first  of  which  must  be  alphabetic.  They  must  contain  no 
special  characters.  The  notation  for  an  identifier  in  a  syntax  description 
is  "id". 

Descriptors  are  character  strings  that  may  contain  up  to  125 
characters.  All  characters  in  the  character  set  of  the  system  are  allowed. 
Descriptors  are  always  enclosed  in  parentheses.  The  notation  for  a  cescrip- 
tor  in  a  syntax  description  is  "string". 
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Parameters  are  usually  numeric  values  chat  define  sites  icr  various 
structural  elements  or  various  format  specifications.  Parameters  are  usually 
abbreviated  by  one  or  two-letter  symbols. 

The  notation  "list"  in  syntax  descriptions  means  a  list  o:  items 
separated  by  blanks  or  commas. 

The  statements  of  the  User  Language,  as  listed  in  Table  1-1,  are 
divided  into  three  groups: 

e  Program  control  statements 

•  Program  structure  definition  statements 

•  Data  manipulation  statements 

These  language  statements  are  presented  in  the  same  order  as  they  are  lifted 
in  Table  1-1. 


SYSTEM 
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3. 1  The  SYSTEM  Statement 


This  is  Che  basic  scatement  establishing  c'ne  control  parameters 
for  the  user-syscem  interface.  The  user-system  interface  is  defined  in  terms 
of  a  set  of  parameters  that  specify  the  following: 

•  Batch  or  on-line  use  of  the  system 

•  Debugging  or  production  use  of  the  system 

•  Width  of  output  page 

•  Length  of  output  page 

e  Unit  numbers  for  off-line  files  on  which  to  save  program 
modules,  from  which  to  load  previously  saved  program 
modules,  from  which  to  read  data  inputs,  or  which  con¬ 
stitute  the  files  within  the  data  base  set 

•  Various  options  for  report  generation 

•  Various  TMS  data  file  options. 

Syntax 

SYSTEM  parameter,  setting 


parameter  is  any  one  of  the  following  keywords :  USER,  DEBUG,  WIDTH, 
LINES,  INPUT,  SAVE,  LOAD,  LIST,  ZERO,  EPS,  SPACE,  DEMO, 
DB1,  DB2,  DB3,  INITIAL,  or  HEADING. 

setting  is  the  particular  parameter  setting  selected  by  the  user. 

Comments 


The  parameter  settings  are  local,  i.e.,  they  remain  in  effect  only 
until  they  are  reset  in  a  subsequent  SYSTEM  parameter  scatement. 


SYSTEM 


The  user  nay  sec  multiple  parameters  with  a  sing.e  SYSTEM 


s  cafement . 


model . 


The  SYSTEM  statement  may  be-issuea  within  or  outside  of  a 


3.1.1  The  USER  Parameter 


[ 


This  parameter  allows  the  user 
with  the  system  from  ics  default  on-line 
vice-versa . 


to  modify  the  mode  of  interaction 
setting  to  the  batch  setting,  or 


Syntax 

A.  For  on-line  use, 

SYSTEM  USER=ON 

B.  For  batch  use, 

SYSTEM  USER-OFF 

The  ON  setting  is  the  default. 


Comments 

Since  the  default  mode  of  interaction  with  the  system  is  ON,  the 
first  statement  of  any  batch  program  must  be: 

SYSTEM  USER-OFF 

This  statement  must  precede  all  program  statements,  including  the  START 
statement. 

The  keyword  USER  itself  may  be  optionally  omitted. 


-jST 


SYSTEM 
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3.1.2  The  DEBUG  Parameter 


This  parameter  allows  the  user  to  find  the  synt-ax  errors  in  his 
program  module  in  a  single  pass  of  the  system.  This  feature  simplifies  the 
debugging  of  complex  batch  programs.  This  parameter  is  relevant  only  when 
the  USER  parameter  is  off.  If  the  USER  parameter  is  off  and  if  the  system 
is  not  in  DEBUG  mode,  then  the  system  will  end  execution  after  a  syntax 
error  is  encountered.  If  the  system  is  in  DEBUG  mode,  it  will  continue 
processing  even  if  errors  are  encountered. 

Syntax 

A.  To  place  the  system  into  DEBUG  mode, 

SYSTEM  DEBUG-ON 


B.  To  take  the  system  out  of  DEBUG  mode, 


SYSTEM  DEBUG-OFF 


The  OFF  setting  is  the  default. 

3.1.3  The  WIDTH  Parameter 


This  parameter  allows  the  user  to  modify  the  default  output  from 
72  characters  per  line  to  the  desired  width.  Given  the  variety  of  printer 
widths  and  terminal  screen  sizes,  the  user  must  be  able  to  control  this 
parameter,  otherwise  his  output  may  be  malformed. 

Syntax 


SYSTEM  WIDTH-n 


n 


is  an  integer  (40  n  £  135)  defining  the  desired  printer 
with  (or  width  of  terminal  screen)  in  characters. 

Default  n-72. 


SYSTEM 
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3.1.4  The  LINES  Parameter 


This  parameter  allows  the  user  to  modify  the  default  output  page 
length  from  55  lines  per  page  to  the  desired  length.  Giver,  the  variety  of 
terminal  screen  sizes  or  printer  page  lengths,  the  user  must  be  able  to 
control  this  parameter,  otherwise  his  output  may  be  malformed. 

Svncax 


SYSTEM  LINES-n 


n  is  a  positive  integer  defining  the  desired  length  in  lines  of 

each  output  page. 

Default  n*55. 

3.1.5  The  INPUT  Parameter 


This  parameter  allows  the  user  to  specify  the  unit  number  assigned 
to  an  off-line  data  file  from  which  subsequent  READ  operations  will  read 
input  data.  This  parameter  allows  the  user  to  read  data  inputs  from  any 
previously  constructed  and  appropriately  "attached"  off-line  data  file,  other 
than  the  default  input  file  of  the  local  computer  facility. 

Syntax 


SYSTEM  INPUT-n 


n  is  an  integer  (l£n^  99)  denoting  the  unit  number  of  the  off-line 

data  file  containing  coded  input  for  subsequent  READ  operations. 


Comments 


The  "attach"  procedures  by  which  the  unit  number  n  is  assigned 
to  a  data  file  are  outside  the  scope  of  this  handbook.  They  are  described 
in  the  local  user  manuals  of  the  computer  facility  where  the  system  resides. 


SYSTEM 
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3.1.6  The  SAVE  Parameter 


This  parameter  allows  the  user  to  modify  the  default  physical 
unit  number,  1,  assigned  to  the  local  file  where  program  modules  may  be 
saved  for  later  loading.  This  local  file  is  referred  to  as  the  "save"  file. 

Svntax 


SYSTEM  SAVE-n 


n  is  an  integer  (l<_n_^_99)  denoting  the  physical  unit  number  of  the 

"save"  file. 

Default  n**l. 


Comments 


The  "save"  file  cannot  be  "cataloged",  i.e.,  saved  permanently, 
unless  previously  "requested".  The  "catalog"  and  "request"  procedures  are 
outside  the  scope  of  this  handbook. 

3.1.7  The  LOAD  Parameter 


This  parameter  allows  the  user  to  modify  the  default  physical 
unit  number,  2,  assigned  to  the  local  file  from  where  previously  saved 
modules  may  be  loaded.  This  local  file  is  referred  to  as  the  "load"  file. 

Syntax 


SYSTEM  LOAD-n 


n 


is  an  integer  (l<_n^_99)  denoting  the  unit  number  of  the  "load" 
file. 

Default  n»2. 


SYSTEM 


•ar2Ke:-2r 


.his  parameter  controls  the  echo  printing  of  inputs  and  is  rele¬ 
vant  only  if  the  USER  parameter  is  OF?,  i.e.,  only  in  baton  mode. 

Svntax 


A.  To  echo  print  all  input  consecutively, 


SYSTEM  LIST-ON 


This  is  the  default  setting. 


B.  To  echo  princ  next  input  line  on  top  of  next  page, 


SYSTEM  LIST-EJECT 


Comments 


When  list  is  set  to  EJECT,  the  system  changes  it  to  ON  once  the 
next  input  line  has  been  printed  at  the  top  of  the  page. 

3.1.9  The  ZERO  Parameter 

This  parameter  allows  the  user  to  control  the  display  of  "zero" 
entries  in  output  tables . 


SYSTEM 


Syntax 
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A.  To  display  "zero"  entries  using  the  character  0  (or  0.00...), 


SYSTEM  ZERO-O 


B.  To  display  "zero"  entries  using  special  characters,  - 


SYSTEM  ZERO (char) 


char  is  a  string  of  up  to  5  characters. 

Comments 


The  default  string  for  the  symbolic  display  of  "zero"  entries 
is  0  (or  0.00...,  depending  on  the  number  of  decimals  specified). 


3.1.10  The  EPS  Parameter 


This  parameter  allows  the  user  to  control  the  display  of  arbitrarily 
small  ("epsilon")  entries  in  output  tables.  Here  "epsilon"  means  "zero  within 
the  level  of  significance  of  the  display". 

Syntax 

A.  To  display  "epsilon"  entries  as  0  (or  0.00...), 

SYSTEM  EPS-0 


B.  To  display  "epsilon"  entries  using  special  characters. 


SYSTEM  EPS (char) 


char 


is  a  string  of  up  to  5  characters. 


er.t 


Mb.L. 


-  _  3  ~ 


The  default  string  for  the  symbolic  display  of  "epsilon" 
is  0  for  0.00...,  depending  on  the  number  of  decimals  specified). 


ries 


3.1.11  The  SPACE  Parameter 


The  report  generator  within  che  system  produces  several  blank 
lines  within  the  heading  of  reports  to  increase  legibility.  At  times  it  is 
convenient  to  suppress  blank  lines  so  as  to  increase  the  number  of  lines 
which  can  be  printed  on  a  given  page.  The  SPACE  parameter  is  used  to  suppress 
blank  lines  used  for  spacing  purposes. 


Svntax 


A.  To  turn  spacing  suppression  on, 


SYSTEM  SPACE-ON 


B.  To  turn  spacing  suppression  off. 


SYSTEM  SPACE-OFF 


This  is  the  default  setting. 

3.1.12  The  DEMO  Parameter 


This  parameter  is  used  to  allow  the  on-line  user  to  interrupt 
the  printing  of  a  long  table  and  to  select  from  the  following  options: 

•  Continue  the  printing 

•  Skip  che  next  page  of  output 

•  Skip  the  remaining  table  entirely  and  issue  the  next  prompt. 


SYSTEM 
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Syntax 

A.  To  allow  printing  interruption. 


SYSTEM  DEMO-ON 


3.  To  print  tables  without  allowing  for  interruptions, 


SYSTEM  DEMO-OFF 


■  This  is  the  default  setting. 

Comments 


When  the  ON  setting  is  in  effect,  the  printing  interruption  takes 
place  at  the  end  of  an  actual  output  page  and  a  slash,  /,  prompt  is  issued. 

At  the  prompt,  if  the  user  enters  a  blank,  then  the  printing  will 
continue  with  the  next  page.  If  the  user  enters  S,  then  the  next  page  is 
skipped,  and  printing  continues  with  the  following  page.  If  the  user  enters 
anything  else,  then  the  remaining  table  is  skipped  entirely. 

This  option  applies  only  in  on-line  mode,  and  it  is  useful  when 
giving  "demonstrations"  of  long  printouts,  hence  the  keyword  DEMO. 

3.1.13  The  DB1,  DB2,  and  DB3  Parameters 

These  parameters  allow  the  user  to  modify  the  default  physical 
unit  numbers  assigned  to  the  three  table  management  system  files  within  the 
data  file  set. 


SYSTEM 


!  SYSTEM  DBl»n  1 

i  [ 

j  SYSTEM  DB2=n  I 
I  SYSTEM  DB2=n  ! 

i — _ 1 

n  is  an  integer  (l<n<99)  denoting  the  physical  unit  number 

of  a  TMS  data  file 

Default  DB1=1,  DB2=2.  DB3=3. 


3.1.14  The  INITIAL  Parameter 


This  parameter  allows  the  user  to  control  the  allocation  of 
working  storage  for  dataset  values. 

Svntax 


A.  To  allocate  all  dataset  values  when  they  are  defined, 


SYSTEM  INITIAL  »  ON 


B.  To  allocated  dataset  values  upon  the  execution  of  the  first 

reference  to  them. 


SYSTEM  INITIAL  =>  OFF 


Comments 


The  default  setting  for  this  parameter  is  ON.  When  on,  all 
datasets  defined  behave  as  though  they  were  fixed  data  sets. 
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3 . 1.15  3"  e  HEADING  Parameter 


This  parameter  allows  the  user  to  tontrol  the  printing  of  page 

headings . 


Svntax 

A.  To  suppress  printing  of  page  headings. 


SYSTEM  HEADING  -  OFF 


B.  To  print  page  headings. 


| SYSTEM  HEADING  ~  ON 


This  is  the  default  setting. 


:-TARI 


me  s..«. 


The  START  statement  starts  a  new  program,  clears  a  working  space 
in  computer  memory-  (core) ,  and  allows  the  user  to  define  a  program  descrip¬ 
tor  and  an  area  descriptor. 

Svntax 


START,  HEADING  (program  descriptor),  AREA  (area  descriptor) 


Comments 


The  START  statement  may  not  be  issued  within  a  model. 
3.2.1  The  HEADING  Option 


The  Program  descriptor  is  an  optional  descriptor  and  it  provides  an  overall 
description  of  the  program  or  a  specific  run.  The  program  descriptor  will 
appear  along  with  the  page  number  at  the  top  right-hand  corner  of  each  page 
of  tabular  program  output. 

Comments 


The  keyword  HEADING  may  be  optionally  omitted. 

3.2.2  The  AREA  Option 

In  regional  analysis  programs,  the  AREA  option  may  be  used  to 
specify  the  name  of  the  area  or  region  being  analyzed.  The  area  descriptor 
will  appear  in  the  title  of  each  page  of  tabular  program  output,  concatenated 
after  the  title. 


i 


STOP 


1-42 


3.3  The  STOP  Statement 


The  STOP  statement  ends  a  program  and  allows  the  user  to  exit 
the  system.  All  information  of  the  program  not  saved  is  lost. 

Svntax 


STOP 


Comments 


Execution  of  a  STOP  statement  within  a  model  causes  both  a 
model  execution  stop  and  an  exit  from  the  system. 


3.4  The  OPEN'  S:a:er?.er.: 

The  OPEN  statement  opens  a  data  file  for  use.  It  r.as  two  forms 
depending  upon  whether  the  file  is  to  be  initialized--! . e .  used  for  the 
first  time  or  whether  the  file  was  initialized  by  a  previous  run  of  the 
system. 


n  is  the  length  of  the  file  control  vector.  If  the  file  is 

to  be  used  as  a  TMS  file  then  this  parameter  equals  the 
number  of  arrays  to  be  placed  on  the  file.  If  the  file 
is  to  be  used  as  a  SAVE  file  then  this  parameter  equals 
four  times  the  number  of  modules  to  be  placed  on  the  file 
plus  one  added  to  the  product. 


size 


is  the  total  size  in  words  of  the  file  to  be  allocated 


Comments 

As  was  discussed  in  Subpart  2.1.9  each  physical  file  may  have  with 
it  a  fixed  maximum  size  which  is  established  either  by  the  operating  system 
directives  used  to  establish  it,  or  which  is  fixed  in  the  system,  or  both. 
The  particulars  depend  upon  the  implementation  of  the  system  at  a  particular 
site.  Before  opening  a  new  physical  file  the  user  should  first  .  _ide  how- 
much  storage  space  he  will  require. 
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Having  made  this  decision  he  should  check  with  his  computer  staff  to  deter 
mine  the  appropriate  operating  system  directives  and/or  program  changes 
needed. 

For  certain  operating  systems  the  opening  of  previously  defined 
files  is  optional.  This  is  not  true  for  all,  however.  Users  at  sites 
where  the  opening  of  old  files  is  optional  are  encouraged  to  use  the  state 
ments  so  as  to  insure  transferability  to  other  operating  systems. 

The  OPEN  statement  may  be  issued  within  or  outside  of  a  model. 


SAVL 


T_„  = 


'.5  The  SAVE  Statement 


The  SAVE  statement  places  the  information  and  instructions  co 
tained  within  a  program  module  or  within  a  segment  of  a  program  module 
on  a  data  base  for  lacer  loading  within  the  same  run  or  within  a  iater 
run.  It  is  through  the  use  of  program  module  segmentation  introduced 
via  the  SAVE  statement  that  the  user  can  define  and  execute  large  model 
within  a  limited  amount  of  core  storage. 


Svntax 


{SAVE  id.  option's"* 

id  is  the  identifier  of  the  program  module  or  program 

module  segment  being  saved 

options  is  a  list  which  may  contain  the  following  two  options: 
XEQ  and  DATA. 


Comments 


The  program  module  information  saved  contains  all  information 
and  instructions  specified  in  all  statements  preceding  the  SAVE  statement 
and  including  either  the  last  START  statement,  if  the  DATA  option  is  not 
used,  or  the  DATA  statement  referenced  by  the  DATA  option.  A  saved  program 
module  is  lost  at  the  end  of  a  batch  run  or  an  on-line  session,  unless  it 
is  saved  permanently  using  the  procedures  of  the  operating  system.  The 
description  of  these  procedures  is  outside  the  scope  of  this  handbook.  , 

The  user  should  consult  his  computer  staff. 

A  saved  module  may  be  used  again  by  using  the  LOAD  statement. 

Normally,  the  saved  module  is  written  on  a  file  whose  physical 
file  number  is  1.  This  file  is  the  save  data  file.  To  change  the  physi- 
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cal  file  number  of  che  file  on  which  modules  are  saved  Che  SYSTEM  SAVE»n 
statement  may  be  used. 

The  SAVE  statement  may  be  used  either  within  a  model  or  outside. 
3.5.1  The  DATA  option 

The  DATA  option  is  used  to  indicate  that  only  a  segment  cf  the 
program  module  is  to  be  saved.  This  is  done  by  specifying  che  dataset 
marking  the  beginning  of  Che  segment  being  saved. 

Svntax 


1  DATA  -  id  | 

id  is  the  identifier  of  the  dataset  whose  information 

represents  the  upper  boundary  of  the  program  module 
segment  being  saved. 

3.5.2  The  XEQ  option 

The  XEQ  option  is  used  to  indicate  that  whenever  the  saved 
program  module  is  loaded,  then  an  indicated  model  within  the  module 
is  to  be  executed  automatically. 

Syntax 


f  XEQ  -  id  j 

id  is  the  identifier  of  the  model  which  is  to  be  executed. 

Comments 


The  keyword  XEQ  may  be  optionally  omitted. 


LOAD 
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3.6  The  LOAD  Statement 

The  LOAD  statement  loads  a  previously  saved  program  module  or 
segment  of  a  program  module. 

Svntax 


LOAD  id,  options 


id  is  the  identifier  of  the  program  module  or  program 

module  segment  being  loaded. 

options  is  a  list  which  may  contain  the  following  two  options: 
XEQ  and  DATA. 


Comments 


The  material  loaded  by  this  statement  may  be  either  an  entire 
program  module  or  a  segment  of  a  program  module.  After  a  LOAD  statement 
for  a  program  module,  the  user  may  continue  interacting  with  the  program 
at  the  same  point  of  interaction  that  he  left  off  when  he  saved  the  pro¬ 
gram.  If  not  preceded  by  a  SAVE  statement,  a  LOAD  statement  permanently 
erases  all  information  currently  in  core  and  replaces  it  with  the  LOAD 
module.  In  this  sense  the  LOAD  statement  behaves  in  the  same  manner  as 
the  STOP  statement. 

In  order  to  be  loaded,  a  module  must  reside  either  on  the  data 
file  whose  physical  unit  number  is  contained  in  the  SYSTEM  SAVE  parameter 
or  whose  physical  unit  number  is  contained  in  the  SYSTEM  LOAD  parameter. 

See  Subparts  3.1.6  and  3.1.7  for  descriptions  of  these  parameters.  If 
a  module  with  the  indicated  identifier  is  present  on  both  data  files  then 
Che  one  from  the  data  file  indicated  by  Che  SAVE  parameter  will  be  loaded. 
If  a  module  with  the  indicated  identifier  is  present  on  neither  data  file, 
Chen  the  system  will  issue  the  message  "UNABLE  TO  LOCATE  REQUESTED  MODULE". 


1 


LOAD 


1-43 


3.6.1  The  DATA  option 

The  DATA  option  is  used  to  indicate  that  only  a  segment  of  a 
program  module  is  to  be  loaded.  This  is  done  by  specifying  the  dataset 
marking  the  beginning  of  the  segment  being  loaded. 

Svntax 


|  DATA  »  id 

id  is  the  identifier  of  the  dataset  marking  the  beginning 

of  the  program  module  segment  being  loaded. 

3.6.2  The  XEQ  option 

The  XEQ  option  is  used  to  indicate  that  a  particular 
model  within  the  program  segment  or  module  being  loaded  is  to  be  executed 
automatically . 

Syntax 

|  XEQ  »  id  1 


id  is  the  identifier  of  the  model  to  be  executed. 

Comments 


If  there  was  an  XEQ  parameter  associated  with  the  SAVE  statement 
which  saved  the  program  module,  the  one  associated  with  the  Load  takes 
dominance . 

The  keyword  XEQ  may  be  optionally  omitted. 


) 


- a  - 


fife 


3.7  The  TIME  Statement 


TIME 


This  command  is  used  in  dynamic  simulation  modeling  to  assign 
values  co  four  key  parameters:  Che  integration  interval  for  time  inte¬ 
grals,  the  beginning  cine  point  or  lower  limit  of  time  integrals,  the 
ending  time  point  or  upper  limit  of  time  integrals ,  and  the  current 
value  for  the  independent  variable  time  itself.  These  four  parameters 
are  stored  by  the  system  in  the  reserved  words  DT,  BEGINNING,  ENDING, 
and  TIME,  respectively. 

Svntax 


TIME (dt .beginning .ending) .option 


dt  is  a  number  greater  than  zero  defining  the  size  of  the  time 

increment  for  subsequent  dynamic  simulations . 

beginning  is  a  number  defining  the  beginning  time  for  subsequent  dynamic 
simulations . 

ending  is  a  number  defining  the  ending  time  for  subsequent  dynamic 
simulations . 

option  is  a  list  that  may  include  the  following  options:  SIZE 

Comments 


r 


The  beginning  and  ending  parameters  define  the  simulation 
period  for  subsequent  runs  of  dynamic  simulation  models.  The  simulation 
wil"  proceed  from  the  beginning  to  the  ending  time  point  in  increments 
of  size  dt. 

Operationally,  the  dt  parameter  is  the  time  interval  used  in 
specifying  the  difference  equations  defining  change  over  time  in  suo- 
sequent  dynamic  models.  Or,  the  dt  parameter  is  the  integration  step 
in  sum  approximations  of  time  integrals  of  continuous  time-dt-punden* 


TIME 


I- 


var_ables.  In  chis  context ,  the  beginning  and  ending  parameters  are 
the  lower  and  upper  limits  of  all  time  integrals  to  be  computed  in 
subsequent  dynamic  models. 

The  constant  values  for  the  dt,  beginning  and  ending  para¬ 
meters  are  stored  in  the  reserved  words  DT,  BEGINNING  and  ENDING, 
respectively.  Also,  the  variable  value  for  time  itself  is  stored 
in  the  word  TIME,  which  initially  is  sec  equal  to  the  value  stored 
in  BEGINNING.  Upon  execution  of  a  dynamic  model,  the  TIME  variable 
takes  the  values 

BEGINNING 

BEGINNING+DT 

BEGINNING+(2*DT) 


BEGINNING+(N*DT) 

where  N  is  the  smallest  integer  to  satisfy  the  following  inequality: 

(N+l) > (ENDING-BEGINNING) /DT+. 5 

The  TIME  command  controls  the  execution  of  subsequent  dynamic 
simulation  models.  Thus,  it  must  be  entered  prior  to  the  execution  of 
any  such  models. 

Before  executing  any  dynamic  simulation  models,  the  DT,  BEGIN¬ 
NING  and  ENDING  parameters  must  have  definite  values.  Also,  since  most 
simulations  move  forward  in  time,  DT  must  have  a  positive  value,  and 
the  ENDING  value  must  be  greater  than  the  BEGINNING  value. 

In  a  program  that  has  a  TIME  command,  all  outputs  are  sub¬ 
scripted  or  concatenated  by  the  current  value  or  range  of  the  TIME 
parameter.  For  fractional  time  values,  use  the  SIZE  option  to  specify 
the  number  of  decimal  places  to  be  shown  with  time  subscripts. 


TIME 
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3 ■ ' . 1  The  SIZE  Option 


The  SIZE  option  is  used  to  specify 
the  TIME  parameter  in  all  outputs  of  dynamic 
the  TIME  parameter  is  displayed  as  a  5-digit 


the  format  for  the  display  of 
simulation  programs.  Normally, 
integer  (with  no  decimal  digits 


Svntax 


SIZE (nc ,nd) 


nc  is  an  integer  specifying  the  total  width  in  characters  in 

all  displays  of  the  TIME  parameter  as  a  subscript  of  simu¬ 
lation  results.  Its  default  value  is  5. 


nd  is  an  integer  specifying  the  number  of  digits  to  the 

right  of  the  decimal  point  in  all  displays  of  the  TIME 
parameter  as  a  subscript  of  simulation  results.  Its 
default  value  is  0. 


ADD 
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3.8  The  ADD  Statement 


The  ADD  statement  adds  table  management  system  arrays  to  data 

files . 

Svntax 

ADD  (db,  table),  options 

db  is  an  integer  value  (l<_dbO)  indicating  which  database 

within  the  database  set  is  to  contain  the  array. 

table  is  an  integer  value  indicating  the  table  number  of  the  array. 

options  is  a  list  that  may  contain  the  following  option:  SIZE 

Comments 


This  statement  places  the  specified  table  on  the  data  file 
whose  physical  unit  number  is  contained  in  the  appropriate  SYSTEM  para¬ 
meter — DB1,  DB2,  or  DB3.  These  are  described  in  Subpart  3.1.4.  The 
values  of  the  array  are  initialized  at  real  zero.  The  table  number 
must  not  exceed  the  length  of  the  file  control  vector  as  specified 
when  the  file  was  opened  for  the  first  time.  See  Section  3.4  for  a 
description  of  the  OPEN  statement. 

Unless  the  size  options  is  used,  the  array  being  added  is 
assumed  to  be  a  scalar — i.e.  it  is  assumed  to  consist  of  a  single 
value . 

The  ADD  statement  may  be  used  either  within  or  outside  of  a 

model. 


The  SIZE  opcion  is  used  Co  indicate  chat  the  array  being  adde 
is  not  scalar.  It  specifies  the  number  of  dimensions  in  the  array  and 
the  number  of  entries  in  each  dimension. 


Svntax 


n 


indicates  the  number  of  dimensions  in  the  array,  (l<n<10) . 


s 


i 


indicates  the  number  of  entries  in,  or  size  of,  the  ith 
dimension,  i=*l,n. 


INDEX 
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3 . 9  The  INDEX  Statement 


The  INDEX  statement  defines  an  index  identifier  and  the  ir.fo- 
maticn  associated  with  it.  3asically,  the  indexes  of  a  program  form,  the 
classification  scheme  for  the  datasets  of  the  program. 


Svntax 


INDEX  id(n),  options 


or 


id 


is  the  index  identifier 


INDEX 

idl  (nl) ,  options 

id2  (n2) ,  options  j 

I 

END 


n 


is  a  positive  integer  defining  the  size  of  the  index 


options  Is  a  list  that  may  include  the  following  options:  LABEL  and 
TIME. 


Comments 


Indexes  are  used  primarily  as  generalized  subscripts  of  data¬ 
sets.  Thus,  index  identifiers  appear  explicitly  in  DATA  statements. 

Ordinarily,  the  index  entries  are  ordered  from  1  to  n  in  the 
following  sequence:  l,2,3,...n. 

The  size  and  relative  ordering  of  the  entries  of  an  index  may¬ 
be  altered  locally  by  the  SET  statement. 

3.9.1  The  LABEL  Option 

This  option  is  used  to  specify  a  descriptor  for  the  index. 


AD-A085  007 

unclassified 

r 1 
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ixqL:; 


Svnca:-: 


LABEL  (label) 


label  Is  a  descriptor  for  the  index  being  defined.  It  may  contain 
up  to  125  characters. 

3.9.2  The  TIME  On  cion 


This  option  is  used  to  specify  a  time  index.  A  time  index  is 
used  in  dynamic  models  to  classify  time-dependent  datasets,  i.e.,  time 
series  variables  whose  values  are  defined  at  the  time  points  specified 
by  the  time  index. 

Syntax 

f  TIME  (list) 


list  is  a  list  of  numbers,  t^.  This  list  may  be  defined  in  two  wavs: 
cl»  c2’  '  ’  '  Cn 

where  n  is  the  size  of  the  index,  or 


in  which  case  n  equally  spaced  time  points  are  generated  implicitly 
by  the  following  formula: 

t.  ,t.  +  h, t.  +  2h,  .  .  .  t 
ill  n 

where  h  -  (t  -  t,)/(n  -  1)  is  the  constant  interval  between  time 
a  1 

points . 
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3.10  The  DATA  Statement 


The  DATA  statement  is  used  to  define  a  dataset.  Basically,  the 
datasets  of  a  program  are  the  storage  arrays  for  the  numeric  information 
that  may  be  manipulated  in  the  equations  of  the  program. 


Syntax 

|  DATA  id(indexes),  options] 


is  the  dataset  identifier 

is  the  list  of  previously  defined  index  identifiers  classifying 
the  dataset.  Up  to  10  different  Indexes  may  be  included  in  this 
list.  Scalar  datasets  do  not  have  any  Indexes. 

is  a  list  that  may  Include  the  following  options:  LABEL,  VALUE, 
and  TMS. 


id 

Indexes 

options 


DATA 

idl  (indexes),  options 
id2  (indexes),  options 


END 


Comments 


If  neither  of  the  options  VALUE  or  TMS  are  used  Chen  the  data¬ 
set  is  assumed  to  be  a  scratch  dataset.  See  Subparc  2.1.3  for  a  discussion 
of  this  notion. 


DATA 


3 . in . 1  The  LABtL  ODCion 
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The  LABEL  option  is  used  to  define  a  label  or  descriptor  for  the 

dataset. 

Svntax 


LABEL  (label) 


label  is  a  descriptor  for  the  dataset. 

Comments 


The  label  of  a  dataset  will  appear  as  the  title  of  all  displays 
of  the  dataset  generated  later  by  SHOW  statements.  Unlabeled  datasets 
are  shown  with  no  titles. 

The  keyword  LABEL  may  be  optionally  omitted. 

3.10.2  The  VALUE  Option 

The  VALUE  option  is  used  to  define  the  numeric  value  or  values 
of  the  dataset. 

Syntax 

VALUE  (values) 


values  is  a  list  of  numeric  values  for  all  storage  cells  of  the  data¬ 
set. 

Comments 


The  list  of  values  for  a  multiply-indexed  dataset  must  be 
entered  in  the  same  order  as  described  explicity  in  the  READ  dataset 
statement,  i.e.,  the  first  index  varies  (increases)  the  fastest,  the 
second  index  varies  the  second  fastest,  and  so  forth. 

When  the  values  of  a  dataset  are  specified  via  a  VALUE  option, 
then  that  dataset  becomes  fixed.  Subpart  2.1.3  discusses  the  notion 
of  fixed  versus  scratch  datasets. 


DATA 
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3.10.3  The  TMS  Option 

The  TMS  option  specifies  chat  the  dataset  is  a  table  management 
system,  or  TMS,  dataset.  As  is  discussed  in  Subparc  2.1.3  TMS  datasets 
have  in  essence  two  locations  —  one  within  a  TMS  array  on  some  data  file 
and  one  in  a  temporary  location  within  the  program  module.  The  purpose 
of  the  TMS  option  is  to  specify  within  which  TMS  array  the  dataset  values 
are  located,  where  within  the  TMS  array  the  dataset  values  begin,  how 
the  dataset  values  are  oriented  within  the  TMS  array,  and  the  logical 
relationship  between  the  dataset  and  the  TMS  array.  Prior  to  discussing 
the  syntax  of  the  TMS  option  itself,  these  concepts  are  discussed  in 
general. 

TMS  arrays  are  discussed  in  general  in  Subparts  2.1.1  and 
2.1.11  and  in  Part  3.8.  This  discussion  is  not  repeated  here.  Basi¬ 
cally,  to  identify  a  TMS  array  two  integer  values  must  be  included  — 
the  number  of  the  data  base  containing  the  array  and  the  table  number 
of  the  array.  These  values  are  always  included  within  the  TMS  option. 

The  specification  of  the  starting  location  of  the  dataset 

value  within  the  TMS  array  is  done  via  a  base  point  specification. 

Let  T  be  a  TMS  array  with  size  specification  (s,  . .  s  )  as  specified 

l  n 

in  Part  3.7.  Now  any  point  P  in  T  can  be  specified  via  an  ordered 
sequence  (e^,...,ea),  where  e ^  is  the  logical  entry  number  of  the 
point  in  the  ith  dimension.  Obviously  lj<e^si  for  all  dimensions  i. 

Now  a  subarray  A  of  T  can  be  defined  via  the  point  in  T  at  which  the 
(l^,...,ln)  point  of  A  occurs  and  the  size  of  A.  The  point  in  T  at 
which  the  (l^,...,!^)  point  of  A  occurs  is  called  the  base  point  of 
A  in  T. 

The  specification  of  the  orientation  of  the  dataset  within  the 
TMS  array  hinges  upon  the  notion  of  transposition.  Assume  A  is  an  N  X  M 
array  and  B  is  the  transpose  of  A.  Then  B  is  an  M  x  N  array.  Further¬ 
more,  it  can  be  said  that  the  rows  of  B  correspond  to  the  columns  of  A, 
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and  Che  columns  of  B  correspond  co  Che  rows  of  A.  Ic  is  Chis  ncCion  of 
correspondence  in  che  cranspose  chat  is  used  co  define  che  oriencacicn 
of  che  dacasee  within  che  T>IS  array. 

Svncax 

| plS  (db,  cable,  opeions) 


db  Chis  inceger  value  (l<db<_3)  indicaces  which  dacabase  wichin 

che  dacabase  sec  concains  che  array 

cable  chis  inceger  value  indicates  che  cable  number  of  the  arrav 

options  includes  any  number  of  the  suboptions  from  che  following 

_list:  BASE.  TRANSPOSE. 

3.10.3.1  The  BASE  Suboption 

The  BASE  suboption  is  used  to  define  the  base  point  of  the 
dataset  within  the  TMS  array. 


DATA 
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Syntax 


BASE  (val, .... ,val  ) 

1 _ n 

n  is  the  number  of  dimensions  in  the  TMS  array  identified  by 

db  and  table 

val^  is  the  ith  coordinate  for  the  base  point.  It  may  either  be 

an  integer  constant,  a  dataset  identifier,  or  one  of  the 
system  parameters  associated  with  time  —  DT ,  TIME, 
BEGINNING  or  ENDING. 

Comments 

If  the  BASE  suboption  is  omitted,  then  the  system  assumes  a 

base  point  specified  as  follows:  BASE  (1,,..,1  ). 

■L  n 

3.1Q.3.2  The  TRANSPOSE  Suboption 

The  TRANSPOSE  suboption  is  used  to  define  the  transposition 
correspondence  between  the  dataset  and  the  TMS  array. 

Syntax 


TRANSPOSE  (ind  . ,ind  ) 

_  A  . n .. 


n  is  the  number  of  dimensions  in  the  TMS  array  identified  by 

db  and  table 

indA  is  either  one  of  the  index  identifiers  defining  the  dataset  or 

an  asterisk.  If  it  is  an  index  identifier,  then  the  dataset 
values  classified  by  that  index  are  stored  along  the  ith  dimen¬ 
sion  oi  the  TMS  array.  7f  Lz  i~  *n  asterisk,  t'ren  the  d'fset 
has  no  values  stored  along  it.',  dimension  of  the  TMS  -r  ay. 
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Comments 

If  the  TRANSPOSE  suboption  is  omitted,  then  the  system  assumes 
a  one-to-one  correspondence  between  the  dataset  and  the  TMS  array.  This 
correspondence  would  be  specified  as  follows:  TRANSPOSE  (indexes),  where 
indexes  is  the  same  list  as  appeared  in  the  dataset  definition. 

If  the  indexes  in  the  TRANSPOSE  suboption  are  few'er  than  the 
dimensions  in  the  TMS  array,  then  the  remaining  dimensions  are  invariant, 
i.e.,  not  transposed,  relative  to  the  dataset. 


t 


sSUMM 


DATA 
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3.10.4  Indirect  Data  Sets 


Another  variation  of  the  DATA  command  allows  the  user  to  define 
an  "indirect"  dataset.  This  may  indirectly  represent  one  of  a  list  of 
other  datasets.  Indexes  may  not  be  defined  with  indirects.  Indirects 
may  not  be  used  in  equations.  (See  the  ELSE  DATA  option  of  the  ASK 
statement  for  examples  of  use.) 

General  Syntax 


DATA  id* 


id  is  the  identifier  of  the  indirect  dataset.  An  indirect 

identifier  may  be  used  as  a  substitute,  or  indirect 
representation,  for  one  of  a  list  of  "direct"  datasets. 


CRA:’ii 
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3.11  The  GRAPH  Statement 


The  GPAPH  statement  is  used  Co  define  a  graph.  Basically,  a 
graph  is  a  sec  of  ordered  pairs  of  numbers  chac  specify  and  x-  and 
y-coordinaces  of  che  points  defining  che  graph. 


Svncax 


|  GRAPH  id  (n) ,  options  ] 


i-  is  che  graph  identifier 

n  is  a  positive  integer  defining  the  size  of  the  graph,  i.e., 

the  number  of  points  defining  the  graph 

options  is  a  list  that  may  include  the  following  options:  X,Y. 

Comments 


or 


GRAPH 

idl  (nl) ,  options 
id!  (n2) ,  options 
!  cND 


Graphs  are  used  in  arithmetic  expressions  in  the  RHS  of  equations 
to  yield  by  linear  interpolation  the  y-value  (or  y-values)  corresponding 
to  some  argument  (or  arguments). 

The  argument  of  a  graph  may  be  an  arithmetic  expression,  in 

general. 


is 


The  syntax  for  the  general  use  of  a  graph  in  arithmetic  expressions 


id  (arithmetic  expres s ion ) 


where  "id"  is  the  graph  identifier  and  an  "arithmetic  expression"  is  an 
expression  of  numbers,  datasets,  graphs,  and  arithmetic  operators. 


GRAPH 
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The  value  of  a  graph  for  an  arbitrary  argument  is  obtained  by 
linear  interpolation  or  extrapolation  between  the  points  defining  the 
graph . 

If  not  specified  by  the  X  and  Y  options,  the  values  of  a  graph 
are  specified  or  modified  by  the  READ  statement. 

Graphs  may  be  used  as  arguments  of  other  graphs. 

The  argument  of  a  graph  need  not  be  a  scalar,  i.e.,  it  may  be 
an  array  of  values  resulting  from  an  indexed  dataset  or  an  arithmetic 
expression  involving  such  datasets. 

3.11.1  The  X  Option 

The  X  option  is  used  to  define  the  x-coordinates  of  the  n  points 
defining  the  graph. 

Syntax 


X(xlt  x2,  .  .  .  ,  x  ) 


X-,  x_,  .  .  ,x  are  the  numeric  values  for  the  x-coordinates  of  the  n 
i  i  n 

points  defining  the  graph. 

3.11.2  The  Y  Option 

The  Y  option  is  used  to  define  the  y-coordinates  of  the  n  points 
defining  the  graph. 

Syntax 


Y(y,»  y2» 


V 


y^i  y2.  •  •  •  yR  are  the  numeric  values  for  the  y-coordinates  of  the  n 
points  defining  the  graph. 


GRAPH 


3.11.3  Examples  of  GRAPH  Statements 

1.  The  Absolute  Value  Graph 

GRAPH  AB-S ( 3 )  ,  X(-1,0,1),  Y(1,0,1) 


A  schematic  represencacion  of  this  graph  is: 


For  an  arbitrary  argument  x  the  value  of  ABS(x)  yields: 

f  x  if  x  >  0 
ABS(x)  *  \ 

-x  if  x  <  0 

2.  The  Step  Function  Graph 

GRAPH  STEP (4) ,  X(-10, -.001,0, 10) ,  Y(0, 0,1,1) 

A  schematic  representation  of  this  graph  is: 


1 

y-STEP(x)J 

‘(0,1) 

,(10.: 

L _ 

(-10,0)  (-.001,0) 

X 

For  an  arbitrary  argument  x,  the  value  of  STEP(x)  yields: 


STEP(x) 


I  o  if  *  <  0 

)  1  if  x  >  0 


The  discontinuity  at  x  ■  0  is  treated  approximately  to  with 


GRAPH 
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3.  A  Graph  of  a  Graph 

Consider  the  following  expression: 

ABS ( (STEP (-20) ) 

What  is  its  value? 

From  Example  2,  the  value  of  STEP(-20)  is  0.  From  Example  1, 
the  value  of  ABS(O)  is  0.  Thus,  the  value  of  ABS (STEP (-20) )  is  also  0. 
Similarly,  the  value  of  STEP (ABS (10) )  is  equal  to  1. 

4.  A  Graph  with  a  Multidimensional  Argument 


Consider  Che  dataset  DR  classified  by  an  index  whose  size 
is  3.  If  the  3  values  of  this  dataset  are:  DR(l)«*-3,  DR(2)*2,  DR(3)“10, 
then  the  expression  ABS (DR)  takes  the  following  3  values:  3,  2,  and  10. 


MODEL 
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3.12  The  MODEL  Statement 


.he  MODEL  statement  is  used  to  introduce  a  model  identifier  and 
the  equations  arc  instructions  associated  with  it  for  later  execution. 

ivr.tax 


id  ts  tr.e  model  identifier. 

The  model  instructions  are  various  statements  that  define  the 
equations  and  other  instructions  of  the  model.  These  instructions  are 
only  interpreted  and  stored.  They  are  not  executed. 

The  END  statement  terminates  the  compilation  of  the  model  and 
returns  the  user  to  command  mode. 

Comments 


The  MODEL  statement  initiates  compile  mode .  Under  compile  mode, 
all  model  instructions  are  interpreted  and  stored,  but  not  executed.  The 
compiled  model  instructions  are  not  executed  until  the  model  is  called 
into  execution. 

Model  execution  is  initiated  by  entering  the  model  identifier. 
Model  execution  proceeds  sequentially  through  the  compiled  instructions 
of  the  model  in  the  order  in  which  they  are  compiled. 

If  the  AREA  option  is  specified  in  the  START  statement,  a  model 
execution  statement  may  contain  an  area  descriptor  characterizing  the 
region  associated  with  this  execution  of  the  model. 

There  are  two  types  of  models  —  static  and  dynamic.  Both 
use  the  same  syntax  when  they  are  initiated.  Dynamic  models  contain 
within  them  the  RATE  and  LEVEL  statements.  These  are  discussed  in  Parts 
3.22  and  3.23. 


I 


MODEL 
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3.12.1  Examples  of  MODEL  Statements 
1.  Model  Compilation 
INDEX  ROW( 3) 

DATA  DR  (ROW) ,  LABEL  (DATA  BY  ROW) 

DATA  TOT,  LABEL  (TOTAL) 

MODEL  SUMR 
READ  DR 

TOT-SUM(R)  (DR(R)) 

SHOW  DR. 2 
SHOW  TOT. 2 
END 

Upon  model  execution  (see  Step  1  below)  the  model  SUMR  will: 

•  Prompt  the  user  to  enter  the  3  ROW  values  for  the 
dataset  DR  (see  Step  2  below) . 

•  Compute  the  value  of  the  dataset  TOT  by  summing 
over  the  3  ROW  values  of  the  dataset  DR 

«  Show  the  values  of  the  dataset  DR,  to  2  decimal 
places  (see  Step  3  below). 

•  Show  the  value  of  the  dataset  TOT,  to  2  decimal 
places  (see  Step  4  below). 
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2.  Model  Execution 
Step 

1  SUMR 

2  DATA  123 

3 

R0W1 
ROW  3 

4  TOTAL  6.00 


DATA  BY  ROW 
1.00  ROW 2 
3.00 


TABLE 
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3.13  The  TABLE  Statement 

The  TABLE  scatemenc  defines  for  later  execution  a  table  identi¬ 
fier,  the  datasets  associated  with  it,  and  its  format. 

To  call  a  table  into  execution,  i.e.,  to  print  it  out  or  to 
display  it,  simply  enter  its  identifier*. 

Svntax 


TABLE  id (indexes) ,  DATA(datasets) ,  TITLE(title) ,  SIZE(sw,cv)  { 


id 


is  the  table  identifier 


indexes  is  a  list  of  previously  defined  index  identifiers.  These  indexes 
define  the  spanner  and  subtitles  of  the  table.  The  order  of 
indexes  in  the  list  is  important  (see  the  comments). 

datasets  is  a  list  of  previously  defined  dataset  identifiers.  This  list 

may  also  contain  slashes,  /  ,  and  strings  enclosed  in  parentheses. 
Upon  table  execution,  the  strings  are  displayed  as  part  of  the 
table  stub  in  the  same  relative  order  as  they  were  specified  in 
this  list.  The  slashes  manifest  themselves  as  blank  lines  or 
rows  in  the  table  in  the  same  relative  order  as  they  were  spec: fied 
in  this  list.  Upon  execution  the  values  of  the  datasets  included 
in  this  list  become  the  table  body. 

title  is  a  descriptor  which  upon  execution  becomes  the  centered  title 
of  the  table 


sw 


is  a  positive  integer  defining  the  table  stub  width  (in 
characters) 


cw 


is  a  positive  integer  defining  the  column  widths  (in 
characters)  for  the  table. 


TABLE 
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Comments 


A  table  execution  statement  may  be  entered  either  in  command 
or  in  compile  (i.e.,  within  a  model)  mode. 

The  order  of  display  of  the  table  indexes  is:  the  missing  index 
classifies  the  table  rows,  the  last  index  classifies  the  table  columns. 

All  other  indexes  classify  the  various  row  by  column  table  sections  in 
a  lef t-to-right  order,  i.e.,  the  first  index  varies  the  fastest,  the 
second  index  varies  the  second  fastest,  etc. 

For  datasets  not  classified  by  the  missing  index,  the  dataset 


label  acts  as  a  stub. 


TABU 
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3.13.1  Examples  of  TABLE  Statements 


1.  Table  Definition 


INDEX  ROW (3) 

INDEX  COL (2) 

DATA  DRC(ROW.COL) ,  LABEL (DATA  BY  ROW  AND  COLUMN) 

DATA  TOTC(COL),  LABEL (TOTALS  BY  COLUMN) 

TABLE  STAB(COL),  DATA (TOTC • 2 , / ,  DRC .  2 ,  /  ,  (NOTE.  THESE  ARE  DUMMY  VALUES)), 
SIZE(20,10) ,  TITLE (A  SAMPLE  TABLE) 

2.  Table  Execution 

DRC*  2 

TOTC (C)=SUM(R)  (DRC(R, C)  ) 

STAB 


A  SAMPLE  PROGRAM  PAGE  1 

A  SAMPLE  TABLE 


COL 

COL 

ONE 

TWO 

TOTALS  BY  COLUMN 

6.00 

6.00 

ROW1 

2.00 

2.00 

ROW2 

2.00 

2.00 

ROW3 

2.00 

2.00 

NOTE.  THESE  ARE  DUMMY  VALUES 


END 


3.14  The  END  Statement 

The  END  statement  is  used  to  terminate  the  compilation  of  a 
model  or  the  section-form  definition  of  a  set  of  indexes,  or  a  set  of 
datasets,  or  a  set  of  graphs,  or  a  set  of  tables.  It  is  also  used  to 
end  the  range  of  a  conditional  statement  such  as  IF  or  DO. 


Svntax 


To  distinguish  between  END  statements  terminating  various 
structured  sections  of  code,  arbitrary  comments  may  be  used  following 
the  word  END  and  at  least  a  blank  or  a  comma. 


READ 
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3.15  The  READ  Statement 


The  READ  statement  is  the  basic  input  operation  of  the  system 
and  is  used  to  read  information  into  the  structure  of  indexes,  datasets, 
or  graphs. 

Svntax 


READ  id,  options 
{information  being  read  in} 


id  is  the  identifier  of  the  index,  dataset,  or  graph  where  the 

information  being  read  in  will  be  stored 

options  these  options  vary  depending  on  whether  information  is  read 
into  an  index,  a  dataset,  or  a  graph.  See  the  READ  index, 
READ  dataset,  and  READ  graph  statements. 

Comments 


To  satisfy  the  READ  statement,  information  must  be  entered  in 
the  cards  or  input  lines  following  the  READ  card. 

The  type  of  information  entered  depends  on  whether  it  is  read 
into  an  index,  dataset,  or  graph. 

The  amount  of  information  entered  depends  on  the  local  dimensions 
of  the  structural  element  identified  in  the  READ  statement. 

The  manner  in  which  the  information  is  entered  on  the  input  cards 
or  lines  depends  both  on  the  implicit  structure  of  the  structural  element 
identified  in  the  READ  statement  and  an  explicit  format  specified  in  the 
options  of  the  READ  statement. 

READ  statements  may  be  specified  either  in  command  or  in  compile 
mode  (except  READ  index).  In  command  mode,  the  information  being  read  in 
must  follow  immediately  the  READ  statement.  In  compile  mode,  the  infor¬ 
mation  being  read  in  must  not  follow  the  READ  statement  which  is  specified 
within  some  model;  rather,  it  must  be  entered  after  a  prompt  issued  by  the 
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model  execution  statement. 

3.15.1  The  READ  Index  Statement 


The  READ  index  statement  is  used  to  read  information  into  the 
structure  of  an  index.  With  this  statement  the  user  nay  introduce  the 
following  information: 

•  The  stub  of  an  index 

•  The  spanner  of  an  index 


Svntax 


READ  id,  options 
{index  entries} 


id  is  an  index  identifier 

options  may  include  format  specifications  for  either  the  stub  or  the 
spanner  entries  of  the  index. 

Comments 


The  index  entries  following  the  READ  statement  must  be  specified 
on  n  entry  cards  or  input  lines,  where  n  is  the  size  of  the  index. 

The  manner  in  which  the  index  entries  are  keypunched  on  the 
cards  or  typed  in  the  input  lines  depends  on  the  index  option  being 
specified. 

The  READ  index  statement  may  be  issued  only  outside  of  a  model, 
and  is  not  sensitive  to  the  SET  command. 

3.15.1.1  The  Stub  Option 

The  stub  option  is  used  with  the  READ  index  statement  to  read  a 
stub  into  the  structure  of  an  index.  The  index  stub  is  the  set  of  descrip¬ 
tors  that  will  become  the  row  descriptors  for  the  displays  of  all  datasets 
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that  are  classified  by  the  index. 

Svntax 


(ic.lc) 

{n  stub  entries; 


ic  is  a  positive  integer  (<80)  defining  the  initial  column  on  each 

stub  card  or  line  where  the  stub  entry  begins 

lc  is  a  positive  integer  (lc  >_  ic)  defining  the  last  column  on  each 

stub  card  or  line  where  the  stub  entry  ends. 

Comments 

Each  of  the  n  stub  cards  or  lines  following  the  READ  statement 
must  contain  a  stub  entry  for  the  index  identified  in  the  READ  index 
statement.  Each  stub  entry  must  be  located  between  the  ic  and  lc  columns 
of  each  card. 

The  number  of  stub  entries  is  n,  where  n  is  the  size  of  the 
index  identified  in  the  READ  index  statement.  These  entries  are  stored 
in  the  order:  1,  2,  3,  .  .  .  ,  n. 

Both  format  parameters  ic  and  lc  may  be  omitted,  in  which  case 
the  default  values  ic  *  1  and  lc  *  20  apply. 

3.15.1.2  The  Spanner  Option 

The  spanner  option  is  used  with  the  READ  index  statement  to  read 
a  spanner  into  the  structure  of  an  index.  The  index  spanner  is  the  set 
of  descriptors  that  will  become  the  column  headings  for  the  tabular  dis¬ 
plays  of  all  datasets  that  are  classified  by  the  index. 
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Svntax 


(ic , lc ,nc) 

{n  spanner  entries} 


ic  is  a  positive  integer  (<80)  defining  the  initial  column  on  each 

spanner  card  or  line  where  the  spanner  entry  begins 

lc  is  a  positive  integer  (lc  _>  ic)  defining  the  last  column  on 

each  spanner  card  or  line  where  the  spanner  entry  ends 

nc  is  a  positive  integer  defining  the  number  of  columns  or  characters 

(including  blanks)  in  each  line  section  of  the  spanner.  The 
following  exact  relationship  must  be  satisfied: 

(lc  -  ic  +•  l)/nc  -  nl 

where  nl  is  a  positive  integer  denoting  the  number  of  lines  in 
each  column  heading  of  the  spanner. 

Comments 


Each  of  the  n  spanner  cards  or  lines  following  the  READ  state¬ 
ment  must  contain  a  spanner  entry  for  the  index  identified  in  the  READ 
statement.  Each  spanner  entry  must  be  located  between  the  ic  and  lc 
columns  of  each  card  or  line,  forming  nl  groups  of  nc  characters  each. 

The  number  of  spanner  entries  is  n,  where  n  is  the  size  of 
the  index  identified  in  the  READ  index  statement.  These  entries  are 
stored  in  the  order:  1,  2,  3,  .  .  .  ,  n. 
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The  format  parameters  ic,  1c ,  and  nc  are  necessary  for  specifying 
the  line  structure  of  each  spanner  and  thus  may  not  be  omitted. 


3.15.2  The  READ  Dataset  Statement 


The  READ  dataset  statement  is  used  to  introduce  values  into  the 
structure  of  a  dataset. 

Svntax 


READ(indexes)  id,  (ic,lc,nc) 
{data  entries} _ 


id  is  a  dataset  identifier 


indexes  is  an  ordered  list  of  the  indexes,  minus  one ,  classifying  the 
dataset.  The  missing  index  plays  an  important  role  (see  the 
comments).  The  ordering  of  the  indexes  in  the  list  is  also 
important  (see  the  comments).  By  definition,  this  list  is  not 
necessary  for  scalar  datasets.  This  list  may  be  omitted,  in 
which  case  the  default  ordering  of  indexes  applies  (see  comments). 

ic  is  a  positive  integer  (<80)  indicating  the  initial  column  on 

each  data  card  or  line  where  the  reading  operation  will  start 

lc  is  a  positive  integer  (lc  _>  ic)  indicating  the  last  column  on 

each  data  card  or  line  where  the  reading  operation  will  end 


nc  is  a  positive  integer  indicating  the  number  of  columns  or 

characters  forming  each  data  entry  on  each  data  card  or  line. 
The  following  exact  relationship  must  be  satisified: 


(lc  -  ic  +  l)/nc  ■  ne 
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where  ne  is  a  positive  integer  denoting  the  maximum  number  of 
entries  that  may  be  entered  on  each  data  card. 


Comments 

The  fixed  format  parameters  ic,  lc,  nc  may  be  omitted,  in  which 
case  the  values  of  the  dataset  may  be  entered  in  free  format.  In  free 
format,  the  values  of  a  dataset  may  be  entered  anywhere  on  an  input  card 
or  line  provided  they  are  separated  from  each  other  by  at  least  one  blank 
or  a  comma. 


3.15.2.1  Examples  of  READ  Dataset  Statements 

The  reading  operation  specified  by  a  READ  dataset  statement  can 
best  be  demonstrated  by  means  of  examples. 

Consider  the  dataset  A(I,J,K)  classified  by  the  indexes  I,  J, 
and  K  whose  sizes  are  3,2,  and  2,  respectively.  I  is  the  first,  J  is  the 
second,  and  K  is  the  third  index  of  the  dataset  A. 


Example  A.  To  satisfy  the  statement 


READ  A 


the  values  of  A  must  be  arranged  on  the  data  entry 


cards  as  follows : 


Card 

1 

2 

3 

4 


A(l,l,l) 

A(l,2,l) 

A(l,l,2) 

A(l,2,2) 


Values 

A(2,l,l) 

A(2,2,l) 

A(2,l,2) 

A(2,2,2) 


A(3,l,l) 
A(3, 2,1) 
A(3,l,2) 
A(3,2,2) 


In  this  arrangement,  the  first  index  of  A,  I,  varies  the  fastest  across 
each  card,  the  second  index,  J,  varies  the  second  fastest,  and  the  third 
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Example  3.  To  satisfy  the  statement 
READ(I,J)  A 
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the  values  of  A  must  be  arranged  as  follows: 


Card 

Values 

1 

A(l, 1 , 1) 

A(l,l,2) 

A(2,l,l) 

A(2,l,2) 

3 

A(3,l,l) 

A(3,l,21 

4 

A(l, 2 , 1) 

A(l, 2, 2) 

5 

A(2,2,l) 

A(2, 2,2) 

6 

A(3,2,l) 

A(3, 2,2) 

In  this  arrangement,  the  missing  index  K  varies  the  fastest  across  each 
card,  the  index  I  varies  the  second  fastest,  and  the  index  J  varies  the 
third  fastest. 

Example  C.  To  satisfy  the  statement 
READ(I,J,K)  A 

the  values  of  A  must  be  arranged  on  the  data  entry  cards  as  follows: 


Card  Values 

1  A(l,l,l) 

2  A(2,l,l) 

3  A(3,l,l) 

4  A(l, 2, 1) 

5  A(2,2,l) 

6  A(3,2,l) 
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Card 

Values 

7 

A(l,l,2) 

8 

A  ( 2 , 1 , 2 ) 

9 

A(3, 1, 2) 

10 

A(l,2,2) 

11 

A(2,2,2) 

12 

A(3 , 2 , 2) 

Free  versus  Fixed  Format.  The  format  parameters  ic,  lc,  nc  are 
optional.  If  specified,  the  dataset  values  must  be  entered  according  to 
the  fixed  format  specified  by  these  parameters.  If  not  specified,  the 
dataset  values  may  be  entered  in  free  format ,  provided  they  are  separated 
by  commas  or  blanks. 

Examples  A,  B,  and  C  are  in  free  format.  The  following  example 
is  in  fixed  format. 

Example  D.  To  satisfy  the  statement 

READ  A,  (11,40,10) 

the  values  of  A  must  be  arranged  on  the  data  entry  cards  in  the  same  order 
as  shown  in  Example  A.  However,  the  entries  on  each  card  must  lie  between 
columns  11  and  40,  and  each  entry  must  not  exceed  10  columns  or  characters. 
Columns  1  to  10  are  ignored  by  the  reading  operation.  However,  they  may 
be  used  for  data  identification  purposes. 

In  this  example  the  format  parameters  ic,  lc,  and  nc  have  the 

values : 


ic  -  11 
lc  -  40 
nc  •  10 


Thus,  the  number  of  entries  per  card  is: 
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ne  =  (lc  -  ic  +  l)/nc  *  (40  -  11  -  D/10  »  3. 

Note.  A  READ  dataset  operation  is  not  satisfied  unless  all  values  of  the 
dataset  are  entered  appropriately. 

Normally,  a  dataset  contains  n  values,  where  n  is  equal  to  the 
product  of  the  sizes  of  the  indexes  classifying  the  dataset.  Thus,  the 
sample  dataset  A  contains  3x2x2=12  values,  where  3,  2,  and  2  are  the  sizes 
of  the  indexes,  I,  J,  and  K. 

However,  the  size  and  relative  ordering  of  an  index  may  be 
modified  by  a  SET  index  statement.  Thus,  the  range  and  structure  of  a 
READ  dataset  operation  will  be  modified  by  those  SET  index  statements  that 
modify  the  indexes  of  the  dataset.  This  is  illustrated  in  the  next  example. 

Example  E.  To  satisfy  the  following  READ  statement: 

SET  1(1,3),  J (2) 

READ  A 

the  values  for  A  must  be  arranged  as  follows: 

Card  _ Values _ 

1  A(l,2,l)  A(3,2,l) 

2  A(l,2,2)  A(3,2,l) 

3.15.3  The  READ  Graph  Statement 


The  READ  graph  statement  is  used  to  introduce  values  into  the 
structure  of  a  graph. 

Syntax 


|  READ  id,  options 
j'n  x-values 
:  n  y-values. 


READ 


1-82 


IS  2 


zzaplx  identifier 


n  is  the  graph  size 

options  include  the  optional  specification  of  fixed  format  parameters. 


Comments 


To  satisfy  the  READ  graph  statement  n  x-values  and  n  y-values 
must  be  entered,  where  n  is  the  size  of  the  graph. 

The  n  x-values  must  be  entered  first  on  as  many  cards  or  input 
lines  as  are  necessary,  followed  by  the  n  y-values. 

Since  the  x-  and  y-values  form  pairs  of  coordinates  on  the  x-v 
plane  for  the  n  points  of  the  graph,  the  n  y-values  must  be  entered  in  the 
same  order  as  the  corresponding  n  x-values . 

In  free  format  the  x-  and  y-values  may  be  entered  in  the  desired 
order  separated  by  commas  or  blanks. 

The  fixed  format  option  may  be  specified  as  follows : 

(ic,lc,nc) 


where 

ic  is  a  positive  integer  (<80)  denoting  the  initial  column  on  each 

data  card  or  line  where  the  reading  operation  will  start 

lc  i_.  a  positive  integer  (lc  >_  ic)  denoting  the  last  column  on 

each  data  card  or  line  where  the  reading  operation  will  end 

nc  is  a  positive  integer  (<80)  indicating  the  number  of  columns 

or  characters  forming  each  x-  or  y-values  on  each  data  card 
or  line.  The  following  exact  relationship  must  be  satisfied: 

(lc  -  ic  +  1) /nc  •  ne 

where  ne  is  a  positive  integer  denoting  the  maximum  number  of 
entries  that  may  be  entered  on  each  data  card. 
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3.15.3.1  Examples  of  READ  Graph  Statements 

1.  The  Absolute  Value  Graph  Specified  in  Free  Format 

GRAPH  ABS(3) 

READ  ABS 
-10  1 
10  1 


2.  The  Step  Function  Graph  Specified  in  Fixed  Formac 

GRAPH  STEP (4) 

READ  STEP,  (6,25,5) 

STEPX  -10  <-.001  0  10 

STEPY  0  011 

Since  the  first  5  columns  of  each  data  card  are  ignored  by  the  READ 
operation,  they  may  be  used  for  identification.  In  this  example, 

ic  *  6 
lc  -  25 
nc  *  5 

and  the  number  of  entries  per  card,  ne,  is: 

ne  ■  (lc  -  ic  +  l)/nc  »  (25  -  6  +  l)/5  *  4, 


the  number  of  points  defining  the  graph  STEP. 
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3.16  The  SET  Statement 


The  SET  statement  is  used  to  set  locally  the  range  of  an  index 
and  the  relative  ordering  of  its  entries. 

The  range  of  an  index  may  be  a  positive  integer  less  than  or 
equal  to  n,  the  size  of  the  index. 

The  normal  ordering  of  the  entries  of  an  index  is:  1,  2,  3, 

.  .  .  ,  n.  The  SET  index  statement  may  modify  this  ordering. 


Syntax 


SET  id,  setting 


id  is  an  index  identifier 

setting  is  the  setting  for  the  index  entries.  This  setting  may  be 

specified  in  two  different  ways : 

(value  list) 

* 

The  *  setting  restores  the  index  to  its  normal  setting  as 
established  in  the  index  definition.  -  i.e.,  range  -  n 
and  entry  ordering:  1,  2,  .  .  .  ,  n 

The  value  list  may  include  up  to  n  positive  values  or  as 
many  nonoverlapping  value  intervals  as  are  necessary  to  cover 
the  set  (1,  2,  .  .  .  ,  n>.  A  value  interval  is  specified  by 
the  two  values  defining  the  interval  separated  by  a  hyphen, 
The  two  values  of  an  interval  may  be  specified  either  in 
ascending  or  in  descending  order. 
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Comments 


A  single  SET  statement  may  be  used  to  set  more  than  one  index., 
as  follows: 


SET  id^(setting  ) ,  id,,  (setting,,) ,  .... 

A 

The  setting  of  an  index  stays  in  effect  until  the  index  is  reset 
by  another  SET  statement  and  it  affects  the  following  instructions:  READ 
dataset,  SHOW  dataset,  and  subscripted  equations. 

Note .  The  SET  index  statement  is  very  powerful  because  it 
affects  subsequent  operations.  Thus,  it  must  be  used  with  caution  and 
only  locally,  i.e.,  once  the  operations  which  should  be  done  under  the 
influence  of  a  SET  statement  are  completed  the  indexes  affected  should 
be  restored  to  their  normal  setting  by  using  the  *  option. 

3.16.1  Examples  of  the  SET  Statement 

INDEX 

R0W(3) 

COL ( 2 ) 

TYPE(5) 

END 


Normally,  the  ROW  index  has  size  3  and  its  entries  are  ordered 
as  follows : 

Entry  Order 

ROW(l)  1 

R0W(2)  2 

R0W(3)  3 

Normally,  the  COL  index  has  size  2  and  its  entries  are  ordered 


as  follows : 


SET 


2.  Setting  Two  Indexes 


SET  ROW(3-l) ,  COL (2) 


SHOW 


3.17  The  SHOW  Statement 
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The  SHOW  statement  is  the  basic  output  instruction  of  the 
language  and  it  is  used  to  show  the  values  of  a  dataset. 

Syntax 


| 

j  SHOW(sw.cw)  id.n,  option 


id  is  a  dataset  identifier 

sw  is  a  positive  integer  denoting  the  stub  width  in  char¬ 

acters  for  the  row  descriptors  of  the  dataset  display 

cw  is  a  positive  integer  denoting  the  column  width  in  char¬ 

acters  for  the  column  headings  of  the  dataset  display 

n  is  a  positive  integer  (<cw)  denoting  the  number  of  decimal 

digits  to  be  shown  with  each  dataset  value 

option  is  one  of  the  optional  parameters  ORDER,  TITLE,  and  TOTAL. 

Comments 

The  stub  width,  sw,  and  column  width,  cw,  parameters  may  be 
omitted.  If  omitted,  the  default  values  in  effect  are:  sw  -  IS 
cw  •  8. 

The  decimal  point  and  the  decimal  digits  parameter,  n,  may 
be  omitted.  If  omitted,  the  default  value  in  effect  is:  n  -  0.  The 
dataset  values  are  shown  as  integers. 

A  scalar  dataset  is  shown  on  a  one-line  display  consisting  of 
the  label  of  the  dataset,  if  specified,  followed  by  the  value  of  the 
dataset. 
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An  indexed  dacasec  is  shown  on  a  centered,  paged,  and  titled 
tabular  display  that  is  properly  described  by  the  entries  of  the  indexes 
classifying  the  dataset.  The  page  heading  is  shown  at  the  cop  right- 
hand  corner  of  the  display  and  it  consists  of  the  program  descriptor 
specified  in  the  START  statement  followed  by  the  page  number.  The 
title  is  the  label  of  the  dataset.  The  row  descriptors  of  the  display 
are  the  stub  entries  of  Che  first  index  of  the  dataset.  The  column 
headings  of  the  display  are  the  spanner  entries  of  the  second  index  of 
the  dataset.  The  stub  entries  for  the  third,  fourth,  etc.,  indexes 
of  the  dataset  are  the  subtitles  of  the  row-by-column  sections  of  the 
display. 

The  tabular  display  of  the  dataset  is  designed  to  fit  properly 
within  an  8-1/2  x  11  inch  page  for  direct  insertion  in  reports. 

Large  dataset  displays  are  segmented  into  sections,  each 
section  fitting  properly  within  an  8-1/2  x  11  inch  page. 

Dataset  values  greater  than  or  equal  to  1,000  are  shown  with 
commas  separating  the  value  into  3-digit  groups. 

3.17,1  The  ORDER  Option 

The  ORDER  option  allows  the  user  to  control  the  order  in 
which  the  display  is  organized — i.e.,  which  index  is  placed  in  the 
rows,  which  in  the  columns,  and  so  forth. 

Syntax 

\ 

ORDER (indexes) 


indexes  is  a  list  of  Indexes  specifying  the  order  of  the  indexes  in 
the  display.  The  first  index  classifies  the  rows  and  the 
second  the  columns  of  the  dataset  being  shown.  All  ether  indexes 
classify  as  subtitles  the  row-by-column  sections  of  the  data¬ 
set  being  shown. 
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3.17.2  The  TITLE  Option 

The  TITLE  option  is  used  to  specify  4  title  for  the  dataset 
being  shown.  If  the  TITLE  option  is  not  used,  the  title  for  the  tabular 
display  is  the  label  of  the  dataset  being  shown . 

Svnta>: 


TITLE  (list) 


list  is  the  title  to  be  shown  for  the  dataset.  This  list  may 

contain  one  or  more  of  the  following  items : 

(1)  a  string  of  up  to  125  characters,  and/or 

(2)  the  special  character  /,  to  skip  to  a 
new  line  for  multiple  line  titles,  and/or 

(3)  id.L,  where  the  notation  "id.L"  refers  to 
the  label  of  the  dataset  whose  identifier 
is  "id". 

If  Item  1  is  not  listed  alone,  then  it  must  be  enclosed  in 
parentheses. 

The  items  In  a  multi  item  list  must  be  separated  by  blanks 

or  commas. 

3.17.3  The  TOTAL  Option 

The  TOTAL  option  allows  the  user  to  specify  that  all  (or  a  sub¬ 
set)  of  the  values  of  a  dataset  be  totaled  and  the  totals  (or  subtotals) 
be  shown  in  the  display  of  the  dataset  values. 


TOTAL (indexes) 


indexes  is  a  list  of  identifiers  specifying  those  indexes  that 
classify  those  values  of  the  dataset  which  are  to  be 
totaled.  If  omitted,  the  values  of  the  dataset  will 
be  totaled  and  subtotaled  over  each  and  all  indexes 
classifying  the  dataset.  If  included,  only  those 
values  will  be  subtotaled  that  are  classified  by  the 
indexes  specified. 


Equacicn- 
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3.18  Equations 


This  is  a  continuation  of  the  discussion  on  equations  given  in 

Part  2.1. 

Syntax 


iid^arithmetic  expression! 

id  is  the  identifier  of  a  previously-defined  dataset  (explicitly 

subscripted  or  not) . 

Arithmetic  expression  is  a  well-formed,  ordered  sequence  of 
numerical  constants,  previously-defined  identifiers,  mathematical  oper¬ 
ators,  and  parentheses. 

Comments 

The  following  are  additional  comments  to  those  given  in  Part 
2.1.  These  comments  address  the  following  topics: 

•  lef t-to-right  hierarchy 

•  dummy  or  symbolic  subscripts 

•  Implicit  subscripts 

•  matrix  equations 

•  division  by  zero 

•  linear  interpolation 

Lef t-to-Right  Hierarchy.  In  the  absence  of  parentheses,  the 
arithmetic  operations  on  the  right-hand  side  of  equations  are  compiled 
and  executed  from  left  to  right. 
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All  arithmetic  operators  +,-,*,/,  and  **  for  addition,  sub¬ 
traction,  multiplication,  division,  and  exponentiation,  respectively, 
are  of  equivalent  rank  in  the  processing  of  equations  by  the  system. 

In  other  computer  languages  and  in  standard  algebraic  notation  this 
is  not  the  case.  In  FORTRAN,  for  example,  there  are  three  levels  of 
hierarchy:  exponentiation  is  processed  first,  then  division  and  muitipl 
cation,  and  then  addition  and  subtraction. 

Hierarchy  rules  for  processing  arithmetic  expressions  are  very 
important.  Without  a  fixed  set  of  rules,  expressions  such  as  7*4+3/2 
could  have  several  meanings  and  values : 

(7*4)+(3/2)=2S  •  5  _ 
or  ( ( 7*4)+3) / 2= 15 . 5 

or  7*(4+3)/2=24. 5 
or  7*(4+(3/2) ) *38 . 5 

Which  is  the  correct  answer?  Algebraic  hierarchy  rules  state 
that  multiplication  and  division  must  be  performed  before  addition. 

Thus,  algebraically  the  correct  answer  is  (7*4)  +  (3/2)=29 . 5 . 

In  the  system  hierarchy,  however,  the  arithmetic  operations 
are  performed  from  lef t-to-right ,  regardless  of  operation  type.  Thus, 
the  arithmetic  expression  7*4+3/2  is  processed  as  follows: 

•  The  multiplication  7*4  is  done  first  yielding  28 

•  The  addition  28+3  is  done  next  yielding  31 

•  The  division  31/2  is  done  next  yielding  15.5 

The  operations  in  an  arithmetic  expression  are  performed 
from  left  to  right  with  each  operator  using  the  preceding  result  as 
the  operand.  Even  though  this  convention  is  different  that  the  stand¬ 
ard  FORTRAN  convention,  it  may  be  easier  to  remember,  especially  by 
beginner  programmers,  and  it  has  the  advantage  chat  the  actual  program¬ 
ming  code  required  to  process  it  is  shorter  and  faster. 


F  i  u  a  t  i  <'  n  s 


To  alter  the  itr.plied  lef t-to-right  hierarchy  in  which  arith¬ 
metic  operations  are  performed,  the  user  may  emplov  parentheses.  The 
meaning  and  use  of  parenthetical  groups  in  arithmetic  expressions  is 
discussed  in  Part  2.1. 

Dummy  Subscripts.  The  following  single  equation: 

DRC  (R,  C)  =■  DR(R)  *  DC(C) 

defined  in  terms  of  the  following  indexes  and  datasets: 

INDEX  ROW (3) 

INDEX  C0L(2) 

DATA  DRC (ROW, COL) 

DATA  DR (ROW) 

DATA  DC (COL) 


is  equivalent  to  the  following  set 


of  6  single-valued  equations: 


DRC (1, 1)  -  DR(1)  *  DC( 1) 
DRC (2,1)  -  DR(2)  *  DC (1) 
DRC (3,1)  -  DR(3)  *  DC(1) 
DRC (1,2)  *  DR(1)  *  DC(2) 
DRC (2, 2)  *  DR (2)  *  DC(2) 
DRC (3, 2)  -  DR (3)  *  DC(2) 


Here,  the  dummy  subscripts  R  and  C  -  symbolic  for  the  indexes  ROW  and 
COL,  respectively  -  take  on  the  values: 


R:  1,2,3 
C:  1,2 

as  dictated  by  the  respective  sizes  of  the  ROW  and  COL  indexes. 

Any  TIE  of  type  2  can  be  used  as  a  dummy  subscript,  provided 
it  has  not  been  used  previously  as  an  identifier.  Usually,  single¬ 
character  type  2  TIEs  make  convenient  dunum-  subscripts. 
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Actual  index  identifiers  may  serve  as  their  own  dummy  sun- 
scripts,  of  course. 

As  can  be  seen  from  the  example  above,  dummy  subscripting 
represents  substantial  simplification  in  notation.  The  power  of  this 
notation  may  be  appreciated  more  fully  when  compared  with  its  FORTRAN 
equivalent,  the  DO  Loop. 

Explicit  subscripts,  i.e. ,  integer  constants  or  dataset 
identifiers,  must  have  nonzero,  positive  values  never  exceeding  the 
sizes  of  the  indexes  which  they  represent. 

A  dummy  subscript  in  the  left-hand  side  of  an  equation  takes 
on  all  the  values  allowed  by  the  size  -  actually,  the  current  range  as 
specified  by  the  last  SET  index  statement  -  of  the  index  it  represents. 

For  a  multiply-subscripted  dataset  the  order  of  the  dummy 
subscripts  must  correspond  to  the  order  of  the  indexes  classifying  the 
dataset.  For  example,  the  dummy  subscripts  I , J , K  in  the  following  data¬ 
set  reference 


DRTC(I,J,K) 

correspond  to  the  indexes  ROW,  TYPE,  and  COL  respectively  that  classify  the 
dataset  DRTC,  as  specified  in  the  dataset  definition: 

.  DATA  DRTC (ROW,  TYPE,  COL) 

An  equation  with  dummy  subscripts  corresponds  to  a  set  of  inde¬ 
pendent,  single-valued  equations.  The  number  of  single-valued  aquations 
represented  by  such  an  equation  is  equal  to  the  product  of  the  current 
ranges  of  the  dummy  subscripts  classifying  the  dataset  in  the  left-hand 
side  of  the  equation.  The  dummy  subscripts  of  the  right-hand  side  of 
an  equation  take  on  the  same  values  as  the  subscripts  with  identical 
symbols  in  the  left-hand  side. 
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Implicit  Subscripts.  In  the  equation 

DR-2 

che  index  ROW  classifying  che  dacasec  DR  is  implicit,  i.e.,  not  specified 
explicitly  as  a  subscript  of  DR.  This  equation  is  equivalent  to  the 
following  3  equations  -  one  for  each  of  the  3  entries  of  the  ROW  incex 
classifying  che  DR  dataset: 


DR( 1)  -  2 
DR( 2)  -  2 
DR(3)  -  2 

This  example  brings  out  the  notational  economy  of  implicit 
subscripting,  which  can  be  considerable  in  equations  involving  multi¬ 
dimensional  datasets  classified  with  large  indexes. 

‘  Matrix  Equations.  Equations  with  dummy  or  implicit  subscripts 
may  be  referred  to  as  matrix  equations.  A  single  matrix  equation  corre¬ 
sponds  to  a  number  of  independent,  single-valued  equations.  The  number 
of  single-valued  equations  represented  by  a  matrix  equation  is  equal  to 
the  product  of  the  current  ranges  of  the  dummy  or  implicit  subscripts 
classifying  the  dataset  in  the  left-hand  side  of  the  equation. 

The  dummy  or  implicit  subscripts  of  the  right-hand  side  of 
a  matrix  equation  take  on  the  same  values  as  the  corresponding  sub¬ 
scripts  of  the  left-hand  side.  A  redundant  subscript  in  the  right-hand 
side,  i.e.,  a  subscript  with  no  corresponding  subscript  in  the  left- 
hand  side,  is  set  at  the  first  sequence  value  in  the  current  setting  of 
the  index  corresponding  to  the  subscript.  For  example,  the  notation 


SET  K(2,l),  L* 

A  (I,  J)  -  B  (I, J ,K,L) 
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is  equivalent  to 


A(  I ,  J)  -  3  (.  I ,  J  ,  2 , 1 ) 

where  2  is  the  first  sequence  number  in  the  current  setting  of  the  index 
K,  and  1  is  the  first  sequence  number  in  the  current  setting  of  the  index 
represented  by  the  dummy  subscript  L.  However  the  notation 

Set  K(?,l) ,  L* 

A  =  B 


is  equivalent  to 


A(I , J)  =  B ( I , J , 1 , 1) 

Dummy  and  implicit  subscripts  must  not  be  mixed.  When  using 
dummy  or  explicit  subscripts,  all  subscripts  of  all  datasets  in  the 
equation  must  be  specified  explicitly  and  in  the  appropriate  order. 

When  using  implicit  subscripts,  all  subscripts  of  all  datasets  in  the 
equation  must  be  implicit. 

If  not  subscripted  with  explicit  or  dummy  indexes,  the  datasets 
in  the  right-hand  side  of  a  matrix  equation  must  be  classified  by  implicit 
indexes  that  match  exactly  the  implicit  indexes  classifying  the  dataset 
in  the  left-hand  side  of  the  equation. 

For  example,  the  equation 

A  =  3 

is  malformed  if  the  dataset  A  is  indexed  by  the  ROW  (size  3)  index  and  the 
dataset  B  is  indexed  by  the  COL  (size  2)  index.  This  can  be  seen  from  the 
following  single-cell  equations  chat  are  equivalent  to  A  *  B : 
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A(l)  -  B(l) 

A(2)  -  B (2) 

A(3)  -  B(3) 

However,  whac  is  che  value  of  B(3)?  The  subscripc  3  is  out  of  range  for 
the  dataset  B.  The  dataset  B  has  only  two  values,  B(l)  and  B(2),  as 
dictated  by  the  size  of  the  COL  index  that  classifies  it.  In  this  case, 
the  system  uses  the  first  value  of  the  dataset  B,  i.e.,  it  interpets  the 
third  equation  above  as  follows : 

A(3)  -  3(1) 

This  is  true  in  general:  a  subscript  out  of  range  in  the 
right-hand  side  of  an  equation  is  set  to  one. 

Division  by  Zero.  Since  all  datasets  are  initialized  at  zero, 
it  is  possible  to  accidentally  introduce  divisions  by  zero  in  the  execution 
of  a  program.  Divisions  by  zero  are  computed  as  zero;  not  infinity.  A 
division  by  zero  causes  neither  a  program  execution  error  nor  a  program 
stop. 


Linear  Interpolation.  Linear  interpolation  is  used  implicitly 
by  the  system  to  evaluate  the  values  of  graphs  for  arbitrary  arguments. 

In  general,  the  method  of  linear  interpolation  is  a  method  for 
approximating  the  value  of  a  function  at  some  arbitrary  point  from  the 
known  values  of  the  function  of  a  finite  set  of  points. 

Consider,  for  example, the  4-point  function  f(x)  depicted  below, 
whose  values  are  given  only  at  four  points: 
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These  four  points  are  defined  by  the  pairs  of  numbers 
(x1>f1) , U2,f2> , (x3>f 3> ,  and  (x^.f,). 


Mint  is  the  value  of  this  function  at  an  arbitrary'  point  x?  If  it  is 
assumed  that  the  four  points  are  connected  by  straight  lines  in  the 
manner  shown: 


then  the  value  of  the  function  at  an  arbitrary  point  x^  is  given  by  the 

point  of  intersection  of  a  vertical  line  drawn  at  x^  and  the  line  connecting 

the  two  known  points  of  the  function  that  are  nearest  to  x  ,  on  either  side. 

P 

In  terms  of  the  diagram  above,  the  value  f  of  the  function  at 
Xp  is  given  by  the  linear  interpolation  formula: 
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f 

P 


E2  +  ^f3  *  f2^*^xp  "  x,)/(x3  -  x2)) 


where  (x^.f^)  and  (x^.f^)  are  the  two  points  that  are  nearest  to  (x 
either  side. 

Similarly,  the  value  f  of  the  function  at  x^  is  given  by 
linear  extrapolation  formula: 


_»f  )  on 

P  P 
the 


f 

q 


f3  +  ( ( f 4  -  f3)*(xq  -  x3)/(x4  -  x3)) 


which  is  given  in  terms  of  the  points  at  x3  ^nd  x^,  which  are  nearest  to  x 
but  lie  to  the  left  of  it. 

In  general,  the  value  of  the  n-point  function,  f (x) ,  at  an 
arbitrary  point  x  is  given  by  the  general  linear  interpolation  formula: 


q 


f(x)  -  f ±-h( Cf i+1  -  fi)*(x  -  *i)/(*1+1  ~  xi)) 

where  x^  and  x^+^  are  those  x  coordinates  which  are  nearest  to  x  on  either 
side  and  and  are  the  corresponding  function  values. 

If  x  is  less  than  xi  or  greater  than  xq,  then  the  following 
linear  extrapolation  formulae  apply: 


f(x  <  xx) 

f(x  >  x  ) 
n 


fl  +  ^f2  ”  “  xi^/^x2  “xi^ 


f  1  +  C(f  “f  J*(x  -  x  )/(x  -  x  .)) 
n-i  n  n-i  n-1  n  n-1 
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3.19  Model  Execution  Statement 


The  model  execution  statement  executes  the  statements  in  oraer 
within  a  previously  defined  user  modal.  It  can  be  used  to  execute  a 
static  model  or  a  dynamic  model.  It  can  also  be  used  to  resume  the 
execution  of  a  dynamic  model. 

Syntax 

A.  Executing  a  model, 


model  I 


B.  Resuming  a  dynamic  model, 


model  »  ending 


model  is  the  identifier  of  a  previously  defined  model 

ending  is  the  new  value  of  ending  time  until  which  the  model 

execution  should  continue. 


Comments 


See  Part  3.7  on  the  TIME  statement  and  Part  3.12  on  the  MODE 
statement  for  descriptions  of  the  manner  in  which  models  are  executed. 
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3.20  Table  Execution  Statement 

The  table  execution  statement  generates  a  tabular  report  as 
specified  by  the  identified  table. 


Svntax 

table 

table  is  the  identifier  of  the  table  specifying  the  tabular  report. 


Comments 

See  Part  3.13  on  the  TABLE  statement  for  details  of  specifying 
tabular  reports. 


RATE 
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3.21  The  RATE  Statement 


The  RATE  statement  is  used  in  the  dynamic  simulation  models  and 
it  has  two  functions :  (a)  it  signals  the  end  of  the  initial  section  of 

the  model  and  the  beginning  of  the  rate  section,  and  (b)  it  declares  that 
the  time  dependent  variables  be  computed  at  each  time  point  of  the  simu¬ 
lation  by  linear  interpolation  or  extrapolation  from  specified  exogenous 
time  series. 

Svntax 


RATE (exogenous  variables  list) 


where  the  exogenous  variables  list  is  a  list  of  correspondences  between 
previously-defined  exogenous  time  series  and  time-dependent  variables 
that  must  be  used  locally  in  the  rate  (and/or  level)  section  of  a  dy¬ 
namic  simulation  model.  This  list  must  be  enclosed  in  parentheses, 
and  each  list  item  must  be  formulated  as  follows: 


exogenous  time  series=local  variable 


where  exogenous  time  series  is  a  previously-defined  dataset  that  is  in¬ 
dexed  by  a  time  index,  and  local  variable  is  a  previously-defined  data 
set  that  is  used  explicitly  in  the  rate  (and/or  level)  section  of  a 
dynamic  model  and  is  computed  at  every  time  point  of  the  simulation. 
Based  on  this  equivalence,  the  values  of  the  local  variable  will  be 
computed  at  the  arbitrary  time  points  of  the  dynamic  simulation  by 
linear  interpolation  or  extrapolation  that  is  based  on  the  fixed  time 
poincs  defining  the  exogenous  time  series. 

Comments 


The  RATE  and  LEVEL  statements  are  mandatory  in  dynamic 
simulation  models.  Each  dynamic  simulation  model  must  have  one  RATE  and 
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one  LEVEL  statement  in  order  to  separate  the  model  equations  into  tr.ree 
sections:  The  initial  section,  the  rate  section,  and  the  level  section. 

The  initial  section  is  the  first  section  of  the  model  and  its 
equations  are  evaluated  at  the  beginning  time  point  (or  interval)  of 
the  simulation  period. 

The  rate  section  is  the  second  section  of  the  model  and  its 
equations  are  evaluated  at  each  time  point  (or  interval)  of  the  simu¬ 
lation  run.  In  contrast  to  level  equations,  both  sides  of  rate  equations 
are  evaluated  at  the  same  time  point  (or  interval) . 

The  level  section  follows  the  rate  section  and  its  equations 
are  also  evaluated  at  each  time  point  (or  interval)  of  the  simulation. 

The  left-hand  side  of  each  level  equation,  however,  is  evaluated  at 
TIME+DT  in  terms  of  the  time  variables  in  the  right-hand  side  which 
are  evaluated  at  TIME,  the  previous  time  point  (or  interval).  It  is 
the  equations  of  the  level  section  which  move  the  dynamic  variables 
through  time. 

Only  those  exogenous  time  series  that  are  used  explicitly 
in  the  rate  or  the  level  section  need  be  included  in  the  exogenous 
variables  list  of  the  RATE  instruction. 

An  exogenous  time  series  must  not  be  used  explicitly  in  the 
equations  of  the  rate  or  level  sections ,  for  its  fixed  time  index  points 
must  be  distinguished  from  the  arbitrary  time  points  of  the  TIME  para¬ 
meter  of  the  simulation  run. 

The  RATE  statement  may  only  be  used  Inside  the  MODEL  state¬ 
ment,  i.e.,  you  must  not  use  it  while  in  command  mode. 
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3.22  The  LEVEL  Statement 


The  LEVEL  scacement  is  used  in  dynamic  s imulat ion  models,  and 
ic  has  two  functions:  (a)  it  signals  the  end  of  the  rate  section  and 
the  beginning  of  the  level  section,  and  (b)  it  declares  that  endogenous 
variables  be  computed  and  stored  at  the  fixed  time  points  of  the  time 
indexes  classifying  the  output  variables  time  series.  The  values  of  an 
output  time  series  at  each  time  point  of  the  time  index  are  set  equal 
to  the  values  of  the  local  endogenous  variable  corresponding  to  the 
nearest  simulation  time  point  plus  or  minus  DT/2. 

Svntax 


LEVEL (endogenous  variables  list) 


where  the  endogenous  variables  list  is  a  list  of  correspondences  between 
output  time  series  and  endogenous  variables  that  are  used  locally  in  <_he 
equations  of  the  level  (and/or  rate)  section  of  a  dynamic  simulation 
model.  This  list  must  be  enclosed  in  parentheses,  and  each  list  item 
must  be  formulated  as  follows: 


output  time  series*endogenous  variable 


where  output  time  series  is  a  previously-defined  data  set  that  is  indexed 
by  a  time  index,  and  endogenous  variable  is  a  previously-defined  data 
set  that  is  used  explicitly  in  the  level  (and/or  rate)  section  of  a  dynamic 
simulation  model.  Based  on  this  equivalence,  the  values  of  the  output 
time  series  will  be  computed  and  stored  at  the  fixed  time  points  of  the 
time  index  classifying  the  series.  The  value  of  an  output  series  at  a 
time  point  T  is  set  equal  to  the  computed  value  of  the  corresponding 
endogenous  variable  that  is  associated  with  the  interval  (T-DT/2. 

T+DT/2).  This  interval  is  closed  at  T-DT/2  and  open  at  T+DT/2:  if  T 
is  the  exact  midpoint  of  the  interval,  then  the  T-DT/2  value  applies. 


4 
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Commencs 


The  LEVEL  statement  is  mandatory  in  each  dynamic  simulation  model, 
for  it  separates  the  rate  from  the  level  section  of  the  model,  and  it 
"moves"  the  model  variables  through  time. 

The  LEVEL  statement  causes  the  TIME  parameter  to  be  incremented 
by  DT  units  from  ics  value  in  the  preceding  rate  section. 

The  variable  in  the  left-hand  side  of  each  level  equation  is 
evaluated  at  TIME+DT,  the  incremented  time  point  or  interval,  in  terms 
of  the  variables  in  the  right-hand  side  of  the  equation  which  are  eval¬ 
uated  at  TIME,  the  previous  time  point  or  interval. 

The  output  time  series  in  the  endogenous  variables  list  of 
the  LEVEL  statement  must  not  be  referenced  explicitly  in  the  equations 
of  a  dynamic  simulation  model. 

Only  those  endogenous  variables  that  you  intend  to  save  for 
later  use  as  time  series  need  to  be  included  in  the  endogenous  variables 
list  of  the  LEVEL  statement. 

The  LEVEL  statement  must  not  be  used  in  command  mode,  i.e.,  it 
may  only  be  used  inside  dynamic  simulation  models . 
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3.23  The  WRITE  Statement 


Though  che  reading  and  writing  of  TMS  array  values  is  normally 
taken  care  of  automatically  via  the  definition  and  manipulation  of  TMS 
datasets,  there  are  occasionally  times  when  the  user  might  simply  wish 
the  treat  the  TMS  array  data  file  as  though  it  were  a  standard  database. 
The  WRITE  statement  allows  the  user  to  explicitly  write  dataset  values 
to  a  TMS  array . 

Syntax 

WRITE  data,  TMS  (pram)‘~j 


data  is  the  identifier  of  the  dataset  containing  the  values 

to  be  written 

pram  contains  a  set  of  TMS  parameters  formed  exactly  as  described 

in  Subpart  3.10.3. 


Comments 


The  keyword  TMS  may  be  optionally  omitted. 


3.24  The  INPUT  Statement 


The  INPUT  statement  allow  the  user  to  explicitly  input  datase 
values  from  a  TMS  array. 

Svntax 


INPUT  data,  TMS  (pram) 


data  is  the  identifier  of  the  data  set  containing  the  values  of 

Che  dataset  to  receive  the  TMS  array  information 

pram  contains  a  set  of  TMS  parameters  formed  exactly  as  described 

in  Subpart  3.10.3* 

Comments 


The  keyword  TMS  may  be  optionally  omitted. 
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3.25  The  PUT  Statement 


The  PUT  statement  is  used  to  control  the  corresponcence  for  7>1S 
dataset  values  between  working  storage  and  off  line  storage. 

Svntax 

A.  To  puc  off  line  values  into  working  storage 


PUT  C0R£  (dataset  list) 

B.  To  put  working  storage  values  into  off  line  storage 


PUT  IMS  (dataset  list) 


dataset  list  contains  the  datasets  whose  values  are  to  be  moved  from  one 
medium  to  the  other. 


Comments 


The  keyword  CORE  may  be  omitted. 
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3.26  The  DROP  Statement 


The  DROP  statement  allows  the  user  to  drop  the  space  associated 
with  nonfixed  datasets  from  working  storage;  thus,  making  that  space 
available  for  other  purposes. 

Svntax 


A.  To  drop  some  values 

DROP (list ) 

list  is  a  list  of  nonfixed  datasets  whose  values  are  to  be 

dropped. 

B.  To  drop  all  values 


DROP* 


3.27  The  IF  Sea tement 


It  is  often  desirable  to  execute  a  certain  section  of  a  model 
only  if  some  particular  condition  is  met.  The  IF  statement  allows  this 
operation. 

Svntax 


Boolean  is  a  Boolean  expression  as  described  below. 

Comments 

The  statements  between  the  IF  and  END  are  executed  only  IF 
the  Boolean  is  true. 
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3.27.1  Boolean  Expressions  and  Variables 


Examples 


EXSIM  LT  1970 

ENSIM  GT  1985 

Boolean  expressions  are  expressions  formulated  with  one-dimen¬ 
sional  data  sets  and  constants  and  relational  and  Boolean  operators .  A 
Boolean  expression  may  be  either  true  or  false.  If  it  is  true,  it  has 
the  value  1.  If  it  is  false,  it  has  the  value  0.  As  such.  Boolean 
expressions  may  be  thought  of  as  variables  that  may  take  the  values  1  or  0. 

The  Boolean  variable  in  the  first  example  has  the  value  1, 
i.e.,  it  is  true,  if  the  value  of  the  dataset  ENSIM  is  less  than  1970. 

If  the  value  of  ENSIM  is  greater  than  or  equal  to  1970,  then  this 
Boolean  variable  is  false  and  it  has  the  value  zero. 

The  Boolean  expression  in  the  second  example  is  true,  if  the 


value  of  ENSIM 

is  greater  than 

1985;  otherwise,  it  is  false. 

The 

system  allows  the 

:  following  six  relational  operators: 

Operator 

Meaning 

LT 

Less  than 

LE 

Less  than  or 

equal  to 

GT 

Greater  than 

GE 

Greater  than 

or  equal  to 

NE 

Not  equal  to 

EQ 

Equal  to 

Relational  operators  are  used  to  formulate  relational  expressions 
between  one  dataset  and  another  or  one  dataset  and  a  constant. 

In  addition,  the  system  allows  the  following  three  Boolean 


operators : 


Ooerator 


Meaning 

and 


AND 

OR  or 

NOT  not 

Boolean  operators  may  be  use’d  to  link  Boolean  variables.  Boolean 
expressions,  therefore,  are  either  single  3oolean  variables  or  strings  of 
such  variables  linked  with  Boolean  operators. 

In  constructing  Boolean  expressions  you  may  link  as  many  3oolean 
variables  as  you  wish. 

Also,  Boolean  expressions  may  be  nested,  the  degree  of  nesting 
being  unlimited,  in  principle.  When  formulating  complicated  nested 
Boolean  expressions,  it  is  advisable  that  you  enclose  in  parentheses  the 
component  Boolean  expressions  being  nested.  As  in  arithmetic  operations, 
nested  parenthetical  groups  of  Boolean  operations  are  executed  in  the 
following  order:  the  most  deeply  nested  Boolean  expression  is  executed 
first.  As  in  every  other  statement,  parentheses  in  Boolean  expressions 
must  be  balanced. 

Boolean  expressions  are  used  primarily  in  IF  statements  where 
they  define  the  various  alternative  conditions  that  must  be  satisfied 
for  execution  to  take  alternative  paths. 

Boolean  variables  may  be  used  in  arithmetic  expressions  where 
they  take  only  two  values:  1  or  0.  A  Boolean  variable  has  the  value 
1  if  it  is  true,  and  it  takes  the  value  0  if  it  is  false. 


DO 
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3 • 28  The  DO-  Sta c ament 


It  is  often  desirable  to  perform  repeatedly  a  given  set  of 
instructions  until  a  given  condition  is  met  or  while  a  given  condition 
is  met.  The  DO  statement  allows  these  operations. 

Syntax 

A.  Do  until  a  condition  is  met, 


DO  UNTIL  Boolean 
{statements ; 

END 


B.  Do  while  a  condition  is  met, 


DO  WHILE  Boolean 
{statements } 

END 


Comments 


See  Subpart  3.27.1  for  a  discussion  of  Boolean  expressions. 


SHOW  PLOTS 
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3.29  The  SHOW  Statement  -  with  XYPLOT.  SCATTER,  or  BAR 


This  statement  permits  the  user  to  plot  either  the  values  or  a 
time  series  dataset  versus  time,  or  the  values  of  one  vector  dataset 
(the  dependent  variable  or  y-coordinate)  versus  the  values  of  another 
vector  dataset  (the  independent  variable  or  x-coordir.ate)  .  In  the  latter 
case,  the  x-coordinates  for  the  points  of  the  plot  correspond  to  the 
values  of  the  time  index  classifying  the  values  of  the  time  series  data¬ 
set  being  depicted.  XYPLOT  displays  a  continuous  line  of  asterisks 
(default)  interpolated  between  values,  whereas  SCATTER  displays  only  a 
single  asterisk  for  each  value  SHOWn  in  the  relationship.  BAR  displays 
a  bar  graph.  Examples  of  each  type  of  plot  are  given  in  3.29.1. 

Syntax 


('SCATTER' 

SHOW, < XYPLOT  S  data  , 

options 

LBAR  J 

(-SCATTER- 

SHOW .^XYPLOT  ,  (data^,  data^)  , 

options 

Lbar  j 

data 


is  che  time-series  dataset  to  be  plotted  versus  time 


data^ 


data^ 


option 


is  the  independent  values  (i.e.,  X-axis  variable);  the 
minimum  and  maximum  values  of  data^  will  appear  as 
numeric  captions  underneath  the  display 
is  the  dependent  values  (ie.  ,  Y-axis  variable);  minimum 
and  maximum  values  of  data^  will  appear  as  numeric 
captions  down  the  left  side  of  the  display 
is  one  of  the  optional  parameters  POINT ,  TITLE,  XLaBEL, 
XRANGE,  Y LAB EL,  Y RANGE  and  SUBTITLE. 


'  SHOW  PLOTS 
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Comments 


The  first  syntax  must  contain  a  dataset  that  is  defined  on^a 
TIME  index.  The  TIME  series  information  (the  independent  variable)  will 
be  displayed  for  each  TIME  point  specified  in  the  TIME  index. 

In  the  second  syntax,  data^  is  the  independent  variable  and  data^ 
the  dependent  variable.  Data^  and  data?  may  be  multiple-dimensioned  data 
sets  defined  providing  they  have  at  least  one  index  that  is  common  to  both 
and  that  the  extra  indices  are  SET  to  the  desired  entries  for  the  display. 

3.29.1  The  POINT  Option 

The  POINT  option  allows  the  user  to  select  the  character  used 
in  the  plot. 

Svntax 


POINT (string) 


string  is  a  single  character  to  be  used  in  the  plots.  Default 
is  an  asterisk. 

3.29.2  The  TITLE  Option 

The  TITLE  option  allows  the  user  to  specify  the  title  for  the 
plot.  The  default  title  for  non-TIME  series  datasets  is  blank  (i.e.,  no 
title) ;  for  a  TIME  series  dataset  the  default  TITLE  is  the  label  of  the 
dataset  followed  by  the  TIME  range  (e.g.,  1970  to  1980). 


SHOW  PLOTS 
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Syntax 


TITLE  (list) 


list  consists  of  one  or  store  of  the  following  in  any  order 

cc  -  carriage  control  characters  (e.g.,/  to  skip  to  the 
next  line) . 

string  -  a  string  must  be  enclosed  within  a  second  set  of 
parenthesis  unless  the  entire  list  consists  only 
of  the  string. 

id.  L  -  id  is  a  dataset  identifier;  id.  L  will  print  the 
label  of  the  dataset  id. 

3.29.3  The  XLABEL  Option 

The  XLABEL  option  allows  the  user  to  supply  a  descriptor  for  the 
X-axis;  default  is  no  descriptor. 

Syntax 


XLABEL  (string) 


string  is  a  string  of  characters  to  appear  below  the  X-axis. 

3.29.4  The  XRANGE  Option 

The  XRANGE  option  is  used  only  with  XYPLOT  and  SCATTER.  It  is 
used  to  supply  the  minimum  and  maximum  values  for  the  X-axis. 

Svntax 

XRANGE  (min,  max ) 


min 


max 


is  the  minimum  value  for  the  X-axis, 
is  the  maximum  value  for  the  X-ari - . 


SHOW  PLO‘>’ 
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Comments 


The  default  value  used  for  min  is  the  actual  minimum  value  of 
the  independent  variable.  The  default  value  used  for  max  is  the  actual 
maximum  value  of  the  independent  variable.  If  the  actual  values  of  the 
independent  variable  lie  outside  the  (min,  max)  range  supplied,  then  the 
resulting  plot  may  be  affected  adversely. 

3.29.5  The  YLABEL  Option 

The  YLABEL  option  allows  the  user  to  supply  a  descriptor  for 
the  Y-axis;  default  is  no  descriptor. 

Syntax 


YLABEL  (string) 


string  is  a  string  of  characters  to  appear  to  the  left  of  the  Y-axis 
in  a  column.  . 

3.29.6  The  YRANGE  Option 

The  YRANGE  option  is  used  to  supply  the  minimum  and  maximum 
values  for  the  Y-axis. 

Syntax 

YRANGE  (min,  max) 

min  is  the  minimum  value  for  the  Y-axis. 

is  the  maximum  value  for  the  Y-axis. 


max 


SHOW  PLOTS 


Comments 


The  default  value  used  for  min  is  the  actual  minimum  value  of 
the  dependent  variable.  The  default  value  used  for  max  is  the  actual 
maximum  value  of  the  dependent  variable.  If  the  actual  values  of  the 
dependent  variable  lie  outside  the  (min,  max)  range  supplied,  then  the 
resulting  plot  may  be  affected  adversely. 

3.29. 7  The  SUBTITLE  Ootion 


Syntax 


SUBTITLE  (string) 


string  is  a  string  of  characters  to  be  positioned  below  the  TITLE 
of  the  display. 


Comments 


The  default  SUBTITLE  is  the  fixed  index  settings  for  multi- 
dimensioned  datasets;  if  singly-dimensioned  vectors  are  used,  the 
default  SUBTITLE  is  blank. 

3.29.9  Examples  of  SHOW  with  SCATTER,  XYPLOT  and  BAR 

1.  If  A  is  a  time  series  dataset  defined  yearly  from  1970  to  1975  then 
SHOW  SCATTER  A 
produces  the  following  plot 


SHOV:  PLOTS 
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DATASET  A,  1970  TO  1975 
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2.  If  A  and  B  are  both  time-series  datasets  defined  yearly  from  1970  to 
1975  then 

SHOW  XYPL0T(A,B) ,XLABEL (DATASET  A) ,YLABEL (DATASET  B) , 

TITLE (DATASET  A  VERSUS  DATASET  B) 
produces  the  following  plot. 


SHOW  PLOTS 
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DATASET  A  VERSUS  DATASET  B,  197C  TO  1975 
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3.  If  A  is  a  dataset  as  in  1,  then 
SHOW  BAR  A 

produces  the  following  plot 


DATASET  A,  1970  TO  1975 
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SHOW  TEXT' 


3.30  The  SHOW  Statement:  -  with  TEXT 


When  used  in  conjunction  with  TEXT,  the  SHOW  statement  vil 
generate  a  textual  display  chat  may  be  used  for  descriptive  purposes 
documentation  etc.  A  combination  of  scored  information  and  textual 
strings  may  be  displayed. 

Svntax 


SHOW  TEXT  (list) 


or 


SHOW  TEXT  =  CENTER  (list) 


SHOW  TEXT  =  LEFT  (list) 


left  justifies  "list" 
centers  "list" 
left  justifies  "list" 


or 


TELL 
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3.31  The  TELL  Statement  -  with  String 


The  TELL  statement  "telLs"  the  user  something,  i.e.,  it  issues 
a  message,  during  execution. 

Syntax 

i - : - 

TELL(string) 


string  is  a  message  of  up  to  125  characters. 


TELL  DATA 


T_ 


3.32  The  TELL  Statement  -  with  DATA 


This  statement  allows  the  user  to  display  a  complete  or  partia 
listing  of  the  datasets  in  a  program.  Each  item  in  the  listing  contains 
the  identifier  and  the  label  of  the  dataset  listed. 

Svntax 


TELL  DATACdata.,  data data  ) 

_ 1 _ 1 _ n 

data.  is  the  identifier  of  the  ith  dataset  to  be  listed.  If  no 

l 

dataset  identifiers  are  specified,  then  all  datasets  in 
the  program  segment  currently  in  working  space  will  be 
listed. 


ASK 
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3.33  The  ASK  Statement:  -  General 


This  statement  "asks"  the  user  a  question  and,  depending  upon 
his  response,  executes  one  of  a  number  of  alternative  procedures. 

Svntax 


ASK(question) .response^ 

procedure^ 

ELSE  response^ 

I  procedure, 

I  1 

1  END  ASK _ 


question  is  a  question  or  statement  eliciting  a  user  response. 

response1  is  a  single,  non  blank  work  (no  special  characters), 

denoting  a  possible  response  to  the  question, 
procedure^  is  the  set  of  instruction(s)  to  be  executed  if  the 
user  enters  "response^". 

ELSE  response^  begins  the  ith  alternative  response/procedure 
branch  of  the  ASK  statement. 

procedure.^  is  the  set  of  instructions  to  be  executed  if  the 
user  enters  "response^". 

Comments 


The  last  response  entry  in  the  ASK  may  be  blank  in  which  case 
the  procedure  following  that  ELSE  statement  will  be  executed  if  any 
response,  other  than  the  previously  specified  ones,  is  given. 


ASK 


3.33.1  The  ASK  Statement  -  with  ELSE  DAlA=indirect 
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This  statement  "asks"  a  question  and  allows  the  user  to  select  a 
single  dataset  from  a  list  of  previously-defined  datasets  to  be  used  in  the 
following  procedure. 

Svncax 


ASK (question) .response^ 

ELSE  DATA=indirect (data  list^) 
procedure^ 


ELSE  DATA*indirect (data  list  ) 

n 

procedure^ 

EOT  ASK 


response^ 
ELSE  DATA= 
indirect 
data  list^ 

procedure^ 

data  list^ 
procedure^ 


passes  control  to  END 

is  used  to  specify  a  list  of  datasets 

is  a  previously-defined  indirect  data  set 

is  any  list  of  previously-defined  data  sets,  excluding 

indirects 

is  the  set  of  instructions  to  be  performed  on  indirect 
(i.e.,  on  the  data  set  specified  by  the  user  from  data 
list^) 

is  any  list  of  previously-defined  data  sets,  excluding 
indirects 

is  the  set  of  instructions  to  be  performed  on  indirect 
in  data  list^. 


Comments 


If  data  list^  is  left  blank  any  defined  dataset  may  be  substituted 


for  the  indirect. 


I- 


ASK 


3.33.2  The  ASK  Statement  -  with  ELSE  INDEX* index 


This  statement  "asks"  a  question  which  allows  the  user  to 
select  a  setting  for  an  index.  This  is  a  conversational  SET  statement. 

Svntax 


ASK(question) , response 


I  ELSE  INDEX=index 

i 

j  procedure 
I  END  ASK 


question 


response 


procedure 


ELSE  INDEX 


index 


is  the  message  asking  the  user  to  select  a  setting 
for  an  index 

is  the  response  to  be  entered  if  the  user  does  not 
wish  to  modify  the  present  index  setting, 
is  the  set  of  instructions  to  be  executed  after  the 
index  setting  operation.  These  instructions  are 
affected  by  this  setting. 

puts  an  unconditional  SET  on  an  index,  using  user 
designated  index  entries.  It  nullifies  and  clears 
any  previous  SET  commands  on  that  index 
is  the  identifier  of  a  previously-defined  index  whose 
setting  is  being  modified. 


4.0  ERROR  MESSAGES 


If  in  formulating  a  statement  the  user  does  not  conform  to  the 
syntax  rules  of  the  User  Language,  an  error  diagnostic  will  be  issued 
immediately  following  the  statement  containing  the  error.  In  interactive 
or  on-line  mode,  the  user  may  correct  the  statement  in  error  and  continue. 
In  batch  mode,  however,  program  execution  will  stop  immediately  after  the 
error  diagnostic,  unless  the  system  parameter  DEBUG  is  on.  If  the  svstem 
is  in  QE3UG  mode,  it  will  continue  processing  the  program  even  if  errors 
are  encountered  and  error  messages  issued.  With  this  option  on,  the  user 
obtains  a  listing  of  all  error  messages  in  a  single  pass  of  the  program 
through  the  system  intepreter. 

The  system  issues  the  following  diagnostics: 

(1)  ILLEGAL  TYPE  n  TIE  IN  COLUMN  m  VALUE  =  v 

n  is  an  integer  code  taking  the  values  1,  2,  3,  or  4 

depending  on  the  type  of  the  TIE  (Textual  Input  Element) 
in  error.  As  discussed  elsewhere  (Part  2.2),  there  are 
four  types  of  TIEs : 

Type 

special  character 
alphanumeric 
integer 
real  number 

m  is  the  column  number  of  the  beginning  character  of 

the  illegal  TIE 

is  the  value  of  the  illegal  TIE 


TIE 

1 

2 

3 

4 


v 
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(21  ILLEGAL  END  OF  INSTRUCTION 

This  diagnostic  is  issued  whenever  the  preceding  state¬ 
ment  is  not  completed  properly.  This  usually  happens 
when  left  and  right  parentheses  are  not  balanced  properly. 
Also,  it  happens  when  the  continuation  of  a  long  state¬ 
ment  is  not  specified  properly  (with  a  comma  or  an  oper¬ 
ator)  at  the  end  of  the  preceding  input  line  or  card. 

(3)  MISSING  STRING  TERMINATOR 

This  diagnostic  is  issued  whenever  a  descriptor  exceeds 
125  characters  in  length. 

(4)  DOUBLY  DEFINED  IDENTIFIER 

This  message  is  issued  after  an  attempt  to  define  an 
identifier  that  has  previously  been  assigned  to  a 
structural  element  (index,  dataset,  model,  graph, 
table,  or  program  module).  Identifiers  must  be  unique 
within  a  structural  element  type.  However,  it  is  good 
practice  to  have  unique  identifiers  across  all  struc¬ 
tural  elements  of  a  program. 

(5)  UNRECOGNIZABLE  COMMAND 

This  message  is  issued  when  the  first  TIE  of  the  pre¬ 
ceding  statement  is  a  Type  2  TIE  but  does  not  belong 
to  the  set  of  primary  statement  words  listed  in  Table 
1-1  or  the  set  of  previously  defined  dataset,  model  and 
table  identifiers. 
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(6)  INDECIPHERABLE  ENTRY 

This  message  is  issued  when  the  first  TIE  of  che  ?re 
ceding  statement  is  not  a  Type  2  TIE.  Note  that  all 
statements  in  a  program  must  begin  with  a  legal  Type 
2  TIE. 

(7)  ILLEGAL  E'.  TRY  DURING  DATA  READ 

This  message  is  issued  whenever  a  dataset  or  graph 
value  entered  during  a  READ  operation  is  in  error. 

A  dataset  or  graph  entry  is  illegal  if  it  is  not 
numeric  (integer  or  real). 

(8)  EMPTY  DATA  BASE 

This  message  is  issued  when  the  file  parameter  of  an 
OPEN  statement  is  misspecif ied ,  i.e.,  is  an  unaccept 
able  physical  unit  number. 

(9)  LOAD  MODULE  NOT  DEFINED  ON  DATA  BASE 

This  message  follows  a  LOAD  statement  attempting  to 
load  from  the  database  a  nonexistent  program  module, 
i.e.,  a  module  not  previously  saved  on  the  database. 

(10)  TABLE  NOT  DEFINED  ON  DATA  BASE 

This  message  is  issued  whenever  reference  is  made 
to  a  TMS  array  number  that  has  not  been  properly  de¬ 
fined  on  a  TMS  database  (using  the  ADD  statement) . 

( 
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(11)  DATA  BASE  SIZE  EXCEEDED 

This  message  is  issued  after  an  ADD  (or  a  SAVE)  state¬ 
ment  attempting  to  save  on  a  database  a  T>IS  array  (or 
a  program  module)  whose  size  exceeds  the  allocated  stor¬ 
age  capacity  of  the  database.  This  usually  happens  be¬ 
cause  the  database  has  not  been  properly  opened  using 
the  OPEN  statement. 

(12)  INSUFFICIENT  SPACE  FOR  DATA  VALUES 

This  message  is  issued  whenever  a  storage  allocation 
statement  (a  READ  or  an  equation)  requires  more  space 
than  the  allocated  storage  capacity  of  a  program. 

The  default  storage  capacity  of  a  program  is  10000 
words . 

(13)  INSUFFICIENT  SPACE  TO  PROCESS  STATEMENT 

This  message  is  issued  whenever  a  storage  definition 
statement,  such  as  a  DATA  statement  defining  a  fixed 
dataset,  requires  more  space  than  the  allocated  stor¬ 
age  capacity  of  the  program.  The  storage  capacity  of 
a  program  is  equal  to  10000  words  (default  value). 

Examples  of  error  messages  are  given  in  Figure  1-1.  The  causes/ 
explanations  of  these  sample  error  messages  are  given  below: 
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READY/STAART(ERROR  MESS ACES  1 

1  Ut;  r.SCOOl* IZABLE  COMMAND 

R EA D'i /START ( ERROR  MESSAGES ) 

READY /INDEX  ROW(3 

2  ILLEGAL  END  OF  INSTRUCTION 
PEADY/INDEX  ROW(3)  ,LA2LE(RG'..-S) 

3  ILLEGAL  TYPE  2  TIE  IN  COLUMN  14  VALUE  =  LAELE 
READY/INDEX  RGV/(  3 ) , LABEL (  RON'S  1 
.READY/INDEX  COL/21 

A  ILLEGAL  TYPE  1  TIE  IN  COLUMN  10  VALUE  =  / 

READY/INDEX  C0L(2) ,44 

5  ILLEGAL  TYPE  3  TIE  IN  COLUMN  14  VALUE  =  44 
READY/INDEX  COL(2.5) 

6  ILLEGAL  TYPE  4  TIE  IN  COLUMN  11  VALUE  =  2.5 
READY/INDEX  CCL(2) 

r.EADY/lNDEX  LONG ( 4 )  ,LAE£L(THIS  IS  A'  CHARACTER  STRING  WHICH  IS  M 
/1 25  CHARACTERS  LONG  AND  WHICH  THEREFORE  CANNOT  EL  ACCEPT! 

7  MISSING  STRING  TERMINATOR 
READY/INDEX  fiOW(4) 

8  DCUBLY  DEFINED  IDENTIFIER 
READY/+DATA  A 

5  INDECIPHERABLE  ENTRY 

READY/DATA  A(ROW),(A  DATASET) 

READY/READ  A 
/I  2  X 

10  ILLEGAL  ENTRY  DURING  DATA  READ 

/I  2  4 

READY/DATA  A(COL) 

11  DOUBLY  DEFINED  IDENTIFIER 
READY/OPEN( 12,23) 

12  EMPTY  DATA  BASE 
READY/OPEf'(  1 ,23) 

READY/SYSTEM.  L0AD=1 
READY/LOAD  PROG 

13  LOAD  MODULE  NOT  DEFINED  ON  DATA  BASE 
READY/INPUT  A,TMS( 1 ,2) 

TABLE  NOT  DEFINED  ON  DATA  EASE 
READY/ADD( 1 ,2) , SIZE (100,10) , 


FIGURE  1-1, 


sample  error  messages 


o  o 


STAART  is  not  a  recognizable  command  word:  START  is. 

The  right  parenthesis  required  in  the  index  size 
definition  is  missing. 

TABLE  is  not  an  acceptable  option  for  the  index  state¬ 
ment;  LABEL  is. 

To  define  the  size  of  index  COL  a  left  parenthesis 
must  be  used,  not  a  /. 

The  integer  44  is  a  spurious  TIE  of  Type  3. 

The  size  of  an  index  must  be  specified  as  a  nonnegative 
integer.  The  real  number  2.5  is  inappropriate  as  an 
index'size . 

This  index  label  is  too  long. 

The  index  ROW  is  already  defined. 

Every  program  statement  must  start  with  a  legal  Type  2 
TIE,  except  in  continuation  statements.  The  +  in  Column 
1  is  unacceptable. 

The  X  is  not  an  acceptable  value  for  the  dataset  A.  The 
number  4  is . 

The  identifier  A  has  already  been  assigned  to  another 
dataset . 

The  physical  unit  number  12  is  too  large  for  the  oper- 
ting  system  processing  this  program.  An  integer  less 
than  10  is  always  acceptable. 

The  program  module  PROG  has  not  been  previously  saved 
(using  the  SAVE  statement). 

Table  or  TMS  array  2  has  not  been  defined  previously 
on  the  database  (using  the  ADD  statement) . 
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5.0  THE  FORTRAN  INTERFACE 


A  major  advantage  of  the  present  approach  is  chat  most  of  the 
information  manipulated  by  programs  written  in  the  user  language  can  be 
stored  on  or  retrieved  from  off-line  files  automatically.  The  values 
of  any  dataset  defined  with  a  TMS  parameter  are  saved  and  can  be  refer¬ 
enced  later  by  other  FORTRAN  programs.  As  a  resulc  of  this  fact,  if 
c’ne  user  wishes  to  interface  other  FORTRAN  programs  with  the  system, 
he  can  do  so  simply  by  manipulating  the  same  TMS  arrays  as  are  manipu¬ 
lated  by  the  system. 

The  only  software  link  to  a  TMS  database  is  a  generalized 
FORTRAN  subroutine  which  has  the  same  capabilities  as  the  TMS  op  cion  on 
datasets.  This  subroutine  may  be  called  directly  from  FORTRAN  programs. 
This  subroutine,  SELNDT,  is  described  in  this  chapter.  This  discussion 
assumes  familiarity  with  the  TMS  concepts  anc  with  FORTRAN. 


5.1  When  to  Use  SELNDT 


Subroutine  SELNDT  can  be  used  whenever  one  needs  to  read  values 
from  or  write  values  to  TMS  arrays  in  existing  TMS  databases  via  a  FORTRAN 
program  external  to  the  system.  SELNDT  assumes  that  the  file  already  exists 
and  that  all  TMS  arrays  to  be  referenced  have  already  been  added  to  the 
database.  These  operations  can  most  easily  be  performed  by  the  host  sys cam 
itself;  therefore,  SELNDT  has  not  been  given  these  capabilities. 

For  those  computer  systems  which  require  that  an  "open  database" 
operation  be  performed  before  any  direct  access  can  be  executed  by  a  FORTRAN 
program,  the  "open  database"  statement  must  be  placed  into  the  FORTRAN  pro¬ 
gram  by  the  user.  Generally,  this  statement  is  a  DEFINE.  Users  who  are 
uncertain  about  these  conventions  on  their  particular  systems  should  con¬ 
sult  with  their  computer  staff. 
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5-2  Parameters  for  SELNDT 

Subroutine  SELNDT  has  the  following  ten  (  13  i  calling  parameters 
or  arguments: 

I0P  Specifies  if  read  (1)  or  write  (0)  is  desired 

FILE  -  Contains  the  physical  unit  number  of  the  file 

TABLE  -  Contains  the  number  of  the  IMS  array 

THR  -  Is  a  scratch  storage  array  used  by  SELNDT 

DAT  -  Contains  or  returns  the  data  values 

BASE  -  Defines  the  base  point 

DIM  -  Defines  the  dimensions  of  DAT 

ENT  -  Defines  the  number  of  entries  in  DAT 

TRANS  -  Defines  the  transposition  to  be  performed 

NDIM  -  Contains  number  of  data  dimensions. 

This  part  discusses  these  parameters  in  detail. 

5.2.1  The  Parameter  TOP 

The  parameter  IOP  is  an  integer  variable.  It  specifies  whether 
values  are  to  be  read  from  the  TMS  array  into  a  core  resident  array  or 
whether  core  resident  values  are  to  be  written  to  the  TMS  array.  A  value 
of  zero  (0)  indicates  write,  while  a  value  of  one  (1)  indicates  read. 

5.2.2  The  Parameters  FILE  and  TABLE 

The  parameters  FILE  and  TABLE  are  both  integer  variables.  They 
specify  the  physical  unit  number  of  the  file  containing  the  TMS  array  and 
the  table  number  of  the  array,  respectively.  Any  file  referenced  by  the 
FILE  parameter  has  to  have  been  previously  created  by  the  host  system.  A 
table  on  that  file  referenced  by  the  TABLE  parameter  has  to  have  been  pre¬ 
viously  added  to  the  file  by  che  host  system. 


5.2.3  The  Parameter  THR 


The  parameter  THR  is  a  one-diner.sional  integer  array  needed  by 
SELNDT  for  scratch  storage.  The  user  must  allocate  at  least  1  -  NO  -r  NV 
words  to  this  array  in  his  calling  program.  ND  represents  the  maximum 
number  of  array  dimensions  associated  with  ar.y  T'.'S  array  to  be  manipulated. 
NV  represents  the  maximum  number  of  virtual  sheets  associated  with  ar.v 
TMS  array  to  be  manipulated.  Usually,  ND=5  ar.b  N  .'=0 ;  thus,  the  dimension 
of  THR  is  7. 

From  the  standpoint  of  the  calling  routine  only  the  first  cell 
of  the  THR  is  important.  If  this  cell  returns  a  nonzero  value,  then  the 
access  was  successful.  If  the  ceil  ret-rns  oerc,  then  the  access  was  not 
successful . 

5.2.4  The  Parameter  DAT 


The  parameter  DAT  is  a  real  array  which  contains  the  core  resi¬ 
dent  values  being  read  or  written.  Its  dimensionality  is  defined  by  the 
DIM  and  NDIM  parameters  in  Che  calling  sequence. 

5.2.5  The  Parameters  DIM  and  NDIM 


The  parameters  DIM  and  NDIM  are  used  to  define  the  dimensionality 
of  the  parameter  DAT.  NDIM  is  an  integer  variable  which  contains  the  number 
of  dimensions  in  DAT.  DIM  is  a  one-dimensional  integer  array  which  contains 
NDIM  entries.  The  values  in  DIM  should  be  identical  both  in  magnitude  and 
order  to  the  values  in  the  DIMENSION  declaration  for  DAT.  Thus,  if  DAT  is 
dimensioned 

DATUO,  20,15) 

.  " : '  would  be  initialized  as  follows 
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5.2.6  The  Piraaecer  ENT 

- - - -  ^ 

The  parameter  ENT  is  a  one-dimensional  integer  array  containing 
NDIM  entries.  It  is  used  to  define  the  actual  number  of  entries  to  be 

processed  for  each  dimension  in  the  DAT  array.  ' 

5.2.7  The  Parameter  BASE 


The  parameter  BASE  is  a  one-dimensional  integer  array  containing 
an  entry’  for  each  dimension  in  the  TMS  array.  This  array  simply  defines 
the  base  point  of  the  core  storage  array  relative  to  the  TMS  array. 

5.2.S  The  Parameter  TRANS 

Tne  parameter  TRANS  is  a  one-dimensional  integer  array  containing 
an  entry  for  each  dimension  in  the  TMS  array.  This  array  defines  the  trans¬ 
position  relationship  between  the  core  storage  array  and  the  TMS  array.  If 
TRANS  (I)  is  zero,  then  the  Ith  dimension  on  the  TMS  array  is  invariant  re¬ 
lative  to  the  core  storage  array  —  i.e.,  it  has  the  value  specified  for  the 
base  point.  If  TRANS  (I)  equals  J,  then  the  Jth  core  storage  array  dimen¬ 
sion  moves  along  the  Ith  TMS  array  dimension. 

Note  that  if  word  one  of  TRANS  simply  contains  a  minus  one,  then 
no  transposition  is  performed. 
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APPENDIX  D.  USER'S  GUIDE 


The  General  Aviation  Dynamics  Computer  System 


Users  Guide 


The  General  Aviation  Dynamics  model  (GAD)  is  implemented  in 
NUCLEUS*,  a  computer  software  system  developed  at  Battelle.  It  resides 
currently  on  the  Battelle  computer  system  and  on  the  United  Computing 
Systems  (UCS)  computer,  and  it  can  be  accessed  remotely  on  either  machine 
via  telephone  from  anywhere  in  the  U.S.  Access  is  possible  via  either  a 
remote  batch  terminal  or  an  interactive  terminal,  providing  that  an 
authorized  user  name  and  password  for  the  particular, computer  system 
being  accessed  are  available.  Procedures  for  accessing  the  GAD  model 
on  the  Battelle  and  UCS  computers  are  given  in  Sections  I  and  II; 

Section  III  shows  an  example  of  an  interactive  session  with  GAD  and 
Section  IV  describes  batch  use  of  the  model. 

Section  I;  Accessing  on  the  Battelle  computer 

The  GAD  model  resides  on  the  Battelle  computer  and  may  be  accessed 
via  telephone  in  either  remote  batch  mode  using  a  remote  batch  terminal  that 
is  appropriately  linked  to  the  Battelle  computer,  or  in  on-line  or  inter¬ 
active  mode  using  an  interactive  terminal.  Details  of  the  actual  access 
procedure  for  interactive  use  follow. 


*  NUCLEUS:  Numerical  Classification  and  Evaluation  System 
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1.1  Login/Logout  Procedures  for  the  Battelle  Computer  System 

via  TYMNET 


(1)  Dial  TYMNET  number  from  the  list  in  Table  1. 

Wait  for  data  tone  from  TYMNET  connection. 

(2)  Place  telephone  headset  in  acoustic  coupler,  or 
initialize  other  type  of  modem. 

(3)  The  message 

PLEASE  ENTER  TERMINAL  IDENTIFIER 
is  printed  or  displayed  on  the  terminal.  On 
some  terminals  this  message  will  appear  garbled. 

If  this  occurs  simply  continue  with  Step  (4). 

(4)  Type  the  appropriate  identifier  for  your  terminal 
from  the  list  in  Table  2. 

(5)  TYMNET  responds: 

PLEASE  LOGIN: 

(6)  Login  by  typing  the  sequence 

(CTRL  H) BCL6400 ; INTERCOM(CR) 

where 

(CTRL  H)  denotes  simultaneous  keying  of 
the  control  and  H  keys,  and  (CR)  denotes 
keying  of  the  carriage  return  (or  cursor  return) 
key. 

No  blanks  (spaces)  are  allowed  in  the  above 
sequence.  Wait  a  few  seconds  (about  10)  for  con¬ 
nection  with  the  Battelle  INTERCOM  system.  The 
TYMNET  username/password  combination  (BCL6400; INTERCOM) 
should  be  treated  with  discretion  and  made  avail¬ 
able  only  to  users  with  valid  Battelle  username/ 
password  combinations. 

(7)  INTERCOM  responds: 

BATTELLE  INTERCOM  4.5 

DATE  MM/DD/YY 

TIME  HH .MM. SS 

PLEASE  LOGIN 

(8)  Enter  the  following  sequence: 

LOGIN , us  e rname , p  as  sword , SUP ( CR) 
where  username  and  password  are  valid  code  words 
Issued  by  the  Battelle  Computer  Center.  No  blanks 
(spaces)  are  allowed  in  the  above  sequence. 

(9)  INTERCOM  responds: 

COMMAND- 

(10)  Carry  out  dialogue  with  INTERCOM  (see  Login/Logout 
Procedures  via  Direct  Dial-Up,  Step  6). 

(11)  To  exit  from  INTERCOM  normally,  enter 

LOGOUT 

after  the  INTERCOM  prompt  COMMAND- 
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(12)  To  exit  from  INTERCOM  abnormally  (i.e.,  abort) 
enter  the  sequence 
Z  A 

after  a  system  pause.  Or,  if  in  the  middle  of  a 
lengthy  response,  enter  the  sequence 
(CRTL  Z)  %  A 

. where 

(CRTL  Z)  denotes  the  simultaneous  keying 
of  the  control  and  Z  keys.  This  keying 
causes  the  system  response  to  halt. 

After  an  abnormal  termination,  the  system  responds 
USER  ABORT 
COMMAND- 

At  this  point,  you  may  log  out  normally. 


(13) 


After  LOGOUT,  the  system  gives  you  an  itemized 
summary  of  your  interaction  with  the  Battelle 
Computer.  This  summary  looks  like  the  following 


CP  X.XXX  SEC. 

10  Y.YYY  SEC. 

CM  Z.ZZZ  SEC 

SS  S.SSS  SEC 

CH  CCC  CHARS . 

CONNECT  TIME  H  HRS.  M  MIN.  MM/DD/YY 
LOGGED  OUT  AT  HH.MM.SS. 


where 


CP  is  the  central  processor  time  of  your  run 
10  Is  a  measure  of  your  use  of  peripheral  equipment 
CM  is  a  measure  of  your  use  of  central  memory 
SS  is  the  combined  amount  of  computer  time  for  which 
you  are  charged 

CH  is  the  number  of  characters  transmitted 
CONNECT  TIME  is  the  time  between  login  and  logout. 


The  cost  for  TYMNET  services  is  approximately  17c  per  connect  minute  which 
is  in  addition  to  the  normal  INTERCOM  processing  cost  incurred.  These 
charges  are  subject  to  change. 


NOTES:  •  To  backspace,  type  CTRL  and  H  simultaneously 

•  To  skip  a  line  of  entry,  type  CTRL  and  X  simultaneously 

•  To  suspend  printing,  type  CTRL  and  Z  simultaneously 

•  To  abort,  enter  %,  then  A 

•  All  user  entries  must  be  terminated  by  (CR) . 
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The  following  TYMNET  messages  may  occur  if  a  TYMNET  connection 
is  not  immediately  available. 


Message 


Meaning 


ALL  PORTS  BUST  OR  HOST  OCT  OF 
PORT 


A  connection  with  the  boat  computer  cannot  be  made  at  thla  time  because  the 
TYMNET  porta  on  that  host  are  all  busy.  Try  again  in  a  few  minutes. 


HOST  NOT  AVAILABLE  THROUGH  NET  Thia  message  can  occur  if:  1)  the  TTMCOM  or  its  nelghbor(s)  is  down,  2)  an 

invalid  host  has  been  requested  (perhaps  mistyped),  or  3)  a  new  Supervisor 
la  taking  over  the  netwok  and  has  not  yet  picked  up  that  host.  Try  again 
in  a  few  minutes. 


HOST  DOWN  OR  HOST  "n*  DOWN  The  network  is  fully  operational  but  the  host  computer  Itself  la  down. 

HOST  SHUT  The  Supervisor  has  been  notified  to  route  no  new  uaers  to  the  host: 

current  users  are  not  affected. 


DROPPED  BY  HOST  SYSTEM  The  user  has  logged  off  and/or  has  been  disconnected  by  the  host  computer. 

If  the  session  is  finished,  hang  up;  if  not,  try  again. 

CIRCUITS  BUSY  All  available  paths  to  the  TYMCOM  are  busy;  try  again  in  a  few  minutes. 

A  new  Supervisor  is  taking  over  the  network  and  a  log-in  will  normally  be 
possible  in  a  few  minutes. 


TRY  AGAIN  IN  5  MINUTES 


TABLE  1:  TYMNET  TELEPHONE  NUMBERS 
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Table  I  (Continued) 


THE  FOLLOWING  LOCATIONS  NOW  HAVE  1200  BAUD  SERVICE  AVAILABLE 

UPDATE  AS  OF  01/05/78 


NAME 

LOCATION 

PHONE  NO. 

MODEM  TYPE* 

SF2 

SAN  FRANCISCO,  CA 

415/421-7121 

1200 

VADIC 

SF2 

SAN  FRANCISCO,  CA 

415/982-3770 

1200 

BELL  212 

SF2 

SAN  FRANCISCO,  CA 

415/397-8461 

1200 

BELL  202S 

HOU3 

HOUSTON,  TX 

713/780-7596 

1200 

VADIC 

HOU3 

HOUSTON,  TX 

713/780-7593 

1200 

BELL  202S 

DALI 

DALLAS ,  TX 

214/688-1117 

1200 

BELL  212 

LA  2 

LOS  ANGELES,  CA 

213/687-8083 

1200 

VADIC 

LA2 

LOS  ANGELES,  CA 

213/629-3001 

1200 

BELL  202S 

LA  2 

LOS  ANGELES,  CA 

213/626-0365 

1200 

BELL  212 

CRP3 

SAN  JOSE,  CA 

408/446-6932 

1200 

BELL  202S 

CRP3 

SAN  JOSE,  CA 

408/446-7001 

1200 

VADIC 

DET3 

DETROIT,  MI 

313/963-2353 

1200 

VADIC 

CRP4 

SAN  JOSE,  CA 

408/446-7309 

1200 

BELL  212 

DET3 

DETROIT,  MI 

313/963-4676 

1200 

BELL  202S 

DET3 

DETROIT,  MI 

313/963-2353 

1200 

VADIC 

WAS  3 

WASHINGTON,  DC 

703/841-9330 

1200 

VADIC 

WAS  3 

WASHINGTON,  DC 

703/525-6290 

1200 

BELL  212 

WAS  3 

WASHINGTON,  DC 

703/527-7106 

1200 

BELL  202S 

SL1 

ST  LOUIS,  MO 

314/621-4660 

1200 

BELL  212 

NY  3 

NEW  YORK,  NY 

212/689-8910 

1200 

VADIC 

NY  3 

NEW  YORK,  NY 

212/689-8850 

1200 

BELL  202S 

BOS  3 

BOSTON,  MA 

617/965-5520 

1200 

VADIC 

BOS  3 

BOSTON,  MA 

617/244-1240 

1200 

BELL  202S 

CHI  3 

CHICAGO,  IL 

312/372-0391 

1200 

VADIC 

CHI3 

CHICAGO,  IL 

312/368-0022 

1200 

BELL  202S 

ATLSR 

ATLANTA,  GA 

404/581-0619 

1200 

BELL  212 

SCOSR1 

HOUSTON,  TX 

713/977-7671 

1200 

BELL  212 

*-WHEN  USING  1200  BAUD  SERVICE  THE  USER  MUST  SELECT  A  PHONE 
NUMBER  EQUIPPED  WITH  THE  SAME  TYPE  OF  MODEM  AS  THE  MODEM 
CONNECTED  TO  THE  USER’S  TERMINAL. 


TABLE  2.  USER  TERMINAL  IDENTIFICATION  CODES 
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1.2  Login/Logout  Procedures  for  the  Battelle  Computer  System 

and  Che  NUCLEUS  System 
via 

Direct  Dial-Up 


( 


(1)  Dial  (614)424-5850.  Telephone  will  ring;  a  high- 
pitched  tone  follows. 

(2)  Place  receiver  in  coupler.  The  green  carrier  light 
will  come  on.  Terminal  is  connected  to  computer. 

(3)  Key  CR  (carriage  return)  to  indicate  the  terminal's 
baud  (300  or  110)  to  the  Battelle  CDC  equipment. 

(4)  A  message  like  the  following  is  displayed  or  printed 
on  the  terminal: 

BATTELLE  INTERCOM  4.5 
DATE  MM/DD/YY 
TIME  HH .MM. SS . 

PLEASE  LOGIN 

(5)  Enter  the  following  sequence: 

LOGIN, username, password, SUP. (CR) 
where  username  and  password  are  valid  code-words  issued 
by  the  Battelle  Computer  Center.  No  blanks  (spaces) 
are  allowed  in  the  above  sequence. 

NOTE:  To  correct  mistakes  In  a  given  line  of 
entry,  backspace  to  the  character  in 
error  by  pressing  the  CTRL  and  H  keys 
simultaneously  as  many  times  as  necessary. 

Then  retype  the  remainder  of  the  line. 

NOTE:  To  delete  line  of  entry  altogether,  press 
the  CTRL  and  X  keys  simultaneously.  Then 
continue  entering  correct  information 

NOTE:  (CR)  denotes  carriage  return  (or  cursor  return). 

(6)  The  Battelle  INTERCOM  system  responds: 

COMMAND- 

To  extend  the  time  allowed  for  the  on-line  interaction, 

enter 


ETL ,  100 


This  ensures  enough  time  for  the  entire  interactive 
dialogue,  calculations,  and  desired  output.  The 
system  will  again  respond: 

COMMAND 

(7)  To  attach  the  GAD  procedure  file  enter 

ATTACH , PROFIL .AVIATION , ID=NUCLEUS ,MR=1 . 

(8)  At  the  system  response,  COMMAND,  begin  execution 
of  the  GAD  procedure  by  entering 

BEGIN, GAD 

The  NUCLEUS  system  will  respond 
READY/ 

To  execute  the  dynamic  model  enter 
LOAD  GAD 

(10)  To  terminate  abnormally  at  any  time  the  interaction 
with  NUCLEUS: 

•  If  at  a  system  pause,  enter  X,  then  A,  then  (CR) . 

•  If  in  the  middle  of  a  lengthy  NUCLEUS,  response, 
enter  CTRL  and  Z  simultaneously,  then  X,  then 

A,  then  (CR) . 

(11)  INTERCOM  responds: 

USER  ABORT 
COMMAND- 

(12)  To  exit  from  INTERCOM,  enter 

LOGOUT 

any  time  you  have  the  prompt  COMMAND 

(13)  After  an  exit  from  NUCLEUS,  you  are  back  In  INTERCOM 
with  the  prompt 

COMMAND 

you  may  enter  the  dynamic  simulation  model  again  or  end 
your  interaction  with  the  system 

(14)  To  exit  from  INTERCOM,  enter 

LOGOUT 

(15)  After  LOGOUT,  don't  forget  to  hang  up  the  phone. 

Examples  of  login/logout  procedures  are  given  in  Table  3. 


TABLE  3.  LOCIN/ LOGOUT  PROCEDURES  -  EXAMPLES 
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Section  II;  Accessing  GAD  on  the  UCS  Computer 

The  GAD  model  also  resides  on  the  UCS  computer  and  may  be  accessed 
via  telephone  in  either  remote  batch  mode  using  a  remote  batch  terminal  that 
is  appropriately  linked  to  UCS,  or  in  on-line  mode  using  an  on-line  or  inter¬ 
active  terminal.  Details  of  the  actual  access  procedure  for  interactive  use 
follow. 


II. 1  Logging  into  the  System 


Your  UCS  representative  must  furnish  to  you  the  local  telephone 
number  which  connects  terminals  in  your  city  with  the  UCS  National  Data 
Center  in  Kansas  City,  Mo.  He  will  also  furnish  you  with  a  unique  seven- 
character  user  number  with  which  you  need  to  identify  yourself  to  the 
computer  each  time  you  log-in.  This  user  nunber  is  assigned  to  you  for 
security  and  for  billing  purposes. 

Log-in  by  following  this  sequence: 

•  If  your  terminal  is  equipped  with  a  HALF/FULL  duplex  switch, 
set  it  to  HALF. 

•  Be  sure  that  the  coupler  is  turned  on  if  it  has  controls 
separate  from  those  of  your  terminal. 

•  Turn  on  the  teletypewriter  if  the  controls  are  separate 
from  the  coupler. 

•  Dial  the  special  phone  number  and  wait  until  you  hear  a 
high-pitched  tone  like  a  whistle.  This  means  that  the 
computer  has  a  free  line  which  it  has  just  dedicated  to 
you. 

•  Place  the  telephone  receiver  in  the  acoustical  coupler 
of  your  terminal. 

•  As  soon  as  the  connection  is  made  between  your  terminal 
and  the  central  site  computer,  you  will  notice  a  signal. 
Depending  on  the  terminal  speed  and  type,  your  terminal 
may  give  a  small  "jump1',  a  button  will  light  up  or  some 
output  will  be  printed.  For  example,  the  10  cps  tele¬ 
type  will  print  L? 

•  Your  response  will  depend  on  the  speed  of  your  terminal 
and  on  the  system  code  you  have  been  Instructed  to  use. 
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If  you  have  any  questions 
concerning  the  appropriate 
response  for  your  particular 
type  of  terminal,  consult  your 
UCS  representative. 


•  The  operating  system  will  print  some  heading  information  con¬ 
sisting  of  date,  time  and  communication  information  and  then 
request  your  user  number.  Type  it  in  and  return  the  carriage. 

•  The  system  will  then  print  an  overscored  area  in  which  you 
must  enter  your  password.  Actually  type  the  password  over 
the  blackened  letters;  this  prevents  anyone  from  seeing  your 
password  if  you  leave  your  terminal.  If  you  are  not  interested 
in  the  security  of  the  overprinting,  you  may  enter  the  pass¬ 
word  after  your  user  number  on  the  same  line.  Simply  type  a 
comma  (or  a  blank)  and  the  password  directly  after  the  user 
number. 

•  If  your  user  number  and  password  are  accepted  by  the  system, 
a  message  prints  and  you  are  ready  to  begin.  If  either  your 
user  number  or  password  is  invalid,  the  message  ILLEGAL  LOG¬ 
IN,  TRY  AGAIN  prints,  and  you  must  enter  both  again.  The 
system  allows  you  five  tries  to  enter  the  proper  combination 
and  then  will  print  ILLEGAL  TERMINAL  and  disconnect  you. 


Line  Speeds 

10 

13.4 

15 

30 

120 

System  Codes 

61 

62 

63 

65 


Log-in  Codes 

? 

3 

8 

T 

P 


L?  Enter  login  code 
UNINET  NODE  12  CHNL  40331 

SERVICE  :  Enter  system  code 

UCS  02/06/79.  10.22.10.  Q113 

USER  NUMBER:  Enter  user  number 
XKAXXXXX  Enter  password 
PROJECT-ID  REQUIRED:  Enter  project- id 
*RDY-F0R* 


iJ-T-A  •  1  •*»>••• 


Notes:  •  Teletypewriters  having  different  keyboards  than  above  will 

usually  perform  the  same  functions  if  the  key  positioned 
in  the  same  place  is  used,  e.g.  for  use  Shift  0. 

•  To  delete  one  character,  type 

•  To  delete  a  line,  type  ALTMODE. 

•  To  abort  while  the  terminal  is  printing,  enter  S. 

•  To  abort  while  the  terminal  is  not  printing,  enter  STOP. 

•  All  user  entries  must  be  terminated  by  RETURN. 
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11. 2  Accessing  GAD  interactively  on  the  UCS  computer 

(1)  After  the  log-in  procedure  the  system  will  respond 

*RDY-FOR* 

Execute  the  GAD  command  file  by  entering 
CMD,  GAD 

(2)  The  system  will  respond 

01/25/79.  15.20.23. 

PROGRAM  NUCLEUS 

READY  ? 

Execute  the  GAD  model  by  entering 
LOAD,  GAD 

11. 3  Log-Off  Command 

When  a  user  is  finished  working  at  the  terminal,  he  should  dis¬ 
connect  -from  the  system  as  soon  as  possible  because  connect  time  is  a 
billable  resource.  This  is  done  by  entering  GOODBYE  (or  BYE). 

The  GOODBYE  Command 

The  GOODBYE  command,  or  the  abbreviation  BYE,  disconnects  the 
terminal  from  the  time-sharing  system  and  prints  the  total  connect  time 
for  the  session  plus  applicable  service  unit  statistics  for  your 
billing  option. 

Command  format :  GOODBYE 
or 
BYE 
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Section  III.  AN  EXAMPLE  OF  THE  GENERAL  AVIATION  DYNAMICS 
MODEL  IN  THE  INTERACTIVE  MODE 


Sensitivity  Analysis 

In  general,  there  are  two  ways  to  use  model  results  or  simulations  - 
individually  as  projections  and  in  pairs  as  sensitivity  measures.  Use  of 
the  model  simply  to  make  projections  is  fraught  with  dangers.  Many  potential 
users  will  not  understand  how  the  projections  were  derived  and  will  expect 
unreasonable  accuracy.  The  model  is  better  used  by  employing  extensive 
sensitivity  analysis  to  evaluate  a  range  of  policies  under  a  range  of 
exogenous  conditions.  This  process  will  identify  the  principal  areas  of 
model  uncertainty  and  those  portions  of  the  model  that  deserve  the 
greatest  additional  research. 

The  logical  structure  of  the  GAD  model  has  been  constructed  such 
that  relative  comparisons  can  be  made  between  the  model  forecasts  from  any 
two  simulations.  In  particular,  during  a  sensitivity  analysis,  absolute 
forecasts  for  each  simulation  are  available,  as  well  as  percent  deviations 
between  the  two  cases.  These  deviations  can  be  displayed  over  time  either 
graphically  or  in  tabular  format. 

A  sensitivity  analysis  can  be  performed  between  any  two  simulations 
which  are  compatible  with  the  model’s  capabilities.  All  GAD  model  output 
data  from  the  first  simulation  are  stored  on  a  separate  file.  This  base 
case  need  not  be  the  "baseline"  forecast  representative  of  expected  future 
conditions,  but  can  be  the  result  of  any  consistent  set  of  conditions  chosen 
by  the  analyst.  Intermediate  absolute  forecast  results  from  this  base  case 
can  be  obtained  by  the  enalyst,  if  desired.  After  obtaining  all  required 
intermediate  output,  the  second  simulation  is  specified  and  run.  Absolute 

results  of  the  second  simulation  are  elso  available  to  the  enalyst. 

•  0 

Sensitivity  results  are  derived  within  the  program  logic  by  subtracting 
the  results  of  the  first  simulation  from  the  second  simulation,  dividing 
by  the  first  simulation,  and  multiplying  by  100  to  convert  differences  to 
percent  deviations  from  the  bese  case;  mathematically , 
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aa(i.j)2-  AA(I,J) 

X  Deviation  - - x  100 

AAd.Dj^ 

where, 

AA(I,J)^  »  the  number  of  active  aircraft  of  type  J  within 
category  I  from  the  first  (base)  simulation 
AACl.J)^  ■  the  number  of  active  aircraft  of  type  J  within 
category  I  from  the  second  simulation 

Values  for  these  parameters  are,  of  course,  obtained  at  the  same  instant  in 
time  during  their  respective  simulations. 

Should  conditions  within  the  second  simulation  not  change  imme¬ 
diately  from  the  base  case,  percent  deviations,  until  the  change  becomes 
effective,  will  be  zero.  Furthermore,  by  continually  computing  these 
deviations  over  time,  the  non-linearity  in  model  response  is  preserved. 

Most  previous  sensitivity  analyses  of  general  aviation  activity  were 
predicated  on  either  linear  or  log- linear  sensitivities. 

An  Example 

The  GAD  model  uses  the  interactive  dialogue  feature  of  NUCLEUS 
to  guide  the  analyst  through  a  series  of  procedures  and  options.  This 
technique  eliminates  the  need  for  preliminary  calculations  by  the  user. 

Simple  yes/no  responses  to  NUCLEUS  questions  establish  the  conditions  of 
the  particular  simulation  to  be  run.  If  the  user  is  uncertain  of  the 
parameter  values  contained  in  the  model,  NUCLEUS  will  display  them.  If 
the  user  desires  to  change  these  values,  NUCLEUS  will  accept  the  new  values. 
Incorrect  (or  unexpected)  responses  to  NUCLEUS  questions  vill  simply  cause  the 
same  question  to  be  repeated. 

•  / 

A  sensitivity  example,  comparing  the  normal  "baseline"  forecast 
to  an  increased  fuel  tax  (effective  January  1,  1979),  is  discussed  below. 

This  example  was  run  on  the  UCS  computer.  Not  all  the  options  available 
for  Input/output  are  displayed;  only  enough  to  illustrate  the  procedures. 

In  this  example,  all  user  entries  are  underlined. 
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Following  the  usual  log-in  procedures  and  LOAD  GADt  the  computer  will 

respond 

YOU  ARE  ENTERING  THE  GENERAL  AVIATION  DYNAMICS  MODEL  CREATED  AT  BATTELLE 
COLUMBUS  LABORATORIES,  WRITTEN  IN  THE  MODELING  LANGUAGE  NUCLEUS.  IN 
THIS  SESSION  YOU  WILL  PROJECT  CERTAIN  LEVELS  OF  GENERAL  AVIATION  ACTIVITY 
FOR  THE  YEARS  1977  TO  1987. 

ENTER  ENDING  YEAR  FOR  SIMULATION 
? 

Any  year  between  197?  and  1987  is  an  acceptable  response;  the  simulation  will 
be  from  1977  through  the  year  specified  and  the  results  will  be  reported  for 
that  range  of  years.  A  response  within  the  acceptable  range 

1980 


allows  the  system  to  continue  with 

Step  1  —  WOULD  LIKE  TO  COMPUTE  THE  FORECAST  WITH  THE  INITIAL  ASSUMPTIONS 
UNCHANGED  (YES  OR  NO)  OR  VIEW  THE  STEPS  OF  THIS  MODEL  (TEACH) 

? 

The  response 


TEACH 


causes  the  steps  of  the  model  to  be  printed  and  then  the  question  to  be  repeated. 
"Teach"  also  causes  the  steps  of  the  model  to  be  printed  out  as  those  steps  are 
executed 

STEP  1  —  COMPUTE  THE  FORECAST  USING  THE  INITIAL  ASSUMPTIONS  UNCHANGED. 

STEP  2  —  DISPLAY  AND/OR  CHANGE  INITIAL  ASSUMPTIONS. 

STEP  3  —  COMPUTE  THE  FORECAST  OF  GENERAL  AVIATION  ACTIVITY. 

STEP  4  —  PRINT  TABLES  OF  RESULTS  OF  THE  FORECAST. 

STEP  5  —  PLOT  THE  RESULTS  OF  THE  FORECAST. 

STEP  6  —  COMPARE  THE  RESULTS  OF  THE  PRESENT  FORECAST  TO  THOSE  OF  A  PREVIOUS 
FORECAST  FOR  SENSITIVITY  ANALYSIS. 

STEP  7  —  PRINT  TABLES  FOR  SENSITIVITY  ANALYSIS. 

STEP  8  --  PLOT  THE  RESULTS  OF  SENSITIVITY  ANALYSIS. 

STEP  9  —  SAVE  THE  RESULTS  OF  THIS  FORECAST  FOR  FUTURE  SENSITIVITY  ANALYSIS. 

STEP  1  —  WOULD  YOU  LIKE  TO  COMPUTE  THE  FORECAST  WITH  THE  INITIAL  ASSUMPTIONS 
UNCHANGED  (YES  OR  NO)  OR  VIEW  THE  STEPS  OF  THIS  MODEL  (TEACH) 

? 
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Now  the  appropriate  response  to  the  question  is  either  YES  or  NO;  a 
response  of 

YES 

causes  the  normal  "baseline"  simulation  to  be  executed.  Since  the  TEACE  flag 
was  set  on  by  the  TEACE  request,  the  steps  are  printed  out  as  they  are  executed 

STEP  3  —  THE  FORECAST  OF  GENERAL  AVIATION  ACTIVITY  IS  BEING  COMPUTED. 

Saving  executed  the  simulation  the  system  prints , 

STEP  4  —  PRINT  TABLES  OF  RESULTS  OF  THE  FORECAST. 

and  then  asks  the  question 

DO  YOU  WANT  TO  SEE  TABLES  OF  RESULTS  OF  THE  FORECAST 
? 

The  response 


YES 


causes  the  system  to  ask 

WHAT  OUTPUT  TABLE  WOULD  YOU  LIKE,  OR  ENTER  LIST  OR  NONE. 

? 

The  unfamiliar  user  will  not  know  the  available  tabular  output  options.  By 
responding 


LIST 


the  following  list  of  output  table  options  will  be  printed. 


IDENTIFIER 

AA1 

AA2 

AIRPORTS 

AIRUTIL 

ECONOMIC 

FIXEDCOST 

VARCOST 

FUEL 

HOURS FLOWN 

OPERATIONS 

PILOTS 

REVENUE 

TOTALS 


DESCRIPTION 

ACTIVE  AIRCRAFT  BY  YEAR 
ACTIVE  AIRCRAFT  BY  USER  CATEGORY 

LOCAL  AND  ITINERANT  OPERATIONS  PLUS  IFR  FLIGHT  PLANS  FILED 
AIRCRAFT  UTILIZATION  RATES 
DPI, (DIP, RAD 
FIXED  COST 

VARIABLE  COST  , 

FUEL  CONSUMED  IN  MILLIONS  OF  GALLONS 
HOURS  FLOWN  IN  THOUSANDS 
TOTAL  OPERATIONS,  IN  THOUSANDS 
SP ,PP ,CP ,ATP ,P ,HP ,TP , IP ,HR, THP 
FEDERAL  TAX  REVENUE 

TOTAL  AIRCRAFT,  TOTAL  HOURS  FLOWN,  TOTAL  OPERATIONS 
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When  the  list  ia  complete,  the  previous  question  will  be  repeated: 

WHAT  OUTPUT  TABLE  WOULD  YOU  LIKE,  OR  ENTER  LIST  OR  NONE. 

? 

A  response  of 

PILOTS 

will  generate  the  following  table 

GENERAL  AVIATION  DYNAMICS  MODEL  PAGE  1 
PILOT  DATA,  1977  TO  1980 


1977 

1978 

1979 

1980 

STUDENT  PILOTS 

188,801 

183,794, 

183,654 

183,176 

PRIVATE  PILOTS 

309,005 

323,821 

335,104 

344,608 

COMMERCIAL  PILOTS 

187,801 

189,699 

192,068 

195,342 

AIR  TRANSPORT  PILOTS 

45,072 

47,784 

50,879 

53,766 

PILOT  SUBTOTAL 

730,679 

745,098 

761,705 

776,891 

HELICOPTER  PILOTS 

4,804 

4,333 

3,940 

3,608 

TOTAL  PILOTS 

735,483 

749,431 

765,645 

780,499 

INSTRUMENT  RATINGS 

211,364 

221,497 

232,437 

243,969 

HELICOPTER  RATINGS 

23,012 

24,395 

25,739 

27,052 

TOTAL  HELIC  RATINGS 

27,816 

28,728 

29,679 

30,659 

Hote  that  the  table  is  printed  for  only  the  requested  years,  1977-1980.  Upon 
completion  of  the  requested  output  table,  the  same  question  is  repeated. 

WHAT  OUTPUT  TABLE  WOULD  YOU  LIKE,  OR  ENTER  LIST  OR  NONE. 

? 

The  user  may  request  as  many  of  the  table  options  as  he  wants.  When  no  more 
tabular  data  is  required  the  response  is 


NONE 
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Since  the  TEACH  flgg  is  on,  the  computer  prints  the  next  step  in  the  model , 
STEP  5  —  PLOT  THE  RESULTS  OF  THE  FORECAST. 
and  then  asks  the  question 

DO  YOU  WANT  TO  SEE  PLOTS  OF  RESULTS  OF  THE  FORECAST 
? 

By  answering 

YES 


the  computer  responds 

WHAT  PLOT  WOULD  YOU  LIKE.  ENTER  THE  VARIABLE,  OR  LIST,  OR  NONE 
? 

Not  being  familiar  with  the  plot  options  the  user  responds 

LIST 

which  will  generate  the  following  list  of  variables: 


IDENTIFIER 

AA 

AASUM 

ATP 

AUR 

CP 

DPI 

FC 

FIX 

FTR 

GNP 

HF 

HFSUM 

HP 

HR 

IP 

OPS 

OPSUM 

P 

PP 

RAD 

SP 

TC 

THP 

TP 

VC 


DESCRIPTION 

NUMBER  OF  ACTIVE  AIRCRAFT 

TOTAL  NUMBER  OF  AIRCRAFT 

AIRLINE  TRANSPORT  PILOTS 

AIRCRAFT  UTILIZATION  RATE  (HRS/AC/YR) 

COMMERCIAL  PILOTS 

DISPOSABLE  PERSONAL  INCOME  (1972  $,  1972-1) 
FUEL  CONSUMED  (MILLION  GALLONS) 

FIXED  COST  INDEX  ($/HR) ,  (1972  $,  1972-1) 
FEDERAL  TAX  REVENUE  (MILLION  DOLLARS) 

GROSS  NATIONAL  PRODUCT  (1972  $,  1972-1) 

HOURS  FLOWN  (THOUSANDS) 

TOTAL  HOURS  FLOWN  (THOUSANDS) 

HELICOPTER  PILOTS 
HELICOPTER  RATINGS 
INSTRUMENT  RATINGS 
OPERATIONS  (THOUSANDS) 

TOTAL  OPERATIONS  (THOUSANDS) 

PILOT  SUBTOTAL 
PRIVATE  PILOTS 

REVENUE  AIRCRAFT  DEPARTURES  (1972,$,  1972-1) 
STUDENT  PILOTS 
TOTAL  COST 

TOTAL  HELICOPTER  RATINGS 
TOTAL  PILOTS 

VARIABLE  COST  INDEX  ($/HR),  (1972  $,  1972-1) 
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followed  by  a  repeat  of  the  question 

WHAT  PLOT  WOULD  YOU  LIKE.  ENTER  THE  VARIABLE,  OR  LIST,  OR  NONE 
? 

Any  variable  identifier  from  the  above  list  can  be  specified.  For  example,  in 
order  to  plot  the  total  number  of  aircraft,  the  user  responds 

AASUM 

Now  the  computer  will  ask 

PLOT  THIS  VARIABLE  AGAINST  TIME  OR  ANOTHER  VARIABLE  OR  LIST 
? 

By  responding 

TIME 

the  following  plot  will  be  displayed , 


YEAR 


Upon  completion  of  the  plot  the  computer  will  again  ask 


WHAT  PLOT  WOULD  YOU  LIKE.  ENTER  THE  VARIABLE,  OR  LIST,  OR  NONE 
Enter  the  next  variable  to  be  plotted 


HFSUM 


Now  the  computer  will  ask 

PLOT  THIS  VARIABLE  AGAINST  TIME  OR  ANOTHER  VARIABLE  OR  LIST 
? 

To  plot  the  total  hours  flown  against  the  number  of  pilots ,  enter 

TP 


The  following  plot  will  be  displayed 

GENERAL  AVIATION  DYNAMICS  MODEL  PAGE  3 
TOTAL  PILOTS  VS  TOTAL  HOURS  FLOWN  (THOUSANDS),  1977  TO  1980 
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Again  the  computer  will  ask 


WHAT  PLOT  WOULD  YOU  LIKE.  ENTER  THE  VARIABLE,  OR  LIST,  OR  HONE 
? 

Entering 


NONE 

will  cause  the  computer  to  print 

STEP  9  --  SAVE  THE  RESULTS  OF  THIS  FORECAST  FOR  FUTURE  SENSITIVITY  ANALYSIS. 
followed  by 

WOULD  YOU  LIKE  TO  SAVE  THE  RESULTS  OF  THIS  SESSION  FOR  LATER  SENSITIVITY 
ANALYSIS 
? 

When  performing  a  sensitivity  analysis ,  any  simulation  run  may  become  the  base¬ 
line  for  future  comparison.  If  the  current  run  is  desired  to  be  a  baseline 
for  comparisons,  enter 

YES 

The  computer  will  save  the  results  of  the  current  run  and  respond 

WOULD  YOU  LIKE  TO  CONTINUE  WITH  ANOTHER  FORECAST 
? 

The  reply 


YES 

will  cause  the  computer  to  print 

ENTER  ENDING  YEAR  FOR  SIMULATION 
? 

Entering  a  valid  year 

1980 

allows  the  system  to  continue  with 
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STEP  1  —  WOULD  YOU  LIKE  TO  COMPUTE  THE  FORECAST  WITH  THE  INITIAL  ASSUMPTIONS 
UNCHANGED  (YES  OR  NO)  OR  VIEW  THE  STEPS  OF  THIS  MODEL  (TEACH) 

? 

Since  we  are  now  familiar  with  the  steps  of  the  model  we  do  not  need  to  set  the 
TEACH  flag  on.  Answer 

NO 

since  we  already  have  stored  the  results  of  the  simulation  with  the  initial 
assumptions  unchanged.  The  computer  will  ask 

ENTER  NAME  OF  VARIABLE  TO  BE  CHANGED,  OR  LIST,  OR  NONE 
? 

The  user  has  already  decided  that  the  second  simulation  will  involve  a  fuel  tax 
increase,  but  does  not  know  how  to  implement  that  in  the  model.  Therefore 
answer 

LIST 

which  will  generate  the  following  list: 
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AIRCRAFT  VARIABLES 
IDENTIFIER  DESCRIPTION 

ADRN  AIRCRAFT  DESTRUCTION  RATE,  NORMALIZED  (AC/YR) 

FC  FIXED  COST  INDEX  (1972  $,  1972-1).  ONE  COMPONENT  OF  FC  IS  THE 

ANNUALIZED  INVESTMENT 
FCINF  FIXED  COST  INFLATION  FACTOR 

VC  VARIABLE  COST  INDEX  ($/HR) ,  (1972  $,  1972-1).  COMPONENTS  OF 

VC  ARE  FT AX  (FEDERAL  FUEL  TAX)  AND  LFEE  (LANDING  FEE). 

VCINF  VARIABLE  COST  INFLATION  FACTOR 

FHRF  FLYING  HOURS  REQUIRED  FACTOR 

ECONOMIC  VARIABLES 

IDENTIFIER  DESCRIPTION 

DPI  DISPOSABLE  PERSONAL  INCOME  (1972  $,  1972-1) 

_  GNP  GROSS  NATIONAL  PRODUCT  (1972  $,  1972-1) 

RAD  REVENUE  AIRCRAFT  DEPARTURES  (1972-1) 

ECONOMIC  DPI,  GNP,  RAD 

FUEL  VARIABLES 

IDENTIFIER  DESCRIPTION 

SFC  SPECIFIC  FUEL  CONSUMPTION 

PILOT  VARIABLES 

IDENTIFIER  DESCRIPTION 

ATrDN  AIRLINE  TRANSPORT  PILOT  DEPARTURE  RATE,  NORMALIZED 

CPDN  COMMERCIAL  PILOT  DEPARTURE  RATE,  NORMALIZED 

IPDN  INSTRUMENT  PILOT  DEPARTURE  RATE,  NORMALIZED 

PCIN  PRIVATE  CERTIFICATES  ISSUED  RATE,  NORMALIZED 

PPDN  PRIVATE  PILOTS  DEPARTURE  RATE,  NORMALIZED 

SPDN  STUDENT  PILOTS  DEPARTURE  RATE,  NORMALIZED 

URIPN  UPGRADE  TO  INSTRUMENT  PILOT  RATE,  NORMALIZED 

SCIX  STUDENT  CERTIFICATES  ISSUED  MULTIPLIER 

PILOT  ATPDN,  CPDN,  IPDN,  PCIN,  PPDN,  SPDN,  URIPN 

Upon  completion  of  the  list ,  the  computer  will  repeat  the  question 

ENTER  NAME  OF  VARIABLE  TO  BE  CHANGED,  OR  LIST  OR  NONE 
? 

The  user  can  see  from  the  above  list  that  the  federal  fuel  tax  is  a  component 
of  VC  the  variable  cost  index  and  so  enters 


VC 
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Tne  computer  responds  by  displaying  the  current  values  of  the  variable  cost 
index  components 

THE  COMPONENTS  OF  THE  VARIABLE  COST  INDEX  ARE  FT  AX,  THE  FEDERAL  FUEL  TAX,  AND 
LFEE,  THE  LANDING  FEE 

THE  CURRENT  VALUES  FOR  THE  VARIABLE  COST  INDEX  COMPONENTS  ARE 

FEDERAL  FUEL  TAX  ($/GAL) 


AV  GAS  JET 

FUEL 

1977 

* 

0.07 

0.07 

1978 

0.07 

0.07 

1979 

0.07 

0.07 

1980 

0.07 

0.07 

LANDING 

FEE  ($ /LANDING) 

SNGL-P 

SNGL-P 

MULTI¬ 

TURBO 

TURBO 

PISTON 

TURBINE 

NON-AER 

AER 

PISTON 

PROP 

JET 

HELIC 

HELIC 

1977 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1978 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1979 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1980 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

and  then  asks 


DO  YOU  WISH  TO  CHANGE  THE  VALUES  OF  FT AX,  LFEE  OR  NONE 
? 

To  change  the  fuel  tax  values  the  user  enters 

FT  AX 


and  the  computer  responds 

WHAT  YEAR  WOULD  YOU  LIKE  THE  NEW  FUEL  TAX  TO  BEGIN 
? 

For  a  January  lt  1979  date  of  effectiveness  enter 

1979 

The  computer  will  ask 

■  e 

WOULD  YOU  LIKE  THE  FUEL  TAX  TO  REMAIN  CONSTANT  FOR  ALL  SUBSEQUENT  YEARS 
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If  the  desired  fuel  tea:  change  is  to  increase  the  fuel  tax  by  64  a  year,  the 
answer  is 


NO 

The  computer  will  ask 

WOULD  YOU  LIKE  THE  FUEL  TAX  TO  CHANGE  AT  A  CONSTANT  RATE  FOR  ALL  SUBSEQUENT 
YEARS 

7 

Again  the  answer  is 


NO 

The  computer  will  respond  by  asking  for  the  fuel  tax  values  to  be  entered 
explicitly  for  each  year,  starting  with  the  year  of  the  change 

ENTER  FUEL  TAX  VALUES,  IN  DOLLARS,  FIRST  FOR  AVIATION  GAS,  THEN  FOR  JET  FUEL. 
ENTER  VALUES  FOR  EACH  YEAR. 

7 

The  1979  fuel  tax  values  are  entered  first 


0.12  0.12 


The  computer  then  asks  for  the  values  for  the  next  year 

7 


and  the  values  for  1980  are  entered 


0.17  0.17 


Since  this  simulation  is  to  end  in  1980  no  more  values  are  required.  The 
computer  then  displays  the  new  values  of  the  fuel  tax 

THE  NEW  VALUES  FOR  THE  FUEL  TAX  ARE 

FEDERAL  FUEL  TAX  ($/GAL) 


AV  GAS 

JET  FUEL 

1977 

0.07 

0;07 

1978 

0.07 

*  0.07 

1979 

0.12 

0.12 

1980 

0.17 

0.17 
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and  then  repeats  the  question 

DO  YOU  WISH  TO  CHANGE  THE  VALUES  OF  FT AX,  LFEE  OR  NONE 

7 

Since  the  new  fuel  tax  has  been  entered  as  desired  and  the  user  does  not  want 
to  impose  a  landing  fee ,  enter 


NONE 

The  computer  then  asks  whether  any  other  variables  are  to  be  changed 

ENTER  NAME  OF  VARIABLE  TO  BE  CHANGED,  OR  LIST,  OR  NONE 
? 

Since  no  other  variables  are  to  be  changed  enter 

NONE 

Since  there  are  no  more  changes ,  the  simulation  is  run.  On  finishing  execution 
of  the  simulation  the  computer  asks 

DO  YOU  WANT  TO  SEE  TABLES  OF  RESULTS  OF  THE  FORECAST 
? 


Responding 


YES 


the  computer  then  asks 

WHAT  OUTPUT  TABLE  WOULD  YOU  LIKE,  OR  ENTER  LIST  OR  NONE. 
? 

Responding  as  for  the  first  simulation , 

PILOTS 

will  generate  the  following  table 


l 
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PILOT 

DATA,  1977  TO 

1980 

1977 

1978  . 

1979 

1980 

STUDENT  PILOTS 

188,801 

183,794 

183,654 

179,129 

PRIVATE  PILOTS 

309,005 

323,821 

335,104 

345,428 

COMMERCIAL  PILOTS 

187,801 

189,699 

192,068 

194,521 

AIR  TRANSPORT  PILOTS 

45,072 

47,784 

50,879 

53,766 

PILOT  SUBTOTAL 

730,679 

745,098 

761,705 

772,844 

HELICOPTER  PILOTS 

4,804 

4,333 

3,940 

3,534 

TOTAL  PILOTS 

735,483 

749,431 

765,645 

776,378 

INSTRUMENT  RATINGS 

211,364 

221,497 

232,437 

243,148 

HELICOPTER  RATINGS 

23,012 

24,395 

25,739 

27,052 

TOTAL  HELIC  RATINGS 

27,816 

28,728 

29,679 

30,586 

Upon  completion  of  the  table  the  computer  will  again  ask 

WHAT  OUTPUT  TABLE  WOULD  YOU  LIKE,  OR  ENTER  LIST  OR  NONE. 
? 

If  no  more  output  tables  are  desired enter 


The  computer  will  ask 

DO  YOU  WANT  TO  SEE  PLOTS  OF  RESULTS  OF  THE  FORECAST 
T 

Responding 

NO 

will  cause  the  computer  to  ask 

DO  YOU  WANT  SENSITIVITY  ANALYSIS,  (THE  PREVIOUSLY  SAVED  FORECAST  IS  THE  BASELINE), 
(YES  OR  NO) 

T 
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Answer 

YES 

The  computer  will  ask 

WHAT  OUTPUT  TABLE  WOULD  YOU  LIKE,  OR  ENTER  LIST  OR  NONE 
? 

The  list  of  output  tables  for  sensitivity  analyses  is  a  subset  of  the  list  for 
absolute  forecasts.  Therefore ,  the  unfamiliar  user  should,  enter 

LIST 

which  will  generate  the  following  list  of  output  tables 

TABLE  VARIABLES  IN  TABLE 

AA1  ACTIVE  AIRCRAFT  BY  YEAR 

AA2  ACITIVE  AIRCRAFT  BY  USER  CATEGORY 

AIRPORTS  LOCAL  AND  ITINERANT  OPERATIONS  PLUS  IFR  FLIGHT  PLANS  FILED 

AIRUTIL  AIRCRAFT  UTILIZATION  RATES 

FUEL  FUEL  CONSUMED 

HOURSFLOWN  HOURS  FLOWN 

OPERATIONS  TOTAL  OPERATIONS 

PILOTS  SP.PP,CP,ATP,P,HP,TP,IP,HR,THP 

REVENUE  FEDERAL  TAX  REVENUE 

TOTALS  TOTAL  AIRCRAFT,  TOTAL  HOURS  FLOWN,  TOTAL  OPERATIONS 

The  computer  will  then  repeat  the  question 

WHAT  OUTPUT  TABLE  WOULD  YOU  LIKE,  Oil  ENTER  LIST  OR  NONE 
? 

By  responding 

PILOTS 

the  following  table  of  percent  deviations  will  be  generated 
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PILOT  DATA,  PERCENT  DEVIATION  PROM  BASELINE,  1977  TO  1980 


1977 

1978 

1979 

1980 

STUDENT  PILOTS 

0.00 

0.00 

0.00 

-2.21 

PRIVATE  PILOTS 

0.00 

0.00 

0.00 

0.24 

COMMERCIAL  PILOTS 

0.00 

0.00 

0.00 

-0.42 

AIR  TRANSPORT  PILOTS 

0.00 

0.00 

0.00 

0.00 

PILOT  SUBTOTAL 

0.00 

0.00 

0.00 

-0.52 

HELICOPTER  PILOTS 

0.00 

0.00 

0.00 

-2.05 

TOTAL  PILOTS 

0.00 

0.00 

0.00 

-0.53 

INSTRUMENT  RATINGS 

0.00 

0.00 

0.00 

-0.34 

HELICOPTER  RATINGS 

0.00 

0.00 

0.00 

0.00 

TOTAL  HELIC  RATINGS 

0.00 

0.00 

0.00 

-0.24 

Note  that  since  the  -pilot  data  reported  in  1979  is  for  the  year  1978  as  reported 
on  January  1,  1979,  the  effects  of  the  changed  fuel  tax  in  1979  are  not  felt 
until  the  1980  pilot  data.  The  computer  will  again  ask 

WHAT  OUTPUT  TABLE  WOULD  YOU  LIKE,  OR  ENTER  LIST  OR  NONE 
7 

If  no  more  sensitivity  output  tables  are  desired *  enter 

NONE 

The  computer  will  ask 

DO  YOU  WANT  TO  SEE  PLOTS  OF  THE  SENSITIVITY  ANALYSIS  RESULTS 
? 

If  sensitivity  plots  are  desired  enter 

YES 


33 


The  computer'  will  ask 

WHAT  SENSITIVITY  PLOT  WOULD  YOU  LIKE.  ENTER  THE  VARIABLE,  OR  LIST,  OR  NONE 
? 

To  obtain  the  list  of  possible  sensitivity  plots  enter 

LIST 

and  the  computer  will  generate  the  following  list 
IDENTIFIER  DESCRIPTION 


AA 

ACTIVE  AIRCRAFT  BY  PRIMARY  USE  DURING  PREVIOUS  YEAR,  AS  REPORTED 

ON  JANUARY  1  OF  DESIGNATED  YEAR,  AS  PERCENT  DEVIATION  FROM  BASELINE 

AASUM 

TOTAL  AIRCRAFT 

ATP 

AIR  TRANSPORT  PILOTS 

CP 

COMMERCIAL  PILOTS 

FC 

FUEL  CONSUMED  DURING  PREVIOUS  YEAR,  AS  REPORTED  ON  JANUARY  1  OF 

DESIGNATED  YEAR,  AS  PERCENT  DEVIATION  FROM  BASELINE 

FTR 

FEDERAL  TAX  REVENUE  DURING  PREVIOUS  YEAR,  AS  REPORTED  ON 
OF  DESIGNATED  YEAR,  AS  PERCENT  DEVIATION  FROM  BASELINE 

JANUARY  1 

HFSOM 

TOTAL  HOURS  FLOWN 

HP 

HELICOPTER  PILOTS 

HR 

HELICOPTER  RATINGS 

IP 

INSTRUMENT  RATINGS 

OPS 

OPERATIONS  (THOUSANDS)  DURING  PREVIOUS  YEAR,  AS  REPORTED 
1  OF  DESIGNATED  YEAR,  AS  PERCENT  DEVIATION  FROM  BASELINE 

ON  JANUARY 

OPSUM 

TOTAL  OPERATIONS 

PP 

PRIVATE  PILOTS 

P 

PILOT  SUBTOTAL 

SP 

STUDENT  PILOTS 

TC 

TOTAL  COST,  AS  PERCENT  DEVIATION  FROM  BASELINE 

THP 

TOTAL  HELIC  RATINGS 

TP 

TOTAL  PILOTS 

Upon 

completion  of  the  list  the  computer  will  repeat  the  question 

WHAT 

SENSITIVITY  PLOT  WOULD  YOU  LIKE.  ENTER  THE  VARIABLE,  OR  LIST, 

OR  NONE 

? 

To  see  the  effect  of  the  changed  fuel  tax  on  the  fuel  consumption,  enter 

PC 
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The  computer  will  respond 

PLEASE  ENTER  EITHER  1  FOR  AVIATION  GAS  OR  2  FOR  JET  FUEL  FOR  THE  PLOT 
? 

For  aviation  gas,  enter 

1 

and  the  computer  will  generate  the  following  sensitivity  plot 

GENERAL  AVIATION  DYNAMICS  MODEL  PAGE  6 

FUEL  CONSUMED  DURING  PREVIOUS  YEAR,  AS  REPORTED  ON 
JANUARY  1  OF  DESIGNATED  YEAR,  AS  PERCENT  DEVIATION 
FROM  BASELINE,  1977  TO  1980 

AVIATION  GAS 

*+******************+******************+******************+* 

0.00+*  *  *  + 


-0.09+ 

P  * 

E  * 

R  * 

C  -0.19+ 

E  * 

N  * 

T  * 

-0.29+ 

* 

* 

* 

-0.39+  *+ 

*+******************+******************+******************+* 

1977  1978  1979  1980 

YEAR 

Upon  completing  the  plot  the  computer  will  repeat  the  question 

WHAT  SENSITIVITY  PLOT  WOULD  YOU  LUCE.  ENTER  TjlE  VARIABLE,  OR  LIST,  OR  NONE 


For  no  further  sensitivity  plots  enter 


+ 


+ 


+ 


NONE 
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Tne  computer  will  ask 

WOULD  YOU  LIKE  TO  SAVE  THE  RESULTS  OF  THIS  SESSION  FOR  LATER  SENSITIVITY 
ANALYSIS 
? 

If  the  previous  run  is  not  to  become  a  new  baseline  enter 

NO 

Tne  computer  will  ask 

WOULD  YOU  LIKE  TO  CONTINUE  WITH  ANOTHER  FORECAST 
? 

If  no  more  forecasts  are  desired,  enter 

NO 

The  computer  will  respond 

YOU  ARE  NOW  LEAVING  THE  GENERAL  AVIATION  DYNAMICS  MODEL. 

and  return  the  user  to  the  interactive  session  with  the  computer  for  log-out. 
On  the  UCS  computer  the  computer  responds 

*END  UNICMD* 

and  the  user  can  log-out  by  entering 

BYE 

Tne  above  interactive  session  was  run  on  the  UCS  computer.  The  only  differ¬ 
ence  between  it  and  a  similar  run  on  the  Battelle  computer,  apart  from  the 
log-in  and  log-out  procedures ,  is  in  the  occurrence  of  a  /  as  a  user  prompt 
instead  of  ?  as  above. 
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Section  IV:  Batch  use  of  the  GAD  model 


The  GAD  model  may  be  run  in  batch  mode  on  both  the  Battelle 
Computer  and  the  UCS  computer 

IV. 1  Batch  mode  on  the  Battelle  computer 

The  control  cards  required  are  as  follows: 

Job  card 

ATTACH , PROFIL , AVIATION , ID-NUCLEUS ,MR-1. 

BEGIN, GAD. 

*EOR 

SYSTEM, USER-OFF, WIDTH-131 
LOAD, GAD 

Appropriate  responses  to  GAD  conversation,  one  per  card. 


*EOF 


IV. 2  Batch  mode  on  the  UCS  computer 

The  control  cards  required  are  as  follows : 

Job  card 
Account  card 

Cards  giving  Instructions  for  mailing  output,  if  any. 
GET, NUCLEUS. 

GET.TAPEl-GADCOMP .TAPE2-GAFILE2 ,TAPE3-GAFILE3. 

NUCLEUS . 

EOR 

SYSTEM , USER-OFF , WIDTH- 131 
LOAD, GAD 

Appropriate  responses  to  GAD  conversation,  one  per  card. 


EOF 


